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ABSTRACT 

Electro-optic  (EO)  devices  that  modulate  optical  signals  by  electric  fields  are  an 
integrative  part  of  the  photonics  industry  and  device  optimization  is  an  important  area  of 
research.  As  applications  move  to  large  bandwidth  and  higher  frequency,  low  electro-optic 
effects  and  the  requirement  for  large  dimension  [2]  become  restrictive  for  microwave-optical 
devices.  Both  experimental  and  computational  evaluations  indicate  that  strain  and  polarization 
distribution  have  a  significant  impact  on  electromagnetic  wave  propagation  resulting  from  a 
resonant  structure;  however,  no  systematic  study  or  fundamental  understandings  are  available. 

This  dissertation  research  has  been  carried  out  to  study  and  further  develop  the  subject  of 
piezoelectric  resonance  enhanced  electro-acoustic-optic  process,  in  order  to  improve  the 
sensitivity  and  efficiency  of  electro-optic  sensors  and  to  explore  novel  applications.  Many  finite 
element  models  have  been  constructed  for  evaluating  the  mechanisms  of  the  phenomena  and  the 
effectiveness  of  the  device  structure.  The  enhancement  in  transmission  is  found  to  be  directly 
related  to  the  strain-coupled  local  polarization.  At  piezoelectric  resonance  oscillating  dipoles  or 
local  polarizations  become  periodic  in  the  material  and  have  the  greatest  impact  on  transmission. 

Results  suggest  that  the  induced  charge  distribution  by  a  piezoelectric  material  at  certain 
resonant  frequencies  is  effective  for  aiding  or  impeding  the  transmission  of  a  propagating  wave. 

The  behavior  of  both  piezoelectric-defined  (or  intrinsic  piezoelectric  materials)  and  engineered 

periodic  structures  are  reported.  The  piezoelectric  response  of  the  surface  displacement  of  samples  is  investigated 

using  an  ultra-high  frequency  laser  Doppler  vibrometer.  A  two 

dimensional  view  of  the  surface  is  obtained  and  the  surface  displacement,  velocity  and 
acceleration  are  compared  to  the  electro-optic  response  under  the  resonant  condition.  A  study  of 
the  acousto-optic  (AO)  effect  in  a  family  of  oxide  crystals  (including  e.g.,  Ti02,  ZnO,  LiNb03, 
and  ferroelectric  perovskites)  has  been  conducted  by  the  finite  element  analysis  method.  This 
study  further  serves  to  show  the  potential  of  optimizing  devices  through  a  consideration  of  their 
directional  dependent  parameters  and  resonant  behavior.  The  acousto-optic  figure  of  merit 
(FOM)  as  a  function  of  the  material's  refractive  index,  density,  effective  AO  coefficient  and  the 
velocity  of  the  acoustic  wave  in  the  material,  is  also  investigated.  By  examining  the  directional 
dependent  velocity,  acousto-optic  coefficients,  and  refractive  index,  the  acousto-optic  FOM  can 
be  calculated  and  plotted  in  all  directions  revealing  the  optimal  crystal  orientation  to  maximize 
coupling  between  the  optical  and  acoustic  waves.  A  finite  element  model  was  developed  to 
corroborate  the  improved  interaction.  The  model  examines  the  diffraction  that  occurs  on  the 
optical  wave  as  it  travels  through  an  acousto-optic  medium.  The  combined  information  gained 
from  commercially  available  multiphysics-based  modeling  platforms  is  shown  to  be  an  effective 
means  of  predicating  acousto-optic  device  functionality. 
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Electro-optic  (EO)  devices  that  modulate  optical  signals  by  electric  fields  are  an 
integrative  part  of  the  photonics  industry  and  device  optimization  is  an  important  area  of 
research.  As  applications  move  to  large  bandwidth  and  higher  frequency,  low  electro-optic 
effects  and  the  requirement  for  large  dimension  [2]  become  restrictive  for  microwave-optical 
devices.  Both  experimental  and  computational  evaluations  indicate  that  strain  and  polarization 
distribution  have  a  significant  impact  on  electromagnetic  wave  propagation  resulting  from  a 
resonant  structure;  however,  no  systematic  study  or  fundamental  understandings  are  available. 

This  dissertation  research  has  been  carried  out  to  study  and  further  develop  the  subject  of 
piezoelectric  resonance  enhanced  electro-acoustic-optic  process,  in  order  to  improve  the 
sensitivity  and  efficiency  of  electro-optic  sensors  and  to  explore  novel  applications.  Many  finite 
element  models  have  been  constructed  for  evaluating  the  mechanisms  of  the  phenomena  and  the 
effectiveness  of  the  device  structure.  The  enhancement  in  transmission  is  found  to  be  directly 
related  to  the  strain-coupled  local  polarization.  At  piezoelectric  resonance  oscillating  dipoles  or 
local  polarizations  become  periodic  in  the  material  and  have  the  greatest  impact  on  transmission. 
Results  suggest  that  the  induced  charge  distribution  by  a  piezoelectric  material  at  certain 
resonant  frequencies  is  effective  for  aiding  or  impeding  the  transmission  of  a  propagating  wave. 
The  behavior  of  both  piezoelectric-defined  (or  intrinsic  piezoelectric  materials)  and  engineered 
periodic  structures  are  reported.  The  piezoelectric  response  of  the  surface  displacement  of 
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samples  is  investigated  using  an  ultra-high  frequency  laser  Doppler  vibrometer.  A  two 
dimensional  view  of  the  surface  is  obtained  and  the  surface  displacement,  velocity  and 
acceleration  are  compared  to  the  electro-optic  response  under  the  resonant  condition.  A  study  of 
the  acousto-optic  (AO)  effect  in  a  family  of  oxide  crystals  (including  e.g.,  TiCF,  ZnO,  LiNb03, 
and  ferroelectric  perovskites)  has  been  conducted  by  the  finite  element  analysis  method.  This 
study  further  serves  to  show  the  potential  of  optimizing  devices  through  a  consideration  of  their 
directional  dependent  parameters  and  resonant  behavior.  The  acousto-optic  figure  of  merit 
(FOM)  as  a  function  of  the  material's  refractive  index,  density,  effective  AO  coefficient  and  the 
velocity  of  the  acoustic  wave  in  the  material,  is  also  investigated.  By  examining  the  directional 
dependent  velocity,  acousto-optic  coefficients,  and  refractive  index,  the  acousto-optic  FOM  can 
be  calculated  and  plotted  in  all  directions  revealing  the  optimal  crystal  orientation  to  maximize 
coupling  between  the  optical  and  acoustic  waves.  A  finite  element  model  was  developed  to 
corroborate  the  improved  interaction.  The  model  examines  the  diffraction  that  occurs  on  the 
optical  wave  as  it  travels  through  an  acousto-optic  medium.  The  combined  infonnation  gained 
from  commercially  available  multiphysics-based  modeling  platforms  is  shown  to  be  an  effective 
means  of  predicating  acousto-optic  device  functionality. 
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CHAPTER  1:  INTRODUCTION 


1.1  Significance  of  the  Subject  of  Study 

This  dissertation  focuses  on  experimental  detennination  of  methods  and  systems  for  the 
enhancement  of  optical  transmission  by  way  of  periodic  structures.  Both  numerical  simulation 
and  theoretical  explanation  are  incorporated  to  better  guide  the  experimental  setup  and  give 
improved  understanding  to  resultant  data.  The  outcome  of  this  research  could  improve  the 
overall  efficiency  of  electro-optic  (EO)  devices.  Higher  EO  coefficient  at  RF  and  microwave 
frequencies  can  lead  to  much  better  performance  with  reduced  device  size.  Greater  modulation 
depth  also  results  in  electro-optic  modulators  (EOM)  with  larger  deflection  angles.  High  EO 
coefficient  means  EOM’s  would  require  less  power  to  perfonn  the  same  function,  more  compact 
and  more  linear  operation.  While  EOMs  do  not  require  much  power  individually,  it  is  desirable, 
if  not  imperative,  to  reduce  the  energy  consumption  of  devices  for  large-scale  implementation. 
Additionally  an  emerging  technology  in  the  computer  world  is  new  processors  which  use  optics; 
power  efficiency  is  a  great  concern  here  especially  when  heat  dissipation  of  new  processors  is  a 
limiting  factor.  The  knowledge  gained  through  this  dissertation  research  will  help  detennine  best 
design  practices  for  better  coupling  of  the  EO  effect  and  electromagnetic  (EM)  transmission. 
Additionally  a  better  understanding  will  be  gained  on  why  resonance  improves  transmission  and 
coupling,  when  it  will  occur,  and  how  to  control  it.  Finally  this  research  will  contribute  to  our 
understanding  of  effective  experimental  and  modeling  practices  underlying  multiphysics 
numerical  methods. 

There  has  been  much  work  on  acousto-optics  and  electro-optics  but  little  on  acousto- 

electro-optics  [1].  There  are  few  papers  [1,  3,  4]  addressing  the  enhanced  EO  coupling  that 

occurs  at  resonance.  There  are  very  few  examples  of  simulation/modeling  of  spontaneous 
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polarization  in  materials  and,  even  more  sparse,  are  examples  of  complex  muliphysics  models 
dealing  with  interactions  between  acoustics,  piezoelectric,  RF,  EM  wave  and  induced  periodic 
polarization  changes.  An  effort  to  model  multi-physics  interactions  such  as  this  might  be 
considered  in  itself  a  significant  research  topic;  the  powerful  combination  of  multiphysics 
modeling  and  advanced  experiments  has  yet  to  be  explored.  There  has  been  much  work  on 
imaging  the  surface  of  materials  using  techniques  like  profilometry  and  atomic  force  microscopy 
and  some  on  imaging  vibrating  materials  using  vibrometry,  but  little  work  has  been  done 
imaging  the  surface  while  monitoring  the  optical  transmission. 

The  interaction  between  light  and  piezoelectric  materials  is  very  important  for  optical 
communications.  All  transparent  piezoelectric  materials  are  also  electro-optic  and  photoelastic. 
Photoelasticity  can  be  expressed  in  terms  of  piezo-optic  or  elasto-optic  effects.  It  has  long  been 
observed  that  the  electro-optic  output  is  strongly  influenced  by  low  frequency  (kHz)  vibrations 
due  to  natural  piezoelectric  resonate  frequencies  [5],  Typically  electro-optic  modulators  are  not 
operated  in  the  region  of  resonance  vibration  to  avoid  non-monotonic  responses;  however  it  has 
been  shown  that  operation  at  these  frequencies  and  at  harmonics  can  actually  be  advantageous  to 
improve  electro-optic  response  [4].  The  response  of  the  electro-optical  signal  tends  to  follow  the 
slope  of  the  magnitude  of  the  electrical  admittance  as  a  function  of  frequency,  this  is  fairly 
unsurprising  as  the  largest  mechanical  vibrations  occur  at  resonance  (typically  between  the 
resonant  and  anti-resonant  frequencies)  and  in  turn  cause  a  large  change  in  refractive  index. 
Further  results  by  Guo  et  al.  [1]  have  shown  that  certain  resonance  modes  with  vibrations  mostly 
in  the  direction  transverse  to  optical  propagation  have  a  larger  influence  on  response  because 
they  are  parallel  to  the  optical  polarization.  Given  that  the  spurious  response  of  these  crystals 
little  attention  has  been  made  to  the  detailed  study  of  their  response  in  conjunction  with  the 
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optical  field  around  resonance,  leaving  it  an  open  and  promising  area  of  research  for  further 
optimization  of  electro-optic  devices. 

1.2  Scope  of  Dissertation 

The  goal  of  this  dissertation  research  was  to  study  and  further  develop  the  subject  of 
piezoelectric  resonance  enhanced  electro-acoustic-optic  processes  in  order  to  improve  the 
sensitivity  and  efficiency  of  electro-optic  sensors  and  to  explore  novel  applications.  Both 
experimental  study  and  numerical  simulations  are  essential  approaches  to  accomplish  the  goal. 

1.3  Organization  of  the  Dissertation 

Following  a  comprehensive  but  selective  review  on  the  background  of  the  research 
summarized  in  Chapter  2,  electro-optic  and  electrical  experiments  performed  to  confirm  the 
enhancement  reported  by  previous  researchers  and  to  investigate  its  limits  are  described  in 
Chapter  3.  The  standard  method  of  crossed  polarizers  with  a  quarter  wave  plate  for  biasing  (the 
Senannont  compensator  technique)  was  used  to  measure  the  linear  electro-optic  Pockels  effect  in 
the  test  samples.  The  results  confirmed  a  very  large  spike  in  the  response  close  to  the  resonance 
that  loosely  correlates  with  the  slope  of  the  Admittance  as  discussed  in  [4], 

To  better  understand  this  effect,  how  to  control  it,  when  it  occurs,  how  to  maximize  it  and 
to  extend  our  ability  to  quickly  test  under  different  conditions,  several  models  were  constructed 
using  Finite  Element  Analysis.  The  details  of  these  models  are  discussed  in  Chapter  4.  Section 
4.2  describes  in  detail  a  model  of  the  transmission  of  a  microwave  EM  signal  through  a 
ferroelectric,  periodically-poled  sample.  Several  poling  configurations  are  used  to  adjust  the 
amplitude  of  the  polarization  as  well  as  the  number  and  direction  of  the  domains.  It  is  observed 
that  there  is  a  polarization  threshold  of  about  1  x  10-12  C/cm2  needed  to  affect  any  change  in 
transmission  (in  present  simulation)  and  after  that  the  polarization  amplitude  can  have  a 
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significant  effect  on  transmission.  The  number  of  polarization  elements  can  be  adjusted  to  tune 
the  transmission  magnitude  provided  that  the  polarization  directions  are  transverse  to  wave 
propagation.  Longitudinal  polarization  in  the  same  direction  as  the  propagating  wave  does  not 
show  any  significant  influence  on  transmission  which  corroborates  with  [4].  A  separate  model 
was  also  developed  to  mimic  the  behavior  of  our  experimental  electro-optic  tests.  This  model 
contains  two  sub  models;  the  piezoelectric  one  and  the  optical  one,  both  were  examined  with  a 
time-harmonic  (frequency  domain)  approach.  The  piezoelectric  model  is  used  to  examine  the 
resonant  behaviors  of  a  piezoelectric  sample  and  the  results  from  this  model  are  coupled  into  the 
optical  model.  The  optical  model  is  used  to  test  the  amount  of  signal  transmitted  through  the 
device  during  resonance.  This  model  has  been  shown  to  accurately  predict  the  increased 
transmission  possible  near  resonance  that  is  consistent  with  experimental  methods.  Additionally 
the  model  used  a  special  coupling  between  the  3D  piezoelectric  and  2D  optical  portions  to 
reduce  the  required  computational  power  and  made  possible  the  simultaneous  evaluation  of  a 
millimeter  sized  device  with  a  propagating  micrometer  EM  wave.  This  kind  of  multiphysics 
study  is  nearly  impossible  without  the  method  derived  and  is  rarely  attempted.  This  model  is  also 
a  useful  tool  for  portraying  the  displacive  current  distribution  near  resonance.  As  proposed  in  [4] 
displacive  current  fonns  a  complete  circuit  near  resonance  and  leads  to  enhancement  of  the 
optical  signal  -  the  model  has  served  to  confirm  this  theory.  In  section  4.4  another  model 
discussed  and  used  to  examine  the  experimental  technique  used  by  Johnson  et  al.  [4,  6]  using  a 
dual  signal  approach  with  ac  biasing  near  piezo-resonance.  This  model  was  simulated  in  the  time 
domain  rather  than  the  frequency  domain  model  discussed  previously;  this  is  because  the  time- 
harmonic  solver  only  allows  for  a  single  excitation  frequency.  Using  the  time  domain  solver  two 
signals  can  easily  be  added.  This  model  does  not  examine  the  optical  response  but  instead 
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calculates  the  piezoelectric  response  and  the  effect  on  the  refractive  index.  The  index  change  is 
appropriate  for  this  study  because  it  is  proportional  to  the  induced  phase  change  in  the  optical 
signal.  This  model  accurately  shows  a  wide  bandwidth  enhancement  of  the  index  change  when 
the  ac  biasing  is  set  to  a  frequency  near  the  material’s  piezo-resonant  frequency.  There  have 
been  many  qualitative  statements  about  the  relative  differences  in  half  wave  voltage  that  can  be 
achieved  with  resonance  and  non-resonant  electro-optic  modulators,  yet  few  examples  offered  a 
quantitative  comparison.  A  finite  element  model  is  discussed  in  section  4.5  that  compares  a 
resonator’s  half  wave  voltage  over  a  range  of  visible  wavelengths.  This  type  of  comparison  is 
used  as  it  is  more  familiar  to  optical  researchers.  It  was  demonstrated  that  piezo-resonant 
modulators  are  significantly  more  efficient  than  even  resonators  assembled  using  an  electrically 
resonant  circuit.  Furthermore  piezo-resonators,  as  predicted  by  the  model  offer  a  half  wave 
voltage  that  can  be  as  much  as  1000  times  less  than  that  of  non-resonant  modulators. 

To  better  understand  the  transversal  lattice  strain  associated  with  piezoelectric  resonance 
modes  that  enhance  the  propagating  optical  signal  several  samples  were  studied  in  Chapter  5 
with  an  ultra-high  frequency  Laser  Doppler  Vibrometer.  This  system  enables  high  resolution 
measurements  of  the  displacement,  velocity,  and  acceleration  of  the  vibrating  sample  at  a  given 
point.  An  array  of  tests  were  conducted  over  the  surface  of  the  sample  to  build  up  a  full  view  of 
the  modes  on  a  given  surface.  Comparing  these  results  to  the  optical  response  indicates  that  even 
in  the  megahertz  regime  where  the  sample  becomes  constrained  and  strain  is  minimal  there  are 
still  spikes  in  the  electro-optic  response  which  is  further  evidence  that  enhancement  is  driven  by 
more  than  high  strains  but  perhaps  a  result  of  the  periodically  induced  polarization  that  can  still 
respond  at  those  frequencies. 
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Chapter  6  extends  the  study  to  acousto-optic  optimization  in  detail.  This  dissertation  is 
about  more  than  just  piezo-resonance,  it  is  in  a  more  general  sense  about  improved  sensitivity  of 
electro-optic  devices.  While  electro-optic  effect  is  predominant  in  the  materials  studied 
vibrations  of  lattices  near  resonance  result  in  high  strains  that  in  turn  affect  the  birefringence  of  a 
crystal  through  the  acousto-optic  effect.  The  best  orientation  for  applying  a  field  for  a  large  index 
change  was  evaluated  using  the  acousto-optic  figure  of  merit  resulting  from  a  full  set  of  data 
including  the  density,  refractive  index,  piezoelectric  constants,  and  elasticity  of  the  material. 
Several  common  materials  were  evaluated  and  it  was  observed  that  the  best  orientation  is  not 
always  along  the  principal  axes. 

The  work  discussed  in  and  the  knowledge  gained  through  this  dissertation  has  a  broad 
range  of  potential  applications  and  opened  up  possibilities  for  many  new  areas  of  research. 
Chapter  7  discusses  proposed  future  research  and  the  attempt  to  design  a  few  possible 
experiments  that  may  be  used  to  investigate  new  applications  for  piezoelectric  resonance 
enhancement.  The  main  interest  described  is  for  an  acoustic  metamaterial  for  RF  and  microwave 
“cloaking”.  Due  to  the  fact  that  resonance  can  influence  the  EM  transmission  in  a  material  it  may 
be  possible  to  make  and  opaque  object,  transparent  (or  vice  versa).  Because  resonance  can  be 
excited,  or  turned  on/off  with  the  applied  signal  it  may  be  uniquely  suited  for  making  a 
dynamically  adjustable  “cloak”.  This  chapter  discusses  a  proposed  research  path  to  realize  such  a 
device.  An  Overall  conclusion  is  given  in  Chapter  8  highlighting  the  main  scientific  and 
technical  advancements  this  dissertation  has  been  able  to  achieve.  Additionally  the  Appendix 
contains  a  great  deal  of  information  on  using  Finite  Element  Analysis  to  calibrate  a  microwave 
waveguide;  this  technique  can  be  useful  for  electro-optic  evaluation  and  frequencies  where 
lumped  methods  are  no  longer  valid. 
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CHAPTER  2:  THEORY  AND  BACKGROUND 


One  of  the  earliest  examples  of  a  significant  change  in  the  electro-optic  coefficient  was 
reported  by  Pisarevski  et  al.  [5].  The  two  samples  used  in  the  experiment  showed  a  drop  in 
permittivity  around  10  kHz.  This  drop  in  permittivity  corresponds  to  a  large  increase  and  sudden 
decrease  in  the  electro-optic  coefficient  in  both  crystals  measured  (See  Figure  2.1). 


Figure  2.1  Permittivity  and  electro-optic  response  of  crystals  near  resonance. 
Figure  from  [5] 


As  examined  by  Johnson  et  al.  [4,  6]  the  EO  effect  can  be  enhanced  not  only  at  resonance 
but  continuously  between  resonance  and  several  harmonics  through  a  method  of  frequency 
mixing.  Experiments  were  carried  out  at  frequencies  below  the  low  MHz  region  using  the 
Senarmont  compensator  technique.  Two  signals  simultaneously  applied  to  the  sample,  one  at  the 
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fundamental  resonance  mode  of  the  sample  and  the  other  which  is  swept  from  well  below 
resonance  to  well  above  it.  The  results  which  are  given  in  terms  of  optical  transmission  (voltage 
monitored  at  the  photodetector)  show  that  certain  resonance  modes  have  enhanced  optical 
transmission  and  furthermore  that  transmission  is  enhanced  over  a  range  of  frequencies  out  to 
even  the  5th  hannonic  of  resonance.  It  is  cited  that  though  enhancement  decreases  after  this  point, 
it  is  likely  due  to  the  frequency  limitations  of  the  measurement  devices;  specifically,  the  lock-in 
amplifier. 

It  is  not  surprising  that  Johnson’s  experiment  gave  more  optical  transmission  output  as 
two  signals  used  simultaneously,  and  presumably  these  signals  used  in  phase  such  that,  when 
equal,  there  was  constructive  interference  and  then,  as  the  frequencies  parted,  there  would  be 
periodic  constructive  and  destructive  interference.  But  this  interestingly  was  not  what  was 
measured  by  Johnson  at  the  photodetector.  In  experiments  that  involve  one  signal  of  constant 
amplitude  and  changing  frequency  the  electro-optic  coefficient  is  certainly  proportional  to  the 
optical  transmission  at  the  photodetector.  However,  due  to  the  frequency  mixing  technique  used 
by  Johnson,  this  is  no  longer  the  case.  It  would  be  more  appropriate  to  monitor  both  the  rms 
signal  at  the  sample  and  at  the  photodetector  and  take  the  ratio  of  the  two  which  will  then  be 
proportional  to  the  EO  coefficient.  It  is  likely  that  the  results  will  still  show  an  improved  EO 
coupling  at  the  piezoresonance  harmonics,  but  perhaps  not  as  profound  an  enhancement  as 
previously  reported. 

The  most  recent  study  on  enhanced  piezoresonance  optical  transmission  by  Guo  et.  al.  [1] 
was  done  on  PZN-8PT  crystal.  This  study  was  done  on  the  fundamental  resonant  modes  d33,  d3i, 
d32  and  it  was  found  that  the  optical  transmission  increased  near  the  d33  and  d3i  modes  where  the 
strain  is  transverse  to  the  propagation  of  light  (See  Figure  2.2).  The  d32  mode  did  not  produce 
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transmission  as  its  strain  was  parallel  to  optical  transmission.  The  authors  cite  that  the  effective 
EO  coefficient  is  largest  when  there  is  a  large  change  in  electrical  admittance  with  frequency. 


U//O0  « 


dT(nf) 

da) 


(2.1) 


Figure  2.2  Comparison  of  experimental  results  of  optical  signal,  phase  angle 
and  admittance  for  PZN-PT  sample  under  electro-optic  modulation 


They  also  note  that  there  is  no  apparent  correlation  between  the  admittance  magnitudes  with  the 
enhanced  optical  transmission.  Using  this  information  another  experiment  was  performed  at  the 
frequency  corresponding  to  the  d33  mode  varying  the  amplitude  of  the  applied  signal  and 
observing  the  optical  transmission.  Results  showed  that  the  signal  strength  at  the  photo  diode  (in 
volts)  was  linearly  proportional  to  the  square  of  the  applied  voltage.  This  is  significant  because 
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under  DC  conditions  the  optical  signal  corresponds  linearly  with  the  applied  voltage.  The  authors 
give  the  refractive  index  change  of  the  crystal  for  each  fundamental  resonance  mode  as: 

=  A  (—)  =  (Pll33  +  P3333)X33  (2-2) 

\n  /1:L  \n  /33 

^31:^(“)  =^(— 2)  =  (Pllll  +  P331l)Xll  (2-3) 

Vn//11  Vn//33 

1  a  1  a 

d32  ■  ^  I  2  )  =^(_ 2/  =  (Pll22  +  P331l)X22  (2-4) 

The  piezoelectric  d  coefficients  indicate  a  coupling  between  the  applied  electric  field  in  the  3 
direction  and  induced  strain  in  the  1,  2  or  3  directions  respectively  for  d3i,  d32,  d33.  The 
photoelastic  p  coefficients  show  the  coupling  between  the  strain  and  refractive  index.  It  was  also 
suggested  that  enhancement  of  d3J  may  be  larger  because  pun  and  p33n  are  of  the  same  sign  and 
thus  additive  corresponding  to  a  larger  index  change,  while  P1122  and  p33n  are  of  a  different  sign 
leading  to  a  smaller  index  change  for  d32. 

In  an  experiment  similar  to  that  of  Shikik  Johnson,  the  Naval  Research  Laboratory  found 
that  the  electro-optic  effect  could  be  enhanced  in  a  Lithium  niobate  crystal  through 
piezoresonance  [3].  The  Senarmont  compensator  technique  was  used  to  measure  the  EO 
coefficient  of  the  crystal  at  frequencies  up  to  700  KHz.  At  the  piezoresonance  of  the  sample  the 
photo  elastic  effect  dominates.  The  sample  was  tested  under  two  different  field  configurations. 
First,  a  sinusoidal  signal  was  used  to  show  that  the  EO  effect  is  enhanced  at  the  resonant  modes. 
It  was  also  shown  that  the  enhancement  can  be  applied  at  frequencies  far  higher  than  just  the 
resonant  frequencies.  Rather  than  a  sinusoidal  input  at  InS,  a  long  pulse  was  applied  with  a 
repetition  rate  equal  to  the  resonant  frequency.  The  authors  provide  the  following  explanation  for 
their  results:  “Our  modeling  of  the  data  suggests  that  the  sensitivity  enhancements  are  produced 
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by  the  interplay  between  photo  elastic  shifts  in  the  refractive  indices  and  the  physical  vibration 
modes  of  the  crystals.”  The  research  group  that  produced  this  paper  has  made  some  investigation 
on  dielectric  enhancing  of  EO  [7]  but  has  not  yet  produced  any  follow  up  papers  specifically  on 
resonance  enhancement. 
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CHAPTER  3:  PIEZORESONANCE  EXPERIMENTS 

3.1  Piezoelectricity 

Piezoelectricity  is  a  property  found  in  certain  materials  that  results  in  an  electrical  charge 
as  a  result  of  mechanical  stress.  It  is  widely  used  in  in  many  devices.  Piezoelectricity  only  occurs 
in  non-centric  materials.  When  a  stress  is  applied  to  a  piezoelectric  material  the  atoms  arrange  so 
there  is  a  net  charge.  Piezoelectric  resonators  are  used  for  timing  pieces  in  digital  clocks. 
Piezoelectric  motors  can  be  used  for  optical  zoom  in  cameras.  Piezoelectric  materials  are  used 
for  detecting  underwater  sound  vibrations  in  submarines  and  it  is  also  a  common  method  of 
producing  a  spark  in  lighters. 

There  are  two  equivalent  types  of  piezoelectricity,  the  direct  effect  that  couples  stress  Xkt 
with  electric  polarization  Pj  and  the  converse  effect  that  the  electric  field  Ek  causes  a  change  in 
strain  xi;-.  The  coefficient  for  the  direct  fonn  is  [C/N]  and  for  the  converse  form  [m/V], 

pj  =  djkixki  Direct  (3.1) 

xij  =  dijkEk  Converse  (3.2) 

The  equations  above  demonstrate  the  tensor  form  of  the  basic  piezoelectric  coupling  relations.  A 
more  detailed  view  of  factors  involved  can  be  seen  in  the  piezoelectric  constitutive  equations  (in 
matrix  fonn)  for  which  there  are  two  equivalent  types;  the  Stress-Charge  and  Strain-Charge 
fonns. 

The  Stress-Charge  defines  the  total  stress  X  as  a  function  of  the  elastic  stiffness  cE,  the 
induced  and  initial  strain  (x-x;),  the  transpose  of  the  piezoelectric  coupling  matrix  e1 ,  the 
electric  field  E,  and  the  initial  stress  Xt .  The  electric  displacement  D  is  a  function  of  the 
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piezoelectric  coupling  e,  the  induced  and  initial  strain  (x-x.) ,  the  free  space  permittivity  e0 ,  the 
relative  permittivity  srS  ,  and  the  remnant  electrical  displacement  Dr. 

Stress-Charge  form: 

X  -  cE{x  - xt)-eT E  +  Xi  (3.3) 

D  =  e(x-xi)-s0srSE+Dr  (3.4) 

The  Strain-Charge  defines  the  Total  strain  x  as  a  function  of  the  elastic  compliance  sE, 
the  induced  and  initial  stress  ( X  -  X; ) ,  the  transpose  of  the  piezoelectric  coupling  matrix  d 1 , 
the  electric  field  E,  and  the  initial  strain  x. .  The  electric  displacement  D  is  a  function  of  the 
piezoelectric  coupling  d,  the  induced  and  initial  stress  (X  -  X( ) ,  the  free  space  permittivity  £0 , 

the  relative  pennittivity  at  constant  temperature  srT ,  and  the  remnant  electrical  displacement. 
Srain-Charge  form: 

x-sE{X-Xi)-dTE  +  xi  (3.5) 

D  =  d(X  —  X;)  —  s{)srJE  +  Dr  (3.6) 

Piezoelectricity  is  a  third  rank  tensor  property,  meaning  that  it  can  only  occur  in  non- 
centrosymmetric  materials  (materials  with  no  center  of  symmetry).  There  are  32  crystal  classes, 
1 1  possess  a  center  of  symmetry  and  thus  are  not  piezoelectric.  One  class,  point  group  432, 
though  it  is  non  centrosymmtric,  still  does  not  exhibit  the  piezoelectric  effect  because  the  charge 
along  the  <1 1 1>  direction  cancels  itself  out  [8].  This  leaves  20  point  groups  with  piezoelectricity 
which  are  listed  below.  The  polar  point  groups  may  have  a  spontaneous  polarization  that  occurs 
even  without  an  applied  electric  field  and  are  also  pyroelectric.  There  are  also  non- 
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crystallographic  point  groups,  called  Curie  groups  that  can  also  be  piezoelectric  but  are  not  listed 
here. 

Table  3.1  List  of  crystal  point  groups  that  allow  for  piezoelectricity 


Polar 

1, 2, 3, 4, 6,  m,  mm2, 3m,  4 mm,  6 mm 

Non-polar 

222, 4, 422, 42 m,  32,  6,  6m2, 622, 43 m,  23 

Piezoelectric  materials  experience  a  resonant  state  when  an  electrical  signal  is  applied  at  certain 
frequencies  as  determined  by  the  density  and  elastic  parameters  of  the  material  and  the 
dimension  and  geometry  of  the  sample.  The  resonant  frequency  (fr)  is  the  point  where  the 
admittance  is  maximum  and  the  anti-resonance  (fa)  is  where  the  admittance  is  at  a  minimum  (See 
Figure  3.1). 
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Log  Frequency  (Hz) 


Figure  3.1  Admittance  frequency  spectrum  of  a  typical  piezoelectric  resonator 
showing  resonance  and  anti-resonant  frequencies 

At  fr  the  admittance  is  at  maximum  (thus  an  impedance  minimum).  If  the  voltage  is  kept  constant 
(constant  amplitude  sinusoidal  signal),  there  will  be  an  influx  of  current  into  the  device.  The 
constant  shape  change  of  the  device  causes  motional  capacitance. 

The  resonant  frequency  fr  is  often  denoted  as  fm  referring  to  the  frequency  on  minimum 
impedance  and  the  anti-resonant  frequency  fa  denoted  as  fn  is  the  frequency  of  maximum 
impedance.  The  electromechanical  coupling  factor  k  is  a  measure  of  how  well  a  piezoelectric 
element  can  convert  electrical  energy  into  mechanical  energy.  There  are  several  types  of 
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coupling  factors  that  are  dependent  on  the  geometry  and  application  direction  of  the  electric 
field.  The  effective  coupling  coefficient  which  applies  to  any  geometry  can  be  calculated  from 
Equation  (3.9).  The  coefficient  k31  is  applied  to  plates  where  the  induced  strain  is  perpendicular 
to  the  polarization.  Conversely  the  coefficient  k33  is  used  for  rods  where  the  induced  strain  is 
parallel  to  the  polarization.  For  both  the  k31  and  k33  coupling  the  applied  electric  field  must  be 
in  the  same  direction  as  the  polarization. 

(Ke„f  =  (3-7) 
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Resonance  in  a  piezoelectric  vibrator  occurs  with  useful  effects  in  the  kilohertz  region  for 
most  materials  and  sample  configurations.  As  the  frequencies  of  agitation  is  increased  the  sample 
cannot  respond  fast  enough  to  the  applied  signal  to  keep  up;  the  result  of  this  is  that  the  crystal 
becomes  more  constrained  as  frequency  increases,  this  is  the  constant  strain  state.  There  are  two 
main  options  around  this  issue,  one  is  use  very  small  samples  that  may  respond  faster  and  the 
other  is  to  induce  vibrations  on  the  air-dielectric  interface;  the  latter  is  the  bases  of  surface 
acoustic  waves  (SAW).  SAW  devices  typically  use  interdigitized  coplanar  electrodes  to  launch 
acoustic  wave  along  the  surface  of  a  piezoelectric  material.  The  wave  stays  along  the  surface  and 
only  a  few  wavelengths  actually  propagate  into  the  material.  SAW  devices  often  are  operated  at 
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hundreds  of  megahertz;  they  are  widely  used  as  electrical  filters  for  cell  phones,  satellites 
communicators  and  other  RF  devices. 


3.2  Crystal  Optics 

There  are  three  main  contributions  to  refractive  index  changes  relevant  to  this  thesis,  the 
photoelastic  effect,  optical  Pockels  and  Kerr  effects.  In  each  of  these  effects  the  refractive  index 
change  is  often  expressed  in  terms  of  the  impermeability  tensor  B,  defined  as: 


1 


(3.10) 


This  impermeability  is  the  difference  between  the  inverse  of  the  square  of  the  natural  material 
index  and  the  inverse  of  the  square  of  the  induced  index.  The  fourth  rank  tensor  photoelastic 
effect  has  two  forms,  piezo-optic  defined  in  terms  of  stress  (X)  and  the  elasto-optic  defined  in 
tenns  of  strain  (x).  We  can  convert  piezo-optic  to  elasto-optic  by  the  matrix  formula 
Pnm  =  71  mPc pn  making  use  of  the  elastic  coupling  c. 


ABij  =  nijklXkl  piezo-optic  (3.11) 

ABij  =  Pijkixki  elasto-optic  (3.12) 

The  electro-optic  effect  comes  in  two  forms,  the  linear  Pockels  effect  (3rd  rank  tensor  r^-fe)  and 
the  quadratic  Kerr  effect  (4th  rank  tensor  Rijki).  Both  effects  can  occur  in  ferroelectric  crystals; 
however  the  Pockels  effect  tends  to  be  far  more  dominant  for  low  field  conditions. 


ABij  =  rijkEk  Pockels 


ABij  =  BijkiEkEi  Kerr 


(3.13) 

(3.14) 


When  the  modulating  frequency  is  low  there  are  two  contributions  to  the  linear  electro¬ 
optic  effect.  This  is  a  combination  of  photoelasticity  and  the  Pockels  effect.  The  true  or  primary 
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EO  is  considered  to  be  under  a  constant  stress  condition  (r fjk),  but  because  materials  with  the 
Pockels  effect  are  also  piezoelectric  there  is  a  physical  defonnation  that  must  also  be  considered 
and  can  be  described  in  terms  of  the  elasto-optic  and  piezoelectric  coefficients  which  is  known  as 
the  secondary  EO  contribution.  Adding  these  two  tenns  gives  the  total  electro-optic  effect  in 
tenns  of  constant  stress. 

Reported  linear  refractive  index  coefficients  are  in  either  constant  strain,  (j?Jk)  clamped 
condition  above  acoustic  resonance  or  constant  stress,  (r fJk)  low  frequency  below  acoustic 
resonance  [9].  The  electro-optic  constant  stress  tensor  is  actually  a  combination  of  the  constant 
strain  tensor,  the  elasto-optic  coefficient  (p ijim)  and  piezoelectric  ( djkt )  coefficient. 

Yijk  =  rijk  +  Pijlmdjkl  (3-15) 

These  effects  are  related  to  the  impermeability  tensor  through  the  electric  field  and  there  is  also  a 
small  contribution  from  the  non-linear  electro-optic  effect  (Rfjki)-  Combining  all  these  effects 
we  have: 

AEij  =  rijkEk  T  EijklEkEl  (3-16) 

ABa  =  rijkEk  +  Pijki*ki  +  RijkiEkEi  (3.17) 

The  COMSOL  multiphysics  platform  [10]  is  a  useful  tool  for  modeling  the  effects 
discussed  in  this  dissertation.  There  is  no  native  support  for  electro-optic  and  photoelastic  effects 
in  COMSOL.  However,  these  equations  can  easily  be  programed  into  the  interface  and  modeled 
with  this  software.  These  calculations  typically  start  with  the  undeformed  optical  indicatrix 
which  relates  the  three  principle  Zn  axes  and  the  refractive  indices  nn. 
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(3.18) 


Including  linear  electro-optic  effects  for  the  most  general  (Triclinic)  crystal  symmetry  group 
yield  the  expanded  equation. 
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+  2(r61E1  +  r62E2  +  r63E3)Z1Z2  =  1 


For  electro-optic  materials  of  3m  symmetry. 
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(3.21) 


If  the  electric  field  is  only  applied  along  the  Z  direction  then  the  equation  reduces  to: 
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Because  there  are  no  cross-coupling  tenns  then  the  indicatrix  does  not  rotate  due  to  the  electric 
field  and  because  nx  =  n2  the  index  change  can  be  described  in  two  equations  (where  optical 
propagation  is  along  the  y  axis).  The  ordinary  refractive  index  nG  is  parallel  to  the  Zi  direction 
while  the  extraordinary  refractive  index  ne  is  parallel  to  the  Z3  direction. 


n; 
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These  equations  are  often  further  simplified: 
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The  equation  can  be  simplified  by  considering  Figure  3.2  (a).  The  function  fl[x]=l/(l+x) 
matches  f2[x]=l-x  very  well  for  very  small  values  of  x.  Thus  if  the  product  of  n^r13Ez  is 
sufficiently  small  the  equation  becomes  [11]: 

n'x  ~  n20{l  -  n2r13Ez)  (3'26) 


And  once  again  if  n2r13Ez  is  sufficiently  small  the  tenn  under  the  radial  can  be  simplified  from 
0[x]=VT^x  to  f4[x]=l  -  ix. 
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Figure  3.2  Common  equation  simplifications  for  calculating  the  electric  field 
induced  Birefringence  in  an  electro-optic  modulator 


This  simple  method  can  be  used  for  the  Z  direction  index  and  is  summarized  below  for  both. 
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3.3  Senarmont  Compensator  Technique 

The  Senarmont  compensator  technique  is  a  method  used  for  determining  the  low 
frequency  (Max  500MHz)  electro-optical  coupling  of  a  crystal.  Electro-optic  crystals  usually 
require  high  fields  to  produce  an  induce  Birefringence  large  enough  to  measure,  which  makes 
them  difficult  to  test  at  higher  frequencies  with  standard  signal  generators.  This  technique 
overcomes  this  issue  and  allows  for  accurate  testing.  A  typical  setup  is  shown  in  the  figure  below 
[11]. 
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1  ft 

Applied  electrical  signal 

Figure  3.3  Schematic  of  experiment  based  on  Senarmont  compensator 
technique 

The  setup  always  starts  with  a  coherent/monochromatic  light  source  which  is  most  commonly  a 
laser.  The  intensity  of  the  laser  light  is  controlled  by  an  attenuator  to  bring  the  power  down  to  a 
safer  but  adequate  few  milliwatts  of  power.  The  light  is  incident  on  the  sample  under  test  by  way 
of  crossed  polarizers  on  either  side  known  as  the  polarizer  and  analyzer  -  though  they  are 
identical  parts.  The  polarizer  allows  only  light  with  a  certain  alignment  to  pass  through  it 
typically  at  45°  from  the  vertical,  while  the  analyzer  is  set  to  -45°  from  the  vertical.  Under  this 
condition  all  light  is  blocked  if  the  light  is  propagating  along  the  optical  axis  or  if  it  is  isotropic. 
When  the  sample  is  introduced  between  these  devices  its  natural  birefringence  causes  a  small 
rotation  of  the  polarized  light  which  in  turn  results  in  light  output  at  the  analyzer.  This 
birefringence  can  be  controlled  by  applying  a  modulating  signal  to  the  sample.  So  application  of 
an  electric  field  causes  an  amplitude  change  of  the  light  which  reaches  the  photodetector.  The 
quarter  wave  plate  is  used  to  turn  elliptically  polarized  light  into  linearly  polarized  light,  the 
angle  of  which  is  dependent  on  the  crystal  birefringence;  this  effectively  adds  an  optical  bias  of 
jt/4  to  the  modulator.  The  figure  below  represents  the  maximum  optical  output  possible  in  the 
crystal.  An  ac  signal  applied  to  the  linear  portion  of  the  transmittance  sine  response  can  give  a 


minimally  distorted  optical  output.  This  allows  us  to  use  a  dc  type  test  for  small  signal  ac 
analysis. 


Figure  3.4  Transmittance  curve  (in  blue)  for  applied  DC  field  on  sample  in 

crossed  polarizer  configuration  and  effect  of  quarter  wave  plate  on 
optical  biasing  for  ac  modulation  around  the  50%  transmittance 
point 


The  transmission  curve  in  Figure  3.4  has  been  derived  in  many  sources  [12,  13]  and  is  shown 
here.  Considering  first  the  more  simple  case  of  just  the  polarizer,  sample  and  analyzer  the  output 
intensity  (Equation  3.30)  is  a  function  of  the  field  amplitude  A  and  phase  retardation  8.  This  can 
reduced  to  a  more  simile  form  in  Equation  (3.32). 


I„  = 


^(e-'S-l)|  =  y  Or1' 5  -!)(«'■ 5-l) 


(3.30) 


Io  =  y(2-(ei5  +  e-i5))  =  y[  2-2 


eiS  +  e  15 


(3.31) 


A1 

/0  =  y(  2 -2  Cos  5) 


(3.32) 
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Using  the  relation  from  the  trigometric  half  angle  formulas. 


Sin2 


Cos  B ) 


h 


A 2 

—  2(1  -  ( Cos  5))  =  2d2  Sin2 


(3.33) 

(3.34) 


2  A2  is  proportional  to  the  input  light  intensity.  The  modulating  signal  and  offset  is  given  by 
V  —  VQ  +  VmSin  cot.  When  Vo  is  set  to  y  (Vn  gives  the  peak  transmittance)  this  makes  an  offset 
equivalent  to  what  a  quartz  wave  plate  adds,  resulting  in  a  transmission  that  is  biased  around  the 
more  linear  y  point.  This  is  key  to  the  measurement  as  very  large  modulation  signals  or  bias 

points  may  result  in  harmonics.  The  normalized  transmittance  on  a  scale  of  zero  to  one  is  a  ratio 
of  the  output  intensity  Io  to  the  input  Intensity  f. 

Y  =  Sin2  =  Sin2  =  Sin2  +  VmSin<ot)\  (3.35) 

V0  =  \  (3.36) 

j;  =  Sin2  i^Y  (y  +  VmSin  =  Sin2  Sin  (3.37) 


Equation  (3.37)  can  be  further  simplified  by  applying  the  identities  in  Equations  (3.38)  and 
(3.39). 

Sin2  (u)  =  -(l  —  Cos(2u))  (3.38) 

Cos  (u  +  — )  =  — Sin(u )  (3.39) 
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If  «  1  then  the  equation  can  be  reduced  to 
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Sin  o)t 


(3.41) 


So  a  modulating  voltage  Vm  will  linearly  change  the  intensity  of  the  signal  at  the  photodetector. 
The  typical  setup  for  this  configuration  is  a  bulk  transverse  electro-optic  modulator  as  shown  in 
Figure  3.5.  The  sample  consists  of  an  electrode  on  the  a-b  faces  so  that  an  electric  field  can  be 
applied  in  the  c  direction.  The  laser  light  travels  in  the  b  direction  and  the  Birefringence  caused 
by  the  3  and  1  components  of  the  refractive  index  induce  a  phase  change  in  the  polarization 
components.  Figure  3.5  is  an  explanatory  arrangement,  but  other  orientations  can  also  be  used. 


c 


A 


Applied 

voltage 


Figure  3.5  Orientation  of  applied  field  and  optical  propagation  for  transverse 
electro-optic  modulator 

The  modulated  signal  can  be  expressed  simply  in  terms  of  the  photodiode  signal  or  as  the 
actual  electro-optic  r  coefficient.  There  are  several  r  coefficients  depending  on  the  type  of 
material  and  configuration.  This  dissertation  mostly  studies  Perovskite  materials  which  have 
three  electro-optic  coefficients  r13,r33  and  r51.  The  configuration  described  above  can  measure  a 
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combination  of  the  r13  and  r33  coefficients  directly,  this  is  the  effective  rc  coefficient  in  Equation 


(3.42). 

rc  =  r33  +  (^)  ri3  (3.42) 

The  rc  coefficient  can  be  calculated  from  the  thickness,  t  of  the  sample,  the  length  of  the  sample 
in  the  light  propagation  direction,  Z0,  the  wavelength  of  the  laser,  A,  and  the  extraordinary 
refractive  index  ne.  The  coefficient  also  requires  three  experimental  parameters;  the  amplitude 
Vapp  of  the  modulating  voltage,  the  measured  voltage,  Vout  at  the  photodiode  (after  the 
transimpedance  amplifier),  and  the  voltage  Vpp.  Vpp  is  obtained  by  rotating  the  quarter  wave  plate 
and  monitoring  the  voltage  at  the  photodiode.  As  the  plate  is  rotated  the  voltage  will  vary 
sinusoidally;  it  is  the  peak-to-peak  voltage  from  this  measurement  that  is  used  to  find  Vpp  in 
Equation  (3.43). 


2A  / 1\ 

(  vout  \ 

nj  u0/ 

\fpp  '  ^app) 

(3.43) 


3.4  Samples 

Three  different  samples  were  used  for  experimental  measurements  in  this  section;  PMN- 
30%PT;  quartz;  and  PZN-4.5%PT.  Both  PMTPT  and  PZNPT  are  rhombohedral  at  room 
temperature  but  macroscopically  tetragonal  4mm  after  poling  and  are  ferroelectric.  Quartz  is 
crystal  class  32,  it  is  piezoelectric  but  does  not  have  reversible  polarization  thus  it  is  not 
ferroelectric.  The  dimension  and  configuration  of  the  samples  is  summarized  below. 
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Table  3.2  Sample  compositions,  poling,  configuration,  and  dimensions  used  in  piezoelectric 
experiments 


Sample 

Composition 

Poling 

dimension  [mm] 

X 

y 

z 

1 

PMN-30%PT 

[001] 

1.25 

3.00 

2.25 

2 

Quartz 

None 

5.64 

2.48 

18.49 

3 

PZN-4.5%PT 

TOO  11 

2.2 

2.46 

2.0 

Table  3.3  Sample  configurations  and  Symmetries  used  in  piezoelectric  experiments 


Sample 

Electroded  face 

2 

Area  [mnr] 

Thickness  [mm] 

Symmetry  @  R.T. 

1 

xy 

2.25 

4mm 

2 

yz 

5.64 

32 

3 

xy 

5.412 

2.0 

4mm 

3.5  Admittance  Spectrum  and  Resonant  Modes 

The  admittance  spectrum  of  the  samples  was  examined  with  the  HP  4194A  to  determine 
the  resonant  frequencies.  The  admittance  spectrum  is  used  to  detennine  the  resonant  frequencies 
of  the  sample.  When  the  magnitude  of  the  admittance  reaches  a  maximum  this  is  the  resonant 
frequency  and  when  minimum  this  is  the  anti-resonance.  The  resonant  frequencies  indicate 
where  there  is  large  electrical-mechanical  coupling  in  the  sample  and  where  the  sample  physical 
displacement  and  strain  is  very  high.  These  parameters  are  very  important  to  the  study  of 
resonant  enhancement  of  the  optical  transmission.  The  resonant  frequencies  fr,  of  the  extension 
modes  of  a  rectangular  rod  can  be  calculated  with  the  formula  below  from  IEEE  standards  on 
Piezoelectricity  [14].  Where  /  is  the  sample  length,  v  is  the  acoustic  velocity,  rho  is  the  density 
and  S  is  the  elastic  compliance.  The  parameter  n  is  a  positive  odd  integer  used  to  calculate  the 
hannonics. 


n  n 

fAn)  =  Tiv  ~  2ijpSll 


n  =  1,3, 5, ... 


(3.44) 
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For  the  quartz  sample 


p  =  2651 


kg 

m 3. 


5fi  =  1.277xlO_11[l/Pa] 


n 

1 

2  *  18.49  [mm] 

A 

2651 

kg 

m3 

1.277xlO_11[l/Pa] 

fr(X)  =  146.971  KHz 
/r( 3)  =  440.914  KHz 
/r(5)  =  734.857  KHz 
/r( 7)  =  1.0288  MHz 


(3.45) 

(3.46) 

(3.47) 


(3.48) 


The  PMN-PT  and  PZN-PT  sample’s  resonance  frequencies  are  not  easily  identified  using  the 
above  formula  due  to  their  low  dimensional  aspect,  to  identify  these  frequencies  finite  element 
simulations  are  used  in  the  next  chapter. 

3.6  Electro-optic  Measurements 

The  above  samples  were  measured  using  a  crossed  polarizer  technique  to  determine  the 
electro-optic  effect.  The  lock-in  amplifier  (Stanford  Research  Model  SR830)  was  used  to 
measure  the  signal  in  the  low  frequency  range  from  1  mHz  to  102  KHz  and  the  higher  frequency 
range  from  25  KHz  to  200  MHz  was  tested  with  a  Model  SR844  (Stanford  Research).  The  high 
frequency  limiting  factor  of  this  measurement  is  the  transimpedance  amplifier  used  to  amplify 
the  photodiode  current  output  to  a  voltage  which  has  a  bandwidth  limitation  that  is  dependent  on 
the  amplifier  gain. 
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The  PMN-PT  and  PZN-PT  samples  are  ferroelectric  and  have  a  reversible  polarization 
while  Quartz  is  on  piezoelectric  and  the  polarization  cannot  be  reversed.  Measurements  of  quartz 
were  done  to  detennine  if  electro-optic  transmission  enhancement  exists  in  piezoelectric 
materials.  The  figure  below  shows  the  measurement  configuration  for  these  tests;  additional 
detail  on  the  optical  setup  is  given  in  the  Senannont  compensator  technique  section  of  this 
dissertation. 


Figure  3.6  Block  diagram  of  electro-optic  measurement  setup 

The  electro-optic  coefficients  of  PMN-PT  and  PZN-PT  are  not  well  reported  on  in  the  literature 
and  are  only  really  available  for  a  few  compositions.  Nevertheless,  the  actual  coefficients  are 
likely  very  similar  to  that  of  known  samples  with  similar  concentrations  of  PT.  The  electro-optic 
coefficients  of  PZN-12%PT  are  reported  to  be  rn=7  [pm/V],  r33=134  [pm/V]  and  r5i=462  [pm/V] 
as  reported  in  [15].  The  electro-optic  coefficients  of  PMN-38%PT  are  reported  to  be  rn=25 
[pm/V],  r33=70  [pm/V]  and  r5i=558  [pm/V]  as  reported  in  [16].  The  electro-optic  coefficients  of 
Quartz  (SiCk)  are  reported  to  be  rn=  -0.47  [pm/V],  rzu=0.20  [pm/V]  for  the  wavelength  range  of 
409  to  605  nm  as  reported  in  [9].  The  electro-optic  matrices  for  the  two  crystal  classes  are  shown 
in  Equation  3.49  below.  These  matrices  relate  the  electric  field  Ei,  E2  and  E3  to  the  six 
components  of  the  impermeability  ABn  (where  n  is  an  integer  from  1  to  6). 
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The  sample  was  measured  in  two  different  configurations;  i.e.,  with  the  optical  propagation  along 
the  long  dimension  (18.49  mm)  of  the  sample  and  the  other  along  the  shortest  dimension  (2.48 
mm)  (See  Figure  3.7). 


Figure  3.7  Electro-optic  test  configuration  used  for  quartz  samples  (shaded 
areas  are  electroded) 

The  electro-optic  response  of  the  quartz  sample  (Sample  #2)  is  shown  in  Figure  3.8.  The  rn 
coefficient  was  calculated  using  Equation  3.43  to  be  approximately  1  [pm/V]  which  is  much 
higher  than  the  reported  value  of  0.47  [pm/V]  which  is  actually  the  low  frequency  (constant 
stress)  response.  Both  of  these  configurations  will  test  the  rn  coefficient,  differences  in  the  two 
are  partly  the  result  of  measurement  error  and  partly  from  the  resonant  nature  of  the  material. 
One  configuration  may  show  higher  values  than  the  other  based  on  the  main  direction  of  strain  at 
that  particular  mode.  Also  shown  in  the  plot  is  the  admittance  spectrum  of  the  sample  (blue 


trace).  The  resonant  modes  tend  to  be  very  small  because  of  the  small  piezoelectric  coupling 
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found  in  quartz,  however  the  electro-optic  response  is  very  complex,  clearly  showing  a  complex 
response  due  to  resonance.  Figure  3.9  displays  the  same  results  of  Figure  3.8  but  zoomed  in  on 
the  strongest  resonance.  The  electro-optic  response  correlates  well  with  the  admittance  response 
and  in  this  case  the  peak  electro-optic  enhancement  seems  to  line  up  with  the  mid-way  point 
between  the  resonance  and  anti-resonance.  This  is  the  point  of  maximum  admittance  slope  and  is 
consistent  with  Equation  2.25  as  proposed  in  [1]. 


Quartz  Electro-optic  and  Admittance  tests 
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Figure  3.8  Results  of  quartz  electro-optic  test  showing  measured  electro-optic 
coefficient  and  admittance  response 
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Quartz  Electro-optic  and  Admittance  tests 
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Figure  3.9  Quartz  electro-optic  test  (from  previous  figure)  zoomed  in  on 
largest  resonance 

The  electro-optic  and  admittance  spectrum  of  PMN-30%PT  (Sample  #1)  is  shown  in 
Figure  3.10.  Below  lMz  there  is  a  strong  signal  from  the  admittance  indicating  strong  vibration 
of  the  sample  and  the  EO  response  is  likewise  strong;  above  1MHz  the  sample  becomes  clamped 
but  clearly  the  large  EO  response  continues  as  indicated  by  several  significant  peaks.  This 
demonstrates  that  the  optical  signal  enhancement  may  be  large  even  without  high  strain. 
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PMN-30%PT,  [001]  poled 


1.0m 


500.0m 
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Figure  3.10  Results  of  PMN-30%PT  electro-optic  test  showing  measured 

electro-optic  voltage  as  measured  at  the  photodiode  and  admittance 
response 


3.7  Optical  Transmission 

The  transmission  of  the  electro-optic  signal  is  known  to  increase  near  resonance  when 
using  the  Senarmont  compensator  technique.  This  technique  is  a  measure  of  the  electric  and 
optical  coupling  but  there  has  been  little  work  examining  a  resonating  electro-optic  crystal  with  a 
simple  light  transmission  measurement. 

During  sample  preparation  ferroelectric  samples  become  more  optically  transparent  to 
visible  light  due  to  the  poling  process.  Under  high  electric  field  the  dipoles  in  the  sample  orient 
in  the  direction  of  the  applied  field  and  the  distance  between  domains  becomes  smaller  resulting 
in  increased  optical  transparency.  This  is  particularly  true  in  samples  such  as  PMN-PT  which  is 
not  a  typical  single  crystal  sample.  The  transparency  of  single  crystals  are  governed  primarily  by 
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Admittance  (S) 


the  intrinsic  band  gap  of  the  material,  but  PMN-PT  has  many  domains  which  typically  causes 
low  transparency.  However,  because  of  the  very  small  distance  between  domains  the  overall 
result  is  transparency. 

There  has  been  some  work  in  the  area  of  biased  transmission  measurement  [17]  though 
the  tests  were  not  done  specifically  at  the  piezoelectric  resonant  frequencies.  The  researchers 
modulated  the  samples  with  high  field  and  observed  the  transmission.  The  tests  were  done  from 
DC  to  150  KHz  and  it  was  found  that  this  modulation  increased  transmission  and  that  higher 
modulation  frequencies  resulted  in  higher  transmission.  The  applied  electric  field  needed  a 
threshold  of  about  100  V/mm  to  improve  transmission.  It  is  believed  that  the  response  of  the 
sample  transparency  is  similar  to  what  occurs  during  the  poling  process,  the  difference  being  a 
lower  field  with  modulation  that  allows  the  crystal  domains  to  arrange  for  greater  transparency. 
It  may  be  possible  to  test  the  transmission  in  the  same  way  and  detennine  if  resonance  can 
improve  the  transmission,  perhaps  with  a  lower  operating  field. 

The  increased  transmission  under  electro-optic  testing  using  crossed  polarizers  has  been 
well  demonstrated  in  the  previous  section  however,  there  is  still  a  question  as  to  whether  AC 
biasing  may  have  an  impact  on  the  optical  transmission,  reflection,  or  absorption  of  the  material. 
The  PMN-30%PT  sample  (Sample  #1)  was  tested  from  500  nm  to  900  mn  in  an  optical 
spectrometer  (See  Figure  3.11).  The  sample  was  tested  under  both  normal  conditions  (no  bias) 
and  under  resonance.  The  strongest  resonance  for  this  sample  is  —338.8  kHz,  however  the  results 
show  no  measurable  change  in  the  response  at  this  frequency  or  at  frequencies  close  to  it.  A  few 
conditions  were  retested  to  verify  repeatability.  The  first  non-biased  condition  shows  a  slightly 
larger  response  at  its  peak  but  further  measurements  verified  there  is  no  repeatable  or  measurable 
change  due  to  resonance. 


34 


Optical  transmission  with  resonance  biasing  PMN-30%PT 


Figure  3.11  Results  of  optical  transmission  test  of  PMN-30%PT  under  piezo¬ 
resonant  ac  bias  using  low  fields  (<  50V/cm) 

Although  this  test  did  not  produce  measurable  difference  due  to  resonance  it  may  still  be  possible 
to  observe  increased  transmission  at  resonance.  The  results  of  work  by  the  previous  author  [17] 
as  mentioned  above  showed  that  the  ac  bias  field  must  be  at  least  100  V/mm  (1000  V/cm);  this  is 
of  course  much  higher  than  the  50  V/cm  used  in  Figure  3.11.  Larger  fields  were  not  used  in  this 
test  because  of  the  lack  of  suitable  test  equipment. 


3.8  Summary 

Optical  and  dielectric  measurements  have  been  performed  to  characterize  the  samples. 
The  admittance  measurements  show  the  resonant  and  anti-resonant  frequencies  of  the  samples. 


Strong  fundamental  resonant  response  is  clear  in  the  100  kHz  to  1MHz  region  for  all  the  samples 
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and  hannonics  continue  into  the  low  megahertz  region.  Samples  with  low  aspect  ratio 
dimensions  tend  to  have  resonant  modes  that  are  not  pure;  that  is,  there  is  mixing  between  the 
modes  leading  to  a  very  complex  admittance  spectrum.  Mode  mixing  becomes  more  pronounced 
for  resonant  harmonics  such  that  individual  modes  are  no  longer  distinguishable.  Enhancement 
of  the  electro-optic  effect  was  observed  near  the  resonant  modes  of  most  samples,  especially  near 
modes  where  strain  is  transverse  to  the  direction  of  light  propagation.  Anti-resonant  modes  do 
not  seem  to  contribute  to  enhancement;  this  is  likely  due  to  an  impedance  maximum  at  anti¬ 
resonance  resulting  in  minimum  current  in  the  device.  When  analyzing  the  transmission  versus 
optical  wavelength  spectrum,  sample  resonance  was  not  shown  to  improve  transmission.  These 
experiments  were  performed  with  a  low  electric  field  of  less  than  50  V/cm;  it  may  be  that  a 
higher  field  will  reveal  enhancement  near  resonance. 
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CHAPTER  4:  PIEZORESONANCE  FINITE  ELEMENT  ANALYSIS 


Section  4.2  McIntosh,  R.,  C.  Garcia,  et  al.  (2011).  "Periodically  Poled  Structure  on  Microwave 
Transmissions  Evaluated  by  Scattering  Parameters. "  Integrated  Ferroelectrics 
131(1):  219-229. 

Section  4.3  McIntosh,  R.,  A.  S.  Bhalla,  et  al.  (2012).  Simulation  of  Enhanced  Optical 
Transmission  in  Piezoelctric  Materials.  Advances  and  Applications  in 
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4.1  Finite  Element  Analysis 

Finite  element  analysis  (FEA)  is  a  computational  method  of  solving  real  world  physics 
problems  described  by  partial  differential  equations  (PDE).  Typically  models  based  on  FEA 
consist  of  at  least  one  domain  and  a  PDE  for  that  domain.  There  are  also  boundary  conditions 
defined  on  the  outer  surface  of  the  domain  or  where  one  domain  meets  another.  Once  these  three 
primary  conditions  are  met  the  domain(s)  are  divided  up  into  a  discrete  number  of  pieces  called  a 
mesh.  These  mesh  elements  are  typically  triangular  or  quadrilateral  in  shape  and  are  smaller  near 
areas  of  greater  flux  or  areas  of  the  domain  that  are  smaller.  FEA  solves  the  PDE  on  each 
element  and  then  recombines  them  all  by  way  of  superposition  to  solve  the  problem  on  the  entire 
geometry. 

4.2  Microwave  Transmission  in  Periodic  Ferroelectric  Domain  Structure 

A  combined  experimental  and  numerical  investigation  has  been  conducted  to  study 
piezoelectric  resonance  process  in  relation  to  the  strain  gradient,  the  local  polarization,  and  the 
displacement  current  in  a  piezoelectric  crystal-resonator.  Experimental  discoveries  showed  that 
at  certain  sub-wavelengths  of  the  fundamental  piezoelectric  resonant  frequencies  of  a  given 
resonator,  electro-optic  coefficients  and  also  the  optical  transmission  in  an  otherwise  opaque 
configuration  are  greatly  enhanced.  This  paper  reports  preliminary  study  by  finite  element 
analysis  modeling  on  microwave  transmission  behavior  when  the  piezoelectric  material  is 
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modulated  by  periodic  polarization  distribution  defined  by  piezoelectric  resonance  or  generated 
by  periodic  poling.  Transmission  parameters  are  simulated  in  wide  frequency  range  to 
characterize  wave  propagation  behavior  of  the  material  in  various  local  polarization  conditions 
related  to  piezoelectric  transverse  or  longitudinal  resonance  modes. 

It  was  reported  that  piezoelectric  resonance  can  have  a  significant  effect  on  the  electro¬ 
optic  (EO)  response  of  a  material  [1,  3,  4,  6].  In  the  resonant  condition  an  improved  coupling  is 
observed  (by  way  of  a  dual  frequency  scanning)  at  frequencies  well  above  the  fundamental 
resonant  frequency.  In  a  previous  experiment  an  electro-optic  crystal  (PZN-PT)  was  measured 
using  the  Senarmont  compensator  technique  for  analyzing  the  electro-optic  r  coefficient  at  radio 
frequencies  (<500MHz).  By  applying  a  signal  at  one  of  the  piezoelectric  resonant  frequencies  of 
the  material  while  measuring  the  EO  coefficient  at  various  frequencies  below,  at,  and  above 
resonance,  improved  optical  transmission  was  obtained  near  the  fundamental  and  at  several 
harmonics  above  the  piezoelectric  resonant  frequencies. 

The  precise  reason  for  this  improved  coupling  is  not  completely  known;  however  it  is 
believed  to  be  a  result  of  increased  coupling  of  the  optical  photons  and  acoustic  phonons  brought 
on  by  the  periodic  electric  displacement  associated  with  the  modulation  of  local  polarization.  The 
displacement  is  the  most  pronounced  while  the  sample  is  under  piezoelectric  resonance.  At 
constant  temperature  in  a  piezoelectric  material  the  electric  displacement  (D)  or  polarization  (P) 
is  a  function  of  the  permittivity  (£0£r),  electric  field  (E),  remnant  polarization  (Pro),  and 
additionally  a  strain  (x)  coupled  piezoelectric  component. 

A  sample  at  resonance  due  to  low  frequency  excitation  induces  a  strain  in  the  material 
and  results  in  a  change  in  the  displacement.  When  this  resonance  is  defined  by  one  of  the 
resonant  modes  of  the  sample  leading  to  a  strong  periodic  electric  displacement  in  the  sample  the 
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electro-mechanical  coupling  can  be  very  effective.  Due  to  this  effect  the  transmission  of  an 
electromagnetic  wave  may  be  increased  (or  decreased)  by  the  induced  local  polarization 
modulation  (or  alternatively  by  a  periodically  polarization-inverted  structure).  Under  resonance 
the  sample  deforms  in  response  to  the  applied  signal.  The  polarization  of  the  material  will 
periodically  change  in  magnitude  and  direction  in  a  sinusoidal  manner  as  shown  in  the  top 
portion  of  Figure  4.1.  Similarly  in  a  ferroelectric  material  where  a  high  field  is  used  to  pole  the 
sample  the  polarization  can  be  reversed  in  a  uniform  fashion  with  ideally  positive  and  negative 
polarizations  of  the  same  magnitude. 

In  this  chapter,  we  report  finite  element  COMSOL  multiphysics  modeling  to  simulate  the 
situation  that  a  periodically  poled  ferroelectric  material  inserted  into  a  two  port  transmission  line. 
One  of  the  scattering  parameters  or  the  transmission  parameters  S21,  defined  as  the  forward 
transmission  gain  with  the  output  port  terminated  in  a  matched  load  (expressed  in  decibels 
5'21(r/5)  =  20-log10(|5'21 1),  is  used  to  examine  broadband  solutions  of  the  system  in  order  to 
characterize  the  behavior  of  the  material  with  various  local  polarization  conditions,  especially  the 
conditions  that  are  associated  with  piezoelectric  resonance  modes. 
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Figure  4. 1  Periodically  polarized  antiparallel  structure  within  a  piezoelectric 
sample 

COMSOL,  a  finite  element  analysis  solver  and  simulation  software,  is  used  to  study  the 
effect  of  a  periodically  poled  structure  on  the  scattering  parameters  of  a  rectangular  waveguide. 
In  COMSOL,  when  using  the  RF  module  and  harmonic  propagation  application  mode  the 
polarization  cannot  be  defined  negatively,  thus  a  few  modifications  were  made  to  the  governing 
equations.  This  modeling  method  has  been  carried  out  based  on  charge  balance  assumptions  and 
typically  polarization  values  are  resultant  of  piezo-resonant  activity.  The  displacement  field  can 
be  described  in  terms  of  permittivity  and  electric  field  Z),  =  SyEj.  We  are  concerned  particularly 

with  the  displacement  in  the  transverse  (z  or  3)  direction  Z>3  =  e3  .E ..  The  remnant  displacement 

that  remains  after  a  ferroelectric  sample  is  poled  and  the  field  has  been  removed  provides  an 
offset  to  the  displacement  field.  Including  the  polarization  offset  (P3r0)  we  have  the  formula  as 

D3  —  <s\  jEj  +  P3r0 . 
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The  basic  waveguide  structure  is  shown  in  Figure  4.2.  The  length,  width,  and  height  of 
the  guide  were  assumed  at  2x1x1  cm' .  The  wave  propagates  in  the  x  direction  and  the  transverse 
electric  field  is  along  the  z  direction.  This  structure  uses  2  ports,  the  wave  is  excited  on  the  yz 
face  farthest  to  the  left  (port  1)  and  the  other  port  is  the  yz  face  farthest  to  the  right  (port  2),  the 
other  4  external  surfaces  are  configured  (with  or  without  metallization)  accordingly  for  each 
simulation. 


Figure  4.2  Simulation  setup  for  sample  inserted  in  microwave  waveguide 

showing  electromagnetic  wave  excited  on  port  1  and  travels  to  port 
2  along  the  x-axis 


Figure  4.3  illustrates  five  polarization  conditions  considered  when  the  direction  of  the 
polarization  is  parallel  to  the  z-axis.  Each  polarization-uniform  area  is  referred  as  an  element. 
The  waveguide  dimensions  are  the  same  regardless  of  the  number  of  elements.  The  two  shades 
distinguish  the  polarization  in  the  positive  or  negative  z-directions;  all  conditions  studied  started 
with  a  polarization  in  +z-direction.  The  piezoelectric  material  is  assumed  to  have  zero  field 
impedance  Z0  and  a  zero  field  permittivity  £r(0). 


41 


1  element 


▲ 


’nL- 


2  elements  3  elements 


4  elements  5  elements 


Figure  4.3  Antiparallel  step  polarization  configuration  for  1-5  elements  in 
microwave  FEA  model 

The  transmission  parameter  S21  in  dB  is  simulated  for  transmitted  signal  in  a  frequency  range  of 
0.1  to  40  GHz  in  steps  of  0.25GHz.  The  effects  of  the  number  of  polarized  elements,  magnitude 
of  the  polarization  and  length  of  the  waveguide  on  the  transmitted  signal  are  studied. 

4.2.1  Resonator  Surfaces  Free  From  Metallization 

All  surfaces  of  the  y-  and  z-  direction  of  the  material,  refer  to  Figure  4.2,  are  free  from  any 
metallization  in  this  simulation.  The  four  surfaces  are  scattering  boundaries  while  the  other  two 
faces  (//x-axis)  are  wave  ports.  This  configuration  is  equivalent  to  that  in  which  the  metallization 
is  removed  after  the  sample  is  (periodically)  poled  along  the  z-direction.  Antiparallel  step 
polarization  configuration  with  one  to  five  elements  are  considered,  as  depict  in  Figure  4.3. 

In  Figure  4.4  we  can  see  the  S21  response  in  a  four  element  configuration  over  a  broad  frequency 
range  from  1  to  40  GHz  as  a  function  of  polarization  strength.  The  waveguide  cutoff  frequency 
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is  at  about  7.5  GHz  for  all  conditions,  attributed  to  the  dimension  and  the  impedance  of  the 
sample.  A  drop  in  transmission  is  observed  at  -32  GHz  for  this  sample  of  four  polarization- 
inverted  elements.  Other  than  that,  an  increased  polarization  amplitude  causes  an  amplification 
of  the  transmission;  however  the  overall  shape  of  the  response  as  function  of  frequencies  remains 
the  same  for  all  polarization  levels  above  0.1  pC/cm  . 

Polarization  perpendicular  to  the  microwave  propagation 


- P=  le-3  (C/cmA2) 

- P=  0.1e-3  (C/cmA2) 

- P=  0.1e-6  (C/cmA2) 

- P=  0.1e-9  (C/cmA2) 

- P=  0.1e-12  (C/cmA2) 

- P=0 


Figure  4.4  Transmission  parameter  as  function  of  frequency  in  4  periodically 
polarization-inverted  elements 

Figure  4.5  shows  the  polarization  at  0.1  pC/cm  up  to  an  artificially  high  value  1  C/cnT 
and  how  the  response  changes  by  varying  the  number  of  inverted  polarization  elements.  When 
the  polarization  level  is  0.1  pC/cm  even  and  odd  number  of  elements  show  two  groups  of 
transmission  behavior  at  below  the  7.5  GHz  cutoff  frequency.  When  the  polarization  level  is  0.1 
nC/crrT  different  number  of  polarization-inverted  elements  has  a  significant  change  in  response. 
The  shape  of  the  S21  for  the  polarization  level  from  0.1  pC/cm  to  1  C/cm"  are  similar  indicating 
the  mechanisms  are  stable  at  the  polarization  level.  Using  5  elements  one  obtains  a  large 
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transmitted  wave  response  from  18  to  25GHz  -  indicating  that  transmission  amplification  over 
certain  bands  can  be  achieved  by  selecting  the  proper  number  of  elements. 


Polarization  perpendicular  to  the  microwave  propagation  direction 


lC/cm2 


— 1  element 
— 2  elements 
— 3  elements 
— 4  elements 
— 5  elements 


O.lxlO'6  C/cm2 


O.lxlO"9  C/cm2 
O.lxlO12  C/cm2 


Frequency  (GHz) 


Figure  4.5  Transmission  parameter  as  function  of  frequency  and  the  number 
of  periodically  polarization-inverted  elements 

Figure  4.6  is  a  visual  representation  of  the  transverse  component  of  the  electric  field  on 
each  boundary  at  20GHz  when  the  polarization  is  0.1  nC/cirT  (refer  to  Figure  4.3).  In  a  dielectric 
waveguide  the  number  of  transverse  electric  modes  along  the  z  direction  will  increase  with 
frequency  but  is  in  general  dependent  on  the  dimension,  permittivity,  and  boundary  conditions. 
However  in  the  figure  the  only  difference  is  the  polarization  with  all  other  conditions  being 
equal.  The  number  of  modes  supported  in  the  propagation  direction  tends  to  be  the  same  as  the 
number  of  polarization-inverted  elements.  The  unpolarized  waveguide  has  approximately  four  or 
five  modes  at  20  GHz  so  that  modes  close  to  this  are  more  readily  supported,  that  is  the  case  for 
three,  four,  and  five  elements. 
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Figure  4.6  A  surface  plot  of  the  transversal  electric  field  strength  at  20GHz 

(P  =0.1  nC/cnr)  that  shows  how  the  polarization-inverted  sections 
tend  to  support  transverse  electric  modes  corresponding  to  the 
number  of  the  elements 

4.2.2  Comparison  of  gradient  and  step  polarization,  surfaces  free  from  metallization 

Most  of  our  simulations  used  a  periodically  inverted  polarization  which  simply  changes 
direction  but  with  the  same  magnitude  throughout  (like  a  square  wave);  however  in  the  case  of 
piezoresonance  the  periodic  polarization  change  would  be  more  sinusoidal  in  nature.  For  this 
reason  we  compare  a  periodical  step  polarization  with  a  sinusoidal  gradient  polarization,  both  of 
which  share  the  same  maximum  amplitude  (yet  the  step  polarization  has  an  RMS  polarization 
value  0.1  nC/cm2  and  the  sinusoidal  is  0.0707  nC/cm2).  From  Figure  4.8  it  is  clear  that  below 
about  25GHz  gradient  polarization  distribution  is  more  effective  in  transmission  amplification 
and  the  opposite  is  true  above  25GHz;  particularly  at  about  37  GHz  gradient  polarization  can  be 
seen  with  improved  transmission 
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Figure  4.7  Gradient  and  step  polarization  in  a  5  element  cavity,  both  with 

2 

magnitude  0.1  nC/cm 


Step  P=  le-3  (C/cmA2) 
Gradient  P=  le-3  (C/cmA2) 
Step  P=  0.1e-3  (C/cmA2) 
Gradient  P=  0.1e-3  (C/cmA2) 

- Step  P=  0.1e-6  (C/cmA2) 

Gradient  P=  0.1e-6  (C/cmA2) 

- Step  P=  0.1e-9  (C/cmA2) 

Gradient  P=  0.1e-9  (C/cmA2) 

- Step  P=  0.1e-12  (C/cmA2) 

Gradient  P=  0.1e-12  (C/cmA2) 
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Figure  4.8  Comparison  of  sinusoidal  and  step  polarization  distribution  in  a 
polarization-inverted  5 -element  waveguide 
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4.2.3  Effect  of  the  sequence  of  the  polarization,  surfaces  free  from  metallization 

Although  unexpected,  by  changing  the  order  of  the  +/-  z-direction  polarization  the 
response  does  change  slightly  within  the  resolution  of  the  simulation.  The  configuration  is 
defined  in  Figure  4.9  using  four  elements.  The  sequence  of  the  polarization  makes  almost  no 
difference  for  large  polarization  (greater  than  or  equal  to  0.1  pC/cm  )  at  high  frequencies; 
however  at  lower  polarization  the  order  can  affect  the  low  frequency  response,  as  illustrated  in 
Figure  4.10. 
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Figure  4.9  FEA  model  setup  showing  periodic  poling  configuration.  “P” 

corresponds  to  polarization  in  the  positive  z  direction  and  “N”  in 
the  negative  z  direction 
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Figure  4.10  Transmission  parameter  as  a  function  of  frequency  for  4 

polarization-inverted  elements  while  the  order  of  antiparallel 
polarization  can  cause  a  significant  shift  in  the  response 

4.2.4  Dielectric  Waveguide  with  Top  and  Bottom  Metallization 

Piezoelectric  resonators  are  typically  electroded  to  allow  the  application  of  electrical 
signals.  Simulation  result  with  top  and  bottom  boundaries  of  the  waveguide  as  perfect  electrical 
conductors  (PECs)  and  the  other  two  sides  as  scattering  boundaries  is  shown  in  Figure  4.11. 
When  the  polarization  is  <0.1  pC/cm  the  most  significant  changes  in  response  due  to 
periodically  polarization-inverted  elements  is  seen  at  below  the  cutoff  frequency.  When  the 
polarization  is  0.1  nC/crrT  or  higher  the  response  is  strongly  dependent  on  the  number  of 
elements.  The  overall  response  of  this  condition  is  similar  to  the  simulation  results  of 
metallization-free  conditions,  which  indicates  that  the  metallization  does  not  screen  the 
polarization  effect  under  the  given  conditions.  These  results  are  consistent  with  results  reported 
earlier  [1,  3,  4,  6]  where  samples  electroded  in  a  similar  manner  displayed  wave  enhancement 
due  to  an  induced  periodic  polarization. 
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Figure  4.11  Transmission  parameter  as  a  function  of  frequency  for  a  number 
of  periodically  polarization-inverted  elements  while  electrodes 
and  polarization  are  perpendicular  to  the  microwave  propagation 
direction 

4.2.5  Polarization  in  line  with  wave  propagation,  surfaces  free  from  metallization 

Considered  here  is  the  effect  of  the  polarization  along  the  axis  of  wave  propagation 
(Figure  4.12).  When  polarization  is  small  it  has  little  to  no  effect  on  the  response  -  thus  it  is  not 
shown  here.  As  seen  in  Figure  4.13,  little  effect  is  seen  when  polarization  is  <0.1  nC/cm  while  it 
causes  significant  changes  on  low  frequency  response.  Polarization  of  0.1  pC/cm  and  higher 
increases  transmission  significantly  over  a  broad  range  (with  exception  at  frequency  near  ~16 
GHz);  however  S21  is  in  general  much  lower  compared  to  the  situation  where  polarization  is 
perpendicular  to  the  wave  propagation  direction. 

Figure  4.14  shows  surface  plots  of  the  electric  field  strength  at  20  GHz  (Pr  =0.1  pC/cnr, 
Pr//x-axis)  for  the  number  of  polarization-inverted  elements  from  0  to  5,  while  polarization  is 
parallel  to  the  microwave  propagation  direction.  It  appears  that  co-linear  polarization  both  with 
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and  against  the  microwave  propagation  has  the  effect  of  confining  the  guided  microwave  thus 
increases  the  transmission. 


Figure  4.12  FEA  microwave  model  setup  of  antiparallel  polarization  in  the 

direction  of  propagation  (positive  x-axis)  and  against  propagation 
(negative  x-axis) 


Figure  4.13  Transmission  parameter  as  a  function  of  frequency  for  a  number 
of  periodically  polarization-inverted  elements  while  polarization 
is  parallel  to  the  microwave  propagation  direction 
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Figure  4.14  Surface  plots  of  the  electric  field  strength  at  20  GHz  (Pr  =0. 1 
uC/cm2,  Pr//x-axis)  for  the  number  of  polarization-inverted 
elements  from  0  to  5,  while  polarization  is  parallel  to  the 
microwave  propagation  direction. 


4.2.6  Polarization  parallel  (with)  and  antiparallel  (against)  propagation,  surfaces  free  from 
metallization 

Using  one  element  and  once  again  the  effect  of  polarization  with  (parallel)  and  against 
(antiparallel)  propagation  is  considered.  The  direction  of  polarization  does  affect  the 
transmission  parameter  however  the  effect  is  appreciable  in  certain  range  of  polarization 
magnitude  (as  seen  in  Figure  4.15,  1  nC/cm  to  1  pC/cm  ).  Respective  plots  of  parallel  and 
antiparallel  polarizations  >1  uC'/cm  are  on  top  of  each  other,  indicating  that  above  a  certain 
threshold  the  polarization  direction  does  not  impede  or  enhance  the  transmission  of  the 
microwave. 


51 


parallel  P=le-5  (C/cmA2) 
anti  para  I  lei  P=  le-5  (C/cmA2) 
parallel  P=le-6  (C/cmA2) 
anti  para  I  lei  P=  le-6  (C/cmA2) 
parallel  P=le-7  (C/cmA2) 
anti  para  I  lei  P=  le-7  (C/cmA2) 
parallel  P=le-8  (C/cmA2) 
antiparallelP=  le-8  (C/cmA2) 
parallel  P=le-9  (C/cmA2) 
anti  para  I  lei  P=  le-9  (C/cmA2) 
- P=0 


Figure  4. 15  Transmission  parameter  as  a  function  of  frequency  and  magnitude 
of  polarization  for  a  single  element  with  polarization  parallel  or 
antiparallel  to  the  microwave  propagation  direction 

4.2.7  Summary 

Finite  element  analysis  method  is  used  to  simulate  microwave  propagation  behavior  in 
terms  of  transmission  parameter  S21,  in  a  ferroelectric  waveguide  with  periodic  polarizations. 
The  periodic  polarization-inverted  sections  (elements)  are  generated  either  dynamically  at  given 
piezoelectric  resonance  frequencies  or  formed  by  patterned  poling.  Analytical  solution  for  such 
nonlinear  polarization-dependent  system  is  not  available  and  numerical  method  permits  detailed 
study.  The  simulation  is  conducted  for  1-40  GHz  and  for  periodic  polarization  in  transverse  or 
longitudinal  of  the  wave  propagation  directions.  Model  ferroelectric  waveguides  containing  up 
to  five  polarization-inverted  sections  were  simulated.  It  is  verified  that  periodically  poled 
structure  in  ferroelectric  waveguide  can  contribute  to  a  significantly  enhanced  transmission  of 
microwave.  At  a  given  polarization  magnitude,  transmission  amplification  over  certain  bands  can 
be  achieved  by  selecting  the  proper  number  of  elements.  Metallization  in  transverse  direction 
was  found  insignificant  in  screening  the  polarization  effect;  the  simulated  transmission 
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enhancement  is  much  stronger  for  polarization  in  transverse,  than  in  longitudinal,  directions  of 
the  propagation  wave. 

It  is  noted  that  the  simulation  results  are  consistent  with  previous  experimental  results 
that  apparent  electro-optic  coefficients  or  optical  transmission  are  enhanced  at  frequencies 
related  to  piezoelectric  resonance  of  certain  modes.  Further  study  will  be  useful  to  extend  the 
simulation  to  optical  wavelength  and  to  include  nonlinearity  of  dielectric  properties  in  the 
simulation  model. 

4.3  Frequency  Domain  Piezoelectric  Model 

4.3.1  The  Sample 

The  material  used  was  (l-x)Pb(Zn1/3Nb2/3)03-(x)PbTi03  where  x  =  0.045  (PZN-4.5%PT)  which 
is  a  single  crystal  poled  along  <00 1>.  At  room  temperature  this  composition  is  rhombohedral, 
however  after  poling  it  is  known  to  be  pseudo-tetragonal  (4mm  symmetry)  [18].  The  sample  is 
electrode  on  its  (001)  faces  (Figure  4.16).  The  electrode  area  =  3.21  mm  x  2.83  mm=9.0843  xlO- 
6  m2,  thickness=6.3x!0-3  m.  The  incident  optical  signal  propagates  along  the  y  axis. 
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Applied 

voltage 


Figure  4.16  PMN-30%PT  sample,  poled  in  [001]  direction.  Electroded  on  top 
and  bottom  {001 }  surfaces  of  the  x  by  y  dimensions 

4.3.2  The  Piezoelectric  Model 

The  piezoelectric  model  is  used  to  identify  the  primary  resonant  modes  of  the  sample  and 
eventually  when  applied  to  the  optical  model  to  calculate  the  refractive  index  changes.  Equation 
(1)  is  the  partial  differential  equation  used  to  describe  the  relationship  between  the  density  (p), 
angular  frequency  (co),  mechanical  displacement  (u)  and  stress  (X)  in  the  piezoelectric  domain 
(for  a  time  harmonic  solver).  Additionally  in  the  domain,  the  divergence  of  the  electric 
displacement  (D)  field  is  equal  to  the  charge  density  ( pv ).  The  electric  field  is  equal  to  the 


gradient  of  the  negated  voltage  and  the  strain  (X)  is  related  to  the  mechanical  displacement. 


u-V  -X  =  Fve1<p 

(4.1) 

VD  =  pv 

(4.2) 

tq 

ii 

i 

< 

(4.3) 
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(4.4) 


[(v«)r 


+  Vu 


This  simulation  uses  the  converse  piezoelectric  effect,  meaning  an  applied  electric  field  causes  a 
change  in  the  strain.  In  this  case  in  equations  (4.5)  and  (4.6)  the  Stress-charge  fonn  is  used. 
Respectively  the  terms  X,  CE,  AC,  e,  s0,  £rS,  Dr  are  the  stress,  elasticity  matrix,  strain,  coupling 
matrix,  vacuum  permittivity,  relative  permittivity,  and  remanent  displacement  field.  The  ACj  tenn 
is  the  initial  stress  in  the  material  and  is  assume  to  be  zero. 


Stress-Charge  form: 


X  -  cE(x-xi)-erE  +  X 

(4.5) 

D  =  e(x- x,. )  -  e0  srSE  +  D 

(4.6) 

The  pertinent  material  parameters  for  a  piezoelectric  material  with  4mm  symmetry  are  shown  in 
Equation  (4.7).  The  symmetry  is  applied  to  the  elasticity  matrix  (cy),  the  piezoelectric  coupling 
matrix  (ey)  and  the  pennittivity  matrix  (sy). 


Table  4.1  Shows  the  material  parameters  used  in  the  piezoelectric  simulation.  Literature  values 
are  from  Zhang  el  al.  [19] 
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Table  4.1  Electrical  and  mechanical  material  parameters  of  the  samples 


Elasticity  (@  constant  E-fded) 

[N/m2]* 

*  *  2 
Coupling  matrix  [C/m"]* 

Permittivity 
(@  constant 
temperature)* 

Density  [kg/m3]* 

Cn=11.7xl010 

e15=13.6 

eU  =  3307 

8038.4 

c12=10.3xl0lu 

e3i=-2.4 

=  1242 

Ci3=10.1xl0lu 

e33=27.1 

c33=10.8xl0lu 

e"  =  1 

c44=7.1x101u 

c66=6.6x101u 

*  Reported  values  for  the  Elasticisy,  coupling,  real  permititiy,  and  density  are  from  Zhang  et  al. 

The  imagonary  permititiy  is  a  measured  values  for  this  sample. 

A  0.5  volt  AC  source  is  used  to  excite  the  sample  and  the  bottom  domain  is  grounded. 
The  remaining  four  boundaries  are  all  under  zero  charge  condition.  The  crystal  is  completely 
unconstrained  on  all  boundaries,  edges  and  points.  In  order  to  calculate  the  response  of  the 
electrical  conductivity  the  real  (Ir)  and  imaginary  (I,)  parts  of  the  current  were  integrated  over  the 
surface  of  the  top  boundary  (the  boundary  of  the  electrode).  Using  these  values  and  combing 
with  the  applied  AC  voltage  (Vapp=0.5)  gives  the  conductance  (G)  in  equation  (4.8)  and 
susceptance  in  (B)  equation  (9).  Then  the  Magnitude  of  the  admittance  can  easily  be  calculated 
in  equation  (10)  and  the  phase  angle  in  equation  (4. 1 1) 


G  =  A/C, 

(4.8) 

(4.9) 

|  Y  |=  <Jg2+B2 

(4.10) 

Phase=  Tan~\B / G) 

(4.11) 
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The  results  of  the  admittance  calculation  are  shown  in  Figure  4.17,  resonance  frequencies  are 
indicated  by  peaks  in  the  |Y|  and  minimums  indicate  anti  resonate  frequencies.  Additionally  the 
sample  was  also  tested  experimentally  and  a  comparison  is  show  in  Figure  4.18.  The  response 
does  not  line  up  well  however  the  amplitudes  are  on  a  comparable  scale  and  tests/simulations  on 
similar  compositions  of  higher  dimensional  aspect  ratio  have  shown  extremely  close  agreement. 


PMN-30%PT  Admittance  Simulation 
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.  Phase 


Frequency  (Hz) 
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Figure  4.17  Simulation  of  Admittance  magnitude  |Y|  in  Siemens  and  the  phase 
angle  in  degrees  over  frequency  range  encompassing  the 
piezoelectric  vibrational  frequencies  below  1  MHz 
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Figure  4. 18  Comparison  of  the  simulation  and  experimental  results  of  the 
admittance  frequency  spectrum  of  PMN-30%PT 

4.3.3  Frequency  Domain  Optical  Model 

The  optical  model  consists  of  two  parts  one  being  a  three  dimensional  model  identical  to 
the  piezoelectric  model  described  above  and  the  other  being  a  two  dimensional  model  used  to 
compute  the  transmission  parameters  of  the  optical  transmission  through  the  center  of  the 
piezoelectric  crystal  (See  Figure  4.19).  Linear  extrusion  coupling  variables  are  used  to  map  the 
refractive  index  changes  from  the  3D  to  the  2D  model 

This  model  uses  a  slightly  different  sample,  with  electrode  area  of  2.26  mm  x  2.52 
mm=5.6952  xlO-6  m2,  and  thickness=2. 05x10-3  m.  This  sample  is  PZN-4.5%PT  and  possesses 
the  same  pseudo  tetragonal  4mm  symmetry  as  the  previous  sample  but  with  higher  opacity 
lending  itself  well  to  optical  testing. 
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Includes  2  models 


Model  1: 3D  piezoelectric 
Model  2: 2D  optic 


Linear  extrusion 
coupling 


This  2D  plane  in  the  3D  model  is  the  exact  same 
dimension,  position  and  orientation  as  in  the  2D  model 


Difficult  to  model  optical  wave  in  3D  because  optical 
path  length  is  mm  compared  to  wavelength  in  nm 


Linear  extrusion  coupling  variable  used  to  apply  the  simulated 
refractive  index  from  the  3D  model  to  the  2D  model 


Figure  4.19  Illustration  of  the  finite  element  optical  model  that  includes  two 
sub  models  to  handle  the  coupling  between  piezoelectric  and 
optical  interaction 

The  Optical  model  is  rectangular  with  the  longer  dimensions  being  the  direction  of  wave 
propagation  (y-axis).  The  side  boundaries  are  perfect  magnetic  conductors  allowing  plain  wave 
propagation.  The  end  boundaries  are  ports,  one  port  is  the  output  and  the  other  is  used  to  launch  a 
5mW  wave  polarized  in  the  x  and  z  directions  (transverse  to  propagation).  The  Propagation 
constant  is  P=(2*7t)/(3xl08[m/s]/2.857xl014  [FIz]),  where  2.857xl014  Hz  is  the  frequency  of  a 
1050  nm  wave  in  free  space.  The  partial  differential  equation  applied  to  the  domain  is: 

c  .  \ 


V  x  ju~x  (V  x  f )  -  k{ 


a>£, 


E  =  0 


(4.12) 


o  7 


Where  [ir  ,  £r  ,£0  ,  k0  ,  a,  a>,  E  are  respectively  the  vacuum  permeability,  relative 
permittivity,  vacuum  permittivity,  free  space  wave  number,  conductivity,  angular  frequency  and 
the  electric  field.  The  refractive  index  changes  due  to  the  applied  electric  field  and  the  shape 


changes  of  the  piezoelectric  sample  are  modeled  by  way  of  the  electro-optic  and  elasto-optic 
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effects  respectively.  The  electro-optic  effect  is  the  change  of  the  refractive  index  due  to  an 
applied  electric  field  (E)  and  at  constant  strain  (the  clamped  condition)  is  represented  by  rx  the 
elasto-optic  effect  is  the  change  in  refractive  index  due  to  a  strain  (x)  and  is  represented  by  p. 
Combining  these  effects  and  defining  in  tensor  form  the  impermeability  tensor  By  is: 

=  rijkEk  +  pmxkl  (4.13) 

A By  can  be  defined  in  terms  of  the  refractive  index  with  zero  applied  field  and  zero  stress  ( nl} ) 
and  the  new  refractive  index  «.. 


A Bj  =  A| 


(l) 

M 

2 

(V 

UJ 

UJ 

l  n«, 

(4.14) 


By  combining  equations  (4.13)  and  (4.14),  converting  to  matrix  form  and  applying  4mm 


symmetry  the  matrix  for  the  six  independent  refractive  indices  is: 
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(4.15) 


Simplifying  gives  the  following  equation  which  is  applied  to  the  COMSOL  model  to  fully 
describe  the  coupling  of  the  refractive  index  values  between  the  piezoelectric  to  the  optical 
models. 
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—  +  r]3E3  +  Pi  ,x, ,  +  pl2x22  +  pux33 

nu 

4"  +  rB^3  +  P\2XU  +  PnX22  +  Pi  3*33 
n22 

2  *33^3  7^3 l*fi  1  Pi lX22  PiiXii  ^  ^ 

nii 

n2i  ~  V ylr5l^2  P \\X 22 

n2l  =l/ylr5lEl+P44X2l 

nU  =1/V/,6  6^12 

The  elasto-optic,  electro-optic  and  refractive  indices  for  the  material  were  taken  from  the 
literature  and  are  listed  in  Table  4.2.  The  refractive  indices  and  elctro  optic  coefficients  are 
actually  for  a  similar  composition  (PZN-12%PT).  The  elasto  optic  coefficients  are  for  PZN-PT 
and  where  actually  converted  from  the  piezo  optic  coefficients  in  the  literature  with  the  equation 

p  —  K  c 

r  mn  mp  pn  • 


Table  4.2  Optical  parameters  for  PZN-PT 


Refractive  indices’ 

Elasto  optic  coefficients# 

ne 

2.57 

Pn 

0.0888 

no 

2.46 

Pl2 

0.0816 

P 1 3 

0.0808 

Electro-optic  coefficients’ 

P31 

-1.032 

r[3 

7  pm/V 

P33 

1.911 

r33 

134  pm/V 

P44 

0 

r5i 

462  pm/V 

P66 

0 

*  from  Lu  et  al.[16] 

#  from  Lu  et  al.  [15] 


The  results  of  the  optical  field  simulation  can  be  seen  in  Figure  4.20  below.  The  plane  wave 
propagates  along  the  entire  length  of  the  sample  domain.  The  launched  wave  is  1050  mn  in  free 
space  but  for  a  refractive  index  of  about  2.5  the  wave  is  420  mn  in  the  material  -  this  confirms 
that  the  meshing  on  the  domain  is  dense  enough  to  represent  the  wave. 
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Figure  4.20  Simulation  results  of  the  optical  plane  wave  propagating  in  the  y 
direction  and  a  plot  of  the  electric  field  (Ez)  in  the  transverse  z 
direction 


In  Figure  4.21  the  transmission  of  the  optical  signal  using  the  transmission  parameter  S21  is 
shown  with  the  admittance  simulation.  The  figure  clearly  shows  an  increase  the  transmission 
near  the  resonant  frequencies,  typically  the  transmission  follows  the  slope  of  the  admittance 
magnitude  and  is  quite  consistent  with  typical  experimental  results  for  a  sample  near  the  resonant 
frequencies!!]. 
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Figure  4.21  Simulation  results  of  the  frequency  response  of  the  optical 

response  over  a  range  of  low  frequency  resonant  modes  showing 
strong  correlation  with  the  slope  of  the  electrical  admittance 

At  resonance  the  current  field  inside  the  material  becomes  periodic  and  leads  to  an  induced 

polarization  change  inside  the  sample  at  integer  multiples  of  the  fundamental  resonate  modes.  As 

proposed  by  Johnson  et  al.  [4]  it  is  this  periodically  changing  field  that  has  a  large  influence  on 

the  transmission  and  our  model  of  the  combined  influence  of  the  electro-optic  and  elasto-optic 

shows  these  effects.  As  illustrated  in  Figure  4.22  the  current  under  certain  resonant  modes  (in 

this  case  at  215  kHz)  promotes  a  refractive  index  gradient  that  is  sufficient  to  cause  a  large 

increase  in  transmission  magnitude.  Both  acoustic  and  optical  phonons  can  exist  in  the  crystal. 

Polarization  is  present  in  the  transverse  optical  phonons  (TO)  which  couples  with  the  displacive 
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current.  When  at  resonance  the  TO  inodes  become  periodic  and  result  in  a  complete  circuit  of 
current  flow  in  the  local  field. 
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Figure  4.22  Current  density  distribution  in  the  sample;  off  piezoelectric 

resonance  at  213  KHz  (Left)  and  near  piezoelectric  resonance  at 
215  KHz  (Right). 


At  certain  resonance  modes  the  current  becomes  periodic  and  leads  to  an  induced  changing 
charge  distribution.  At  resonance  more  of  the  current  is  co-propagating  with  the  optical  wave  and 
enhancement  is  observed.  The  y-axis  is  the  direction  of  propagation  and  z-axis  the  wave 
polarization  direction 

It  has  been  shown  that  under  certain  conditions  of  piezoelectric  resonance  the  optical 
response  of  the  material  can  be  tuned.  Experimental  results  on  electro  optic  crystals  show  high 
correlation  between  crystal  vibrational  modes  and  electro  optical  transmission  signal.  A 
COMSOL  model  was  developed  to  describe  the  low  frequency  crystal  admittance  response,  the 
combined  effects  of  electro-optic  and  elasto-optics  and  transmission  of  an  optical  wavelength 
plane  wave.  The  COMSOL  model  is  in  good  agreement  with  experimental  results  for  both 


vibration  modes  and  optical  response.  Previously  proposed  models  used  to  explain  the  influence 
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of  bound  charge  distribution  on  wave  transmission  have  been  given  further  credence  by 
multiphysics  modeling  of  the  electric  displacement  current  influence  on  the  propagating  optical 
wave. 

At  resonance  by  definition  the  impedance  of  the  sample  is  at  a  minimum,  meaning  with  when 
using  a  voltage  source  there  is  a  large  influx  of  current  at  resonance  and  conversely  a  large 
impedance  at  anti-resonance  gives  a  very  low  current. 

4.4  Frequency  Domain  Piezoelectric  Model  -  High  Aspect  Ratio 

Section  4.3  modeled  a  real  piezoelectric  sample  that  had  a  geometry  similar  to  a  cube. 
Due  to  this  low  aspect  ratio  between  the  dimensions  individual  vibrating  modes  are  not  pure.  To 
observe  these  modes  a  high  aspect  ratio  bar  of  PMN-30%PT  was  modeled.  The  sample  was 
configured  with  its  x  and  y  dimension  electroded  giving  it  an  area  of  20  x  2  mm  and  a  thickness 
of  0.5mm.  This  is  a  typical  configuration  for  a  k3i  bar.  The  sample  parameters  are  listed  in  the 
table  below. 


Table  4.3  Sample  electrical  and  mechanical  material  parameters  from  [20] 


Compliance  (@  constant  E-filed) 

10'12  [m2/N] 

Coupling  matrix 

10'12  [CN]* 

Permittivity 
(@  constant 
temperature) 

Density  [kg/m3] 

sn=  52 

O- 

1l 

SO 

o 

e\  4  =  3600 

8040 

s12=  -18.9 

d3j=  -921 

£33  =  7800 

s13=  -31.1 

d33=  1918 

s33=  67.7 

s44=  14.0 

s66=  15.2 

The  approximate  resonate  frequency  and  hannonics  of  the  bar  can  be  calculated  fonn 
equation  (3.1),  they  are  summarized  in  Table  4.4.  Figure  4.23  shows  the  response  of  the 
admittance  magnitude  and  phase  angle  frequency  spectrum.  The  figure  shows  the  fundamental 
resonance  at  38.750  KHz  and  harmonics  at  115.250,  189.750  KHz.  The  harmonics  are  easily 
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identified  because  they  are  successive  odd  multiples  of  the  fundamental  i.e.  38.75  KHz  *  3  = 
116.25  KHz,  38.75  KHz  *  5  =  193.75  KHz,  38.75  KHz  *  7  =  271.25  KHz.  The  slight  deviation 
from  this  trend  is  due  to  a  mixing  of  resonance  mode  vibration  with  other  vibration  directions. 
As  compared  in  Table  4.4  the  simulation  and  calculated  values  agree  well.  There  are  also  anti¬ 
resonance  frequencies  associated  with  each  resonance  that  are  a  slightly  higher  in  frequency.  By 
using  the  fundamental  resonant  frequency  as  fm=38.75KHz  and  anti-resonance  at  fn=43.0  KHz 
the  coupling  coefficient  k3i  can  be  calculated  from  Equation  (3.8)  Using  these  values  k3i  =0.482 
and  calculating  directly  from  the  material  parameters  gives  k3i  =  0.486  which  matches  very  well. 


PMN-30%PT  piezoelectric  resonant  fundamental  and  harmonics 


Figure  4.23  Admittance  frequency  spectrum  of  typical  piezoelectric  resonator 
showing  resonance  and  anti-resonance  frequencies  of  fundamental 
and  harmonics 
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Table  4.4  Comparison  of  the  piezoelectric  resonate  frequencies:  calculated  using  a  formula  and 
derived  by  finite  element  analysis 


n 

Formula  (KHz) 

COMSOL  FEA  (KHz) 

1 

38.664 

38.750 

3 

115.993 

115.250 

5 

193.322 

189.750 

freq(76)=38750  Surface:  Total  displacement  (m) 

▲  2.6002X10'6 
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▼  6. 0163xl0'8 
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Figure  4.24  Total  displacement  of  simulated  PMN-30%PT  k3i  bar  at 
fundamental  resonant  frequency  38.750  KHz 
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freq(382)=1.1525e5  Surface:  Total  displacement  (m) 


5 


▼  1.568xl0‘7 


Figure  4.25  Total  displacement  of  simulated  PMN-30%PT  k3i  bar  at  2nd 
harmonic  of  resonant  frequency  1 15.250  KHz 
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freq(680)=1.8975e5  Surface:  Total  displacement  (m) 
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Figure  4.26  Total  displacement  of  simulated  PMN-30%PT  k^i  bar  at  3rd 
harmonic  of  resonant  frequency  189.750  KHz 

4.4.1  Summary 

The  plots  in  Figure  4.24,  Figure  4.25,  and  Figure  4.26  are  each  dual  graphs  of  the  physical  shape 
change  and  the  total  displacement.  The  physical  shape  is  an  exaggerated  representation  of  the 
sample  geometry  under  sinusoidal  modulation  at  the  specified  frequency.  The  total  displacement 
is  indicated  the  by  the  color  on  the  color  legend  and  is  in  meters.  The  total  displacement  is  the 
distance  from  any  given  point  on  the  non-modulated  sample  to  the  new  point  with  the  applied 
signal.  It  is  clear  from  the  figures  that  the  modes  are  longitudinal  and  that  higher  frequencies  will 
increase  the  nodal  points.  Using  this  model  in  COMSOL  multiphysics  results  in  reliable 
prediction  of  the  resonate  modes. 
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4.5  Time  Domain  Piezoelectric  model 


In  order  to  model  AC  biasing  of  the  EOM  and  get  a  full  view  of  its  response  it  is 
necessary  to  model  the  system  in  the  time  domain  which  brings  its  own  set  of  challenges.  This 
model  uses  only  the  COMSOL  piezoelectric  application  mode  to  determine  the  phase  shift  due  to 
the  induced  Birefringence  of  the  device.  The  optical  propagation  through  the  device  is  not 
modeled  here  because  of  the  extremely  large  computational  cost.  The  model  is  setup  as 
illustrated  in  the  below  figure.  The  sample  modeled  in  PZN-4.5%PT  with  electroded  on  the  xy 
faces  with  the  assumed  optical  path  along  the  y  axis.  Two  sinusoidal  voltage  sources  are  used. 
Vsignal  is  swept  over  a  range  of  frequencies  but  Vbais  is  set  to  one  constant  frequency. 


Figure  4.15  Schematic  setup  of  dual  frequency  (piezo-resonant  ac  biased) 
electro-optic  modulator  simulation 
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The  refractive  index  changes  are  calculated  in  the  same  manner  as  the  frequency  domain.  The 
phase  change  due  to  the  induced  birefringence  is  calculated  by  first  finding  the  average  refractive 
index  in  the  z  direction  along  the  optical  path  through  the  very  center  of  the  sample.  This  average 
index  value  can  then  be  compared  to  the  results  with  no  applied  voltage  to  find  the  induced  phase 
shift. 


Birefringence  due  to  primary  and  secondary 
electro-optic  effects  and  AC  biasing 
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Figure  4.16  Effect  of  piezo-resonant  ac  signal  biasing  on  Birefringence  of 
electro-optic  modulator 

The  near  the  resonant  frequency  of  the  sample  a  beat  frequency  is  often  apparent  and  in  the  form 
of  an  additional  amplitude  modulation.  The  frequency  of  the  beat  is  a  function  of  the  modulation 
voltage  and  the  resonant  frequency  of  the  sample. 
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(4.17) 


\fbeat  I  —  I  f modulation  fresonancel 

Because  of  fbeat  there  is  no  steady-state  response  but  there  is  a  constant  envelope  steady-state 
response  (CESS).  This  CESS  refers  to  the  fact  that  the  amplitude  of  the  envelope  is  sinusoidal 
and  predictable  and  it  is  only  necessary  to  simulate  to  the  point  where  the  envelope  is  maximum 
in  order  to  find  the  full  response.  Furthermore  because  as  the  modulation  gets  close  to  the 
resonance  frequency  the  bean  frequency  becomes  smaller  thus  more  cycles  are  required  to  reach 
the  CESS.  When  modulating  at  the  exact  frequency  the  real  CESS  would  actually  take  an  infinite 
about  of  time  to  determine  but  can  still  be  accurately  estimated. 


lim  CESS  =  oo 

f  mocL^fres 
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Figure  4.17  FEA  time  domain  simulation  of  electro-optic  modulator  showing 
additional  modulation  near  resonance 


This  extra  amplitude  modulation  on  top  of  the  phase  modulated  signal  is  a  well-known 
disadvantage  of  using  modulation  near  the  natural  piezoelectric  resonance.  This  effect  however 
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could  possibly  be  compensated  for  through  either  real  time  processing  or  an  alternative  optical 
configuration.  Assuming  a  posteriori  knowledge  of  the  device  resonant  modes  it  real  time  FPGA 
processing  of  the  signal  could  be  done  to  correct  the  modulation  amplitude.  Clearly  the  unwanted 
modulation  corresponds  to  the  resonant  frequency  as  demonstrated  above  thus  a  foreknowledge 
of  the  resonant  modes  could  be  used  to  amplitude  correct  using  high  speed  electronics. 
Alternatively  this  correction  could  also  be  using  two  identical  crystals  as  shown  in  Figure  4.18. 
A  common  light  source  is  incident  upon  both  crystals  and  the  driving  modulation  signal  is 
directly  connect6ed  to  one  crystal  while  the  other  is  delayed  by  rc/2  before  reaching  the  seconds 
crystal.  When  the  modulated  signals  are  again  combined  into  a  common  light  signal  the  resultant 
signal  will  be  the  addition  of  the  signals  giving  a  corrected  output. 


Figure  4. 18  Proposed  method  of  correction  for  unintended  amplitude 
modulation  near  sample  resonance 
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4.6  Comparison  of  Bulk  EOM’s 

There  are  many  types  of  electro-optic  modulators  configured  for  modulation  of 
amplitude,  phase,  frequency,  etc.  and  each  can  also  be  configured  as  a  non-resonant  of  resonant 
modulator.  Resonant  modulators  typically  use  an  inductor  in  series  with  the  mostly  capacitive 
parallel  plate  configuration  of  the  EO  device.  The  addition  of  this  component  allows  for  a  higher 
voltage  than  that  applied  by  the  frequency  source  at  one  specific  resonant  frequency,  effectively 
lowering  the  half  wave  voltage  of  the  device.  Typically  these  devices  (both  resonant  and  non¬ 
resonant)  are  operated  far  from  the  natural  mechanical  resonance  so  are  to  avoid 

A  distinction  must  be  made  here  between  the  electrically  resonant  modulators  and  the 
mechanically  resonant  modulator  discussed  in  this  thesis.  Electrically  resonant  modulators  use 
additional  electrical  components  to  form  and  electrical  tank  circuit  resulting  in  a  low  half  wave 
voltage  at  one  frequency  over  a  very  small  bandwidth  of  operation  (high  Q).  Mechanically 
resonant  modulators  use  high  strains  resultant  from  the  natural  acoustic  vibration  mode  of  the 
material  to  give  higher  phase  changes  effectively  lowering  the  half  wave  voltage.  The  drawback 
of  these  devices  is  an  often  unwanted  amplitude  modulation  (in  addition  to  the  desired  phase 
modulation).  The  addition  of  AC  biasing  allows  for  broadening  the  effective  usually  bandwidth 
of  these  devices. 


Table  4.5  Qualitative  comparison  of  half-wave  voltage  and  bandwidth  of  common  electro-optic 
modulator  types _ _ _ _ 


Type 

V, 

Bandwidth 

comments 

Non-resonant 

high 

high 

Electrically  resonant 

low 

low 

Mechanically  resonant  (no  AC  bias) 

low 

low 

Unwanted  amplitude 
modulation 

Mechanically  resonant  (with  AC  bias) 

low 

high 
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The  series  resonant  circuit  depicted  in  Figure  4.19  consists  of  a  voltage  driver  (VI),  a  resistor  R1 
which  represent  the  transmission  line,  a  resistor  R2  witch  is  a  terminator  designed  to  minimize 
reflections,  an  inductor  LI  and  a  capacitor  Cl  which  models  the  actual  EOM  parallel  plate 
configuration. 


R1 


R2 


LI 


.ac  oct  10  100  100e6 


Figure  4.19  Series  resonant  electro-optic  modulator  circuit 


The  capacitance  of  the  device  is  Cl=(so*sr*As)/d  and  the  inductance  of  LI  is  chosen  to  give 
equal  but  opposite  reactance  at  a  specific  resonant  frequency. 
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4  2 

Assuming  an  EOM  of  thickness  d=10[mm]  and  length=25[mm]  the  area  As=2.25xl0‘  m  and 
sr=29.16  gives  Cl=6.45pF.  For  a  resonate  frequency  at  250  KHz,  Ll=62.83  mH.  Using  these 
parameters  a  1  volt  source  voltage  can  easily  give  841  volts  at  the  sample.  The  figure  below 
shows  the  response  of  the  circuit  in  Decibels;  clearly  there  is  a  very  large  signal  at  resonance 
over  a  very  narrow  band. 


Figure  4.20  Voltage  at  Sample  (in  dB)  resultant  from  series  resonant  circuit 


4.6.1  Finite  Element  Model  Comparing  Modulator  Types 

To  compare  these  modulator  types  a  finite  element  model  was  made  for  each.  Staring 
with  the  formulas  for  calculating  the  Pockels  effect  for  3m  symmetry  materials  as  derived  in 
Chapter  3. 

As  derived  in  chapter  1  for  3m  Symmetry  and  assuming  Ei  =  E2  =  0 

1 

n'x  ~  no  -  2  nor13Ez 


n'y  ~  n0  -  -n30r?3Ez 

1 

n'z~ne-  ~n3er*3Ez 


(4.23) 
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The  Birefringence  is 
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(4.24) 


The  phase  for  a  transverse  modulator  of  thickness  t  and  Optical  propagation  length  L: 
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The  half  wave  voltage  is  the  point  when  induced  phase  change  =  n 
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(4.31) 


To  induce  photoelastic  effects  the  constant  strain  coefficient  r~  can  be  replaced  by  the  constant 
stress  coefficient  r~  which  include  the  constant  strain  coefficient,  elasto-optic  and  piezoelectric 


coefficients. 

ri3  =  ri3  +  d-3iPn  +  d31p12  +  d33p13 
r33  =  V13  +  2d31p31  +  d33p33 


(4.32) 

(4.33) 
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(4.34) 


Vn  = 
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Alternatively  when  fonning  the  equations  for  the  FEA  model  slightly  less  simplified  method  can 
be  used.  For  3m  Symmetry  and  assuming  Ei  =  E2  =  0 


—f  +  rnE 3  +  PuXU  +  PuX  22  +  P\3X33  +  P\4X23 


+  ruE}  +  pnxu  +  pnx22  +  px 3X33  -  pux23 

/  n -  (4-35) 
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«31  =  l/  'jP44X31+P4lXl2 

nn  =  l/  yl  Pl4X3\+  P66X\2 


The  total  phase  shift 
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The  induced  phase  shift  can  be  found  by  subtracting  out  the  natural  Birefringence 

A(p  =  —r(-(n'z  W  ~  n*M)  -  (ne  [d]  -  n0[A])) 


At  14  Acp  =  n 
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(4.37) 
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(4.39) 


The  perturbed  refractive  index  n'  can  be  found  in  the  model  by  taking  the  average  refractive 
index  through  the  center  of  the  sample  in  the  Y  direction.  The  refractive  index  with  no  electric 
field  applied  (n[A])  is  a  function  of  the  wavelength,  the  data  used  is  from  the  handbook  of 
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Optical  materials  [9].  This  particular  geometry  has  a  piezoelectric  resonance  at  about  72  KHz  but 
no  resonant  effects  at  60  KHz  so  a  non-resonant  modulator  and  mechanically  resonant  modulator 
can  be  easily  compared  simply  by  modulating  on  and  off  resonance.  The  linear  electro-optic 
coefficient  can  be  assumed  to  be  constant  because  is  changes  little  with  wavelength  [9]. 


Comparison  of  EO  modulators 


Figure  4.21  Comparison  of  the  half-wave  voltage  of  non-resonant,  electrically-resonant, 
and  piezo-resonant  electro-optic  phase  modulators 


These  results  show  that  the  COMSOL  simulation  matches  equation  very  well  confinning  the 
validity  of  the  model.  The  model  accurately  calculates  phase  modulator  half  wave  voltage.  The 
two  resonant  EOMs  give  significantly  improved  VV  In  this  case  Piezo-resonant  performs  better 
than  the  electrically  resonant  modulator,  however  it  is  possible  to  observe  the  opposite  and  the 
results  are  largely  dependent  on  the  chosen  piezoelectric  resonant  mode  used  and  the 
piezoelectric  coefficient.  These  results  demonstrate  improvement  in  half  wave  voltage  for  Piezo¬ 
resonant  modulators  and  significant  advantage  over  electrically  resonant  modulators.  It  is  also 
possible  to  combine  mechanically  and  electrically  resonant  to  get  very  low  half  wave  voltage, 
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this  would  over  limit  the  wide  bandwidth  potential  of  the  device.  It  may  be  that  this  sensitivity 
will  one  day  allow  for  modulation  directly  from  wireless  signals  without  the  need  for  active 
electrical  components. 

4.7  Summary 

Several  Finite  Element  models  were  constructed  to  study  electro-optic  enhancement,  to 
more  fully  understand  experimental  observations  and  to  examine  methods  that  would  be  difficult 
to  realize  experimentally.  The  model  of  the  combined  interaction  of  a  piezoelectric  sample  with 
optical  plane  wave  propagation  has  been  shown  to  give  results  consistent  with  experimental 
observations  of  piezoelectric  induced  enhancement.  The  model  in  section  4.5  used  to  compare 
electro-optic  modulator  types  produces  a  response  that  matches  the  electro-optic  equations. 
Additionally  this  model  clearly  demonstrates  the  usefulness  of  piezo-resonant  enhancement  in 
terms  of  low  half  wave  voltage.  The  time  domain  model  of  dual  signal  ac  biased  configuration 
was  able  to  show  wide  bandwidth  enhancement  when  biased  at  the  resonant  frequency  and 
matches  well  with  experimental  observations  from  [6].  Microwave  transmission  in  periodically 
poled  structures  was  modeled  in  section  4.2.  Poling  can  have  a  significant  effect  on  transmission 
particularly  when  the  poling  direction  is  transverse  to  the  wave  propagation.  The  transmission 
spectrum  can  be  tuned  by  adjusting  the  number  of  polled  elements  in  a  structure. 
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CHAPTER  5:  MECHANICAL  VIBRATION  MEASUREMENTS 
5.1  Laser  Vibrometry 

Laser  vibrometers  are  used  to  make  single  point  physical  displacement  measurements  and  can 
sense  these  movements  at  high  frequencies.  They  can  also  sense  velocity  and  acceleration.  The 
two  main  principles  involved  are  the  Doppler  Effect  and  interferometry.  The  Doppler  Effect  is 
the  observed  shift  in  frequency  of  a  reflected  wave  due  to  the  velocity  of  the  sensed  object. 
Vibrometers  measure  velocity  of  the  object  by  determining  the  frequency  shift  with 
interferometry.  The  PolyTec  UHF-120  vibrometer  used  for  the  research  in  this  dissertation  uses  a 
heterodyne  interferometer  (See  Figure  5.1). 


Figure  5.1  Photo  of  PolyTec  UHF-120  Vibrometer  and  measurement  fixture 
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Heterodyne  interferometry  uses  two  laser  frequencies  as  opposed  to  Homodyne  which  uses  one. 
As  illustrated  in  Figure  5.2  both  laser  signals  start  from  the  same  location;  they  are  both  linearly 
polarized  with  polarization  orthogonal  to  each  other  [21].  Both  signals  reach  a  polarizing  beam 
splitter  which  directs  the  F2  signal  over  to  a  set  of  mirrors  that  steer  the  beam  back  left  to  the 
splitter  to  the  direction  of  the  source.  The  Fi  signal  passes  straight  through  the  beam  splitter  and 
is  incident  on  the  sample  being  tested.  After  experiencing  a  Doppler  shift  from  the  sample 
vibration  the  signal  is  frequency  becomes  Fi  +/-  5F  and  is  reflected  back  towards  the  direction  of 
the  source.  The  detector  and  control  electronics  are  used  to  measure  the  beat  frequency  that 
results  from  the  F2  reference  signal  and  Doppler  adjusted  Fi  +/-  5F  to  find  the  velocity  and  the 
direction  of  displacement. 
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Figure  5.2  Heterodyne  Interferometry  configuration 

The  measured  frequency  shift  fD  is  twice  the  velocity  v  divided  by  the  wavelength  of  the  emitted 
wave  X. 


(5.1) 
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The  phase  shift  of  the  interfering  waves  A<p  is  determined  by  the  path  length  difference  between 
the  waves  and  is  proportional  to  the  number  of  interference  fringes  N.  This  is  how  the 
displacement  is  determined. 


4tt  d. 

A</>  =  =  2nN 

A 


(5.2) 


The  total  intensity  at  the  photodetector  is  determined  by  the  intensity  of  the  two  individual  beams 
/j  and  1 2  the  wavelength  and  the  path  difference  between  the  beams  rx  —  r2.  The  total  intensity 
may  vary  from  zero  to  four  times  one  single  intensity. 

hot  =  h  +  h  +  2y[Tj2Cos(2n(r1  -  r2)/X)  (5.3) 


A  photodetector  is  used  to  pick  up  the  optical  signal  from  the  interferometer.  The  current  from 
the  detector  is  shown  in  Equation  (5.4)  and  has  a  DC  component  iDC  and  AC  component  f  [22], 
The  AC  component  is  governed  by  the  Bragg  frequency  fB,  the  modulation  phase  angle  cpm  and 
the  zero-phase  angle  90.  The  displacement  of  the  sample  s(t)  will  determine  the  modulated  phase 
angle  as  in  Equation  (5.5).  The  derivative  of  both  sides  of  this  equation  using  the  relationships  in 
Equations  (5.6)  and  (5.7)  and  the  results  is  Equation  (5.8).  This  represents  the  Doppler  frequency 

idetOO  =  be  +  iCos(2nfBt  +  (pm  +  90)  (5.4) 
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The  Vibrometer  can  handle  a  range  of  scanned  input  frequencies.  The  device  records  the 
signal  from  the  source  and  measurement  head  in  the  time  domain  for  a  time  period  determined 
by  the  user.  This  time  period  must  be  long  enough  to  capture  all  of  the  scan  frequencies.  The 
vibrometer  software  then  does  a  transformation  of  the  signal  to  the  time  domain.  Because  of  the 
leakage  effects  of  the  transformation  used  the  resulting  amplitude  of  the  spectrum  may  not 
represent  the  actual  signal.  The  software  uses  the  Fast  Fourier  Transform  (FFT)  for  converting  to 
the  frequency  domain.  The  FFT  uses  a  discrete  spectrum  meaning  that  if  a  time  signal  frequency 
does  not  match  the  FFT  line  it  will  be  distributed  over  several  lines.  There  are  various 
windowing  techniques  used  to  reduce  to  reduce  these  affects.  A  window  function  is  multiplied 
by  the  time  domain  signal;  most  window  functions  start  and  stop  at  zero  in  order  to  reduce  jumps 
in  the  signal  near  the  time  window  edge  [23].  For  this  reason  the  ratio  of  the  source  and 
vibrometer  signals  are  used  in  this  dissertation.  Thus  for  example  instead  of  displacement  in 
meters,  the  ratio  meters/volt  is  used. 

The  enhancement  of  the  electro-optic  signal  is  in  part  results  of  the  high  strained  inside 
the  sample  during  resonance.  These  strains,  physical  displacement,  velocity  and  acceleration  can 
all  be  examined  and  understood  more  fully  by  use  of  laser  vibrometry.  A  high  frequency 
vibrometer,  the  Polytec  UHF- Vibrometer  was  used  to  test  the  mechanical  deformation  of  several 
electro-optic  materials  under  sinusoidal  agitation.  This  vibrometer  offers  only  single  point 
measurements  but  can  scan  the  response  of  the  crystal  over  a  range  of  frequencies.  In  order  to 
give  a  full  view  of  the  sample  surface  the  measurement  was  repeated  over  each  sample  in  a  grid 
patterned  layout. 
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5.2  Measurement  Setup 

The  vibrometer  measurement  setup  is  shown  in  Figure  5.3  below.  The  vibrometer  signal 
is  digitized  by  an  oscilloscope  that  is  controlled  by  a  computer.  The  sample  is  excited  by  an 
arbitrary  function  generator,  the  output  of  the  source  is  split  by  an  RF  splitter,  one  end  goes  to 
the  sample  and  the  other  to  one  channel  of  the  scope.  This  allows  the  software  to  measure  the 
actual  signal  applied  to  the  sample  and  the  resulting  vibration  signal  form  the  sample  vibration. 
The  Sync  output  of  the  source  which  is  a  TTL  level  signal  corresponding  to  the  frequency  at  the 
Source  output  is  feed  into  the  scope  so  that  the  excitation  frequency  can  be  measured.  The 
Polytec  VibroSoft  software  is  used  to  analyze  the  source  and  vibrometer  signals.  A  custom 
Visual  Basic  for  Applications  (VBA)  program  was  written  to  control  the  entire  vibrometer 
system.  The  program  accepts  the  user  input  for  the  grid  to  be  measured  and  the  automatically 
prompts  the  user  for  where  to  move  the  positioning  stage  and  then  controls  the  vibrometer 
software  and  then  saves  the  displacement  data  to  a  common  file.  Additionally  a  link  is  made  to 
MatLab  for  simultaneous  plotting  of  the  sample  surface.  Alternatively  a  program  was  also 
written  to  sweep  through  a  range  of  frequencies  and  collect  time  domain  data  at  a  single 
measuring  point.  The  program  does  this  by  setting  the  Source  frequency  through  GPIB  control. 
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O-Scope 


Figure  5.3  Vibrometer  measurement  configuration 

5.3  Samples  and  Configuration 

To  build  up  a  representation  of  the  sample  surface  a  2D  raster  scan  of  the  surface  is  made.  The 
figure  below  shows  a  representation  of  the  Quartz  and  PMN-PT  samples  that  were  measured. 
The  samples  were  both  constrained  on  the  same  surfaces  were  they  are  electroded,  using  a 
special  fixture  to  hold  them  from  one  end  to  allow  for  overall  vibration  of  the  sample.  The  quartz 
sample  was  measured  over  a  total  of  1 16  points  and  the  PMN-PT  sample  which  is  smaller  was 
measured  over  30  points.  The  vibrometer  is  able  to  measure  and  analyze  surface  vibration  is  a 
variety  of  ways.  Typically  the  source  generator  makes  a  scan  over  a  range  of  frequencies  and  the 
resulting  signal  from  the  device  is  converted  to  the  frequency  domain  with  a  Fast  Fourier 
Transform  (FFT).  When  a  cyclic  signal  source  is  used  the  sample  deformation  naturally  changes 
from  expansion  to  contraction  at  some  constant  offset  phase  from  the  source.  Given  this  response 
the  FFT  takes  the  maximum  absolute  value  of  this  displacement  thus  on  any  surface  plot  of  the 
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displacement  is  not  a  true  plot  of  the  displacement  at  a  given  point  in  time  but  of  the  maximum 
displacement  at  over  the  whole  time  cycle.  This  is  illustrated  Figure  5.5  in  A  where  is  it  shown 
that  the  sample  surface  is  not  constant  but  changes  with  time;  part  B  shows  that  the  vibrometer 
displacement  data  would  record  this  vibration  as  the  envelope  of  the  shift. 


Quartz 

2.48  x  5.64  x  18.49  mm 


PMN-30%PT  fOOll 

1.25  x  2.25  x  3  mm 


Figure  5.4  Multipoint  measurement  for  two  dimensional  representation  of 
sample  surface  using  vibrometer  test  data 
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Surface  geometry  is 
exp  ending-  c  ontrac  ting 


Yibrometer  measures 
maximum  displacement 
at  any  given  point 


B 


Figure  5.5  Surface  displacement  measurement  by  vibrometer  gives  amplitude 
of  displacement  while  actual  sample  is  periodically  changing  as 
per  the  applied  sinusoidal  field 

5.4  Surface  Plot  Results 

The  surface  displacement  of  Quart  and  PMN-30%PT  were  measured  over  a  range  of 
frequencies.  The  displacement  of  these  samples  was  also  simulated  under  identical  conditions 
using  Finite  Element  Modeling.  In  the  plots  below  the  yellow  section  indicates  the  measured  area 
with  displacement  at  zero.  The  black  section  is  the  clamped  portion  of  the  sample  and  the  brown 
is  both  unmeasured  and  unclamped.  The  Figure  below  indicates  the  Quartz  surface  displacement 
at  the  strongest  resonance  indicating  good  agreement  between  simulation  and  experiment. 
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Figure  5.6  Quartz  (sample  #2)  comparison  of  COMSOL  Simulation  and 
vibrometer  surface  displacement 

The  PMN-30%PT  sample  was  measured  both  in  the  time  domain  and  frequency  domain. 
The  frequency  domain  resulting  are  show  in  the  next  section.  The  time  domain  results  are  shown 
here  at  329  KHz  and  the  displacement  is  given  at  one  specific  time.  The  displacement  has  two 
main  factors,  a  purely  up  and  down  z  direction  motion  and  partly  a  twisting.  The  measurement 
results  match  well  with  the  COMSOL  eigenfrequency  solver  at  368  KHz  that  demonstrates  the 
good  agreement  for  the  resonate  mode. 
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PMN-30%PT  Simulation  vs.  Experiment 


Comsol  Simulation 
368  KHz 


Vibrometer  experiment 
329  KHz 


Figure  5.7  PMN-30%PT  (sample  #1)  COMSOL  simulation  and  vibrometer 
surface  displacement 

In  Figure  5.8  the  PMN-PT  sample  plots  are  shown  at  the  major  resonant  frequencies,  all 
plots  are  on  the  same  scale.  Also  included  in  the  figure  for  comparison  are  the  admittance  and 
electro-optic  measurement  results.  It  is  clear  from  the  results  that  the  largest  sample  displacement 
does  not  always  match  the  largest  electro-optic  response.  The  excitation  frequency  of  328750  Hz 
has  a  EO  response  of  about  10.4  uVrms  while  the  relatively  small  displacement  at  891875  Hz 
gives  the  largest  EO  response. 
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Figure  5.8  PMN-30%PT  displacement  plots  at  major  resonant  frequencies 

5.5  High  Frequency  Displacement 

Typically  photoelastic  effects  are  only  apparent  up  to  the  low  MHz  region  however 
because  it  is  believed  that  piezoresonant  displacement  currents  can  couple  with  transverse  optical 
phonons  (which  can  exists  up  to  GHz)  it  may  be  the  enhancement  can  be  much  higher  than 
previously  measured.  Thus  vibrometer  experiments  have  been  perfonned  in  to  the  high  MHz 
rejoin  to  find  the  typical  surface  vibrations  of  crystals  in  this  range  and  how  this  response  may 
effect  enhancement. 

PMN-30%PT  d33=198 1  [C/N],  d3i=  d32=-921  [C/N]  (5.4) 

[C/N]  =  [m/V]  (5.3) 


The  PMN-30%PT  sample  was  tested  both  for  its  electro-optic  response  and  for  its 

vibration.  The  vibration  was  tested  along  its  major  X,  Y,  Z  axis  directions.  Figure  5.9  shows  the 

testing  configurations;  silver  paste  electrodes  were  used  on  the  XY  surfaces  allowing  for  electric 

92 


field  in  the  Z  direction  (also  the  poling  direction).  The  sample  was  rotated  in  three  different 
directions  to  measure  the  d3 1,  d.32  and  d33  piezoelectric  coupling.  The  vibrometer  is  suited  to  high 
frequency  vibration  measurements  while  piezoelectric  coefficients  are  typically  measured  at  very 
low  frequencies.  Piezoelectric  coefficients  are  reposted  in  terms  of  Coulombs  per  Newton  (C/N) 
or  meters  per  Volt  (mV)  which  are  equivalent  units.  PMN-30%PT  has  an  identical  d3i  and  d32 
coefficient  of  -0.921  [nm/V]  and  a  d33  of  1.981  [nm/V]  [19].  The  coefficients  cannot  be  expected 
to  be  exactly  the  same  when  measured  at  high  frequencies  but  it  is  clear  from  the  results  of 
Figure  5.10  that  the  measurement  is  very  consistent  with  reported  values.  Note  that  this  particular 
vibrometer  test  only  considers  the  magnitude  of  the  response,  not  the  phase,  thus  it  is  the 
absolute  value  to  the  piezoelectric  coefficient  that  is  measured. 
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Figure  5.9  Configuration  for  PMN-30%PT  sample  for  electro-optic  and 
vibrometer  tests 


The  results  of  this  test  are  shown  in  Figure  5.9  and  there  is  clearly  a  large  increase  in  the  electro¬ 
optic  signal  near  the  resonant  frequencies  as  indicated  by  the  vibrometer  results.  This  correlation 
is  particularly  true  for  the  fundamental  modes  below  1MHz  and  is  expected  due  to  the  high  strain 
coupled  to  the  refractive  index  through  the  photoelastic  effect.  Above  1MHz  the  crystal  becomes 
more  clamped  meaning  there  is  less  displacement  however  at  the  same  time  there  is  still  a  rather 
high  spike  in  the  electro-optic  response.  This  coupling  seems  to  be  due  to  more  than  just  the 
photoelastic  and  electro-optic  response  and  is  perhaps  a  result  of  the  induced  polarization  in  the 
material.  This  induced  polarization  has  less  of  frequency  limitation  as  compared  to 
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photoelasticity  but  can  still  result  in  the  sample  type  of  resonant  induced  improvement  of  the 
signal. 


PMN-30SPT  Electro  opoc  arc  v Beaton  response 


* 


Figure  5.10  PMN-30%PT  electro-optic  and  vibrometer  tests  indicating  that 
electro-optic  enhancement  persists  in  to  the  MHz  region  where 
sample  strain  is  minimal 


In  addition  to  the  single  point  measurements  a  full  view  of  the  surface  displacement  for  PMN- 
30%PT  is  presented  in  Figure  5.11.  In  the  figure  the  red  block  represent  the  sample  at  rest  with 


no  displacement,  the  blue  dots  are  the  actual  measurement  points  on  the  sample  and  the  other 
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EO(Mttb) 


area  is  the  interpolation  grid  between  the  test  points.  Three  scan  tests  where  done  to  test  out-of¬ 
plane  displacement  in  the  x,  y,  and  z  directions  and  each  was  combined  in  to  one  figure  for 
primary  resonance  at  324375  Hz.  The  actual  displacement  was  multiplied  by  10  to  give  an 
agitated  view  of  the  sinusoidally  agitated  geometry.  The  results  show  that  the  largest  change  in 
sample  dimension  is  in  the  z  direction  (same  as  the  applied  field). 
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Figure  5.11  Full  view  of  PMN-30%PT  sample  surface  displacement 
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5.6  Velocity  and  Acceleration 

The  enhanced  optical  signal  is  due  to  more  than  just  the  static  strain  in  the  sample  but 
rather  the  dynamic  nature  of  the  sample  under  resonance.  Because  of  this  relationship  it  may  be 
interesting  to  determine  the  association  of  the  surface  velocity  and  acceleration  under  applied 
field  and  association  with  the  optical  signal  enhancement.  A  single  point  measurement  was 
performed  on  the  PMN-30%PT  sample  and  the  displacement,  velocity  and  acceleration  were 
recorded  (See  Figure  5.12).  The  electric  field  was  applied  in  the  z  direction  and  the  measured 
parameter  was  in  the  same  direction.  The  actual  measurement  point  was  near  the  center  of  the 
sample  on  the  surface.  The  all  three  measurements  display  peak  near  the  resonant  frequencies 
particularly  near  330,  570  and  880  KHz;  all  three  measurements  also  have  lower  amplitude  at 
each  of  these  successive  frequencies.  The  acceleration  results  however  do  stand  out,  in  that  they 
not  decrease  nearly  as  much  as  the  displacement  results  do. 
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Vibrometer  measurement  of  PMN-30%PT 
Field  applied  in  z  direction  and  measured  in  z  direction 


m/s/V  m/s  /V 


■  4x10 


-3x10 


-2x10 


1x10 


Figure  5.12  Vibrometer  measurement  of  displacement,  velocity,  and 
acceleration  of  a  PMN-30%PT  sample 


5.7  Summary 

The  physical  dimension  change  of  the  piezoelectric  samples  was  studied  with  an  ultra- 
frequency  vibrometer.  Point-by-point  out-of-plane  data  was  reconstructed  and  interpolated  to 
give  a  representation  of  the  sample  surface  under  sinusoidal  agitation.  The  ratio  of  applied  signal 
divided  by  the  displacement  in  the  measured  direction  was  recorded  effectively  determining  the 
high  frequency  piezoelectric  coefficient.  The  piezoelectric  coefficient  measured  is  in  general 
smaller  than  that  reported  in  literature  because  the  crystal  cannot  respond  fast  enough  to  the  field 
however  at  resonance  the  coupling  is  much  higher  and  comparable  to  the  low  frequency  values. 
The  velocity  and  acceleration  of  the  sample  surface  was  also  examined  and  it  was  observed  that 
the  acceleration  of  the  surface  does  not  decrease  as  much  as  the  displacement  when  compared  at 
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major  resonate  frequencies.  These  results  suggest  that  the  acceleration  have  a  role  in  high 
frequency  enhancement  given  the  both  the  acceleration  and  electro-optic  signal  may  show 
significant  peaks  when  evaluated  at  higher  frequencies. 
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CHAPTER  6:  ACOUSTO-OPTIC  FIGURE  OR  MERIT  OPTIMIZATION 


Section  6  McIntosh,  R.,  A.  S.  Bhalla,  et  al.  (2012).  " Finite  element  modeling  of  acousto-optic 
effect  and  optimization  of  the  figure  of  merit. "  849703-849703. 

6.1  Introduction 

A  study  of  the  acousto-optic  (AO)  effect  in  a  family  of  oxide  crystals  (including  e.g., 
TiCF,  ZnO,  LiNb03,  and  ferroelectric  perovskites)  as  well  as  semiconductors  has  been  conducted 
by  finite  element  analysis  method.  In  addition,  the  acousto-optic  figure  of  merit  (FOM)  as  a 
function  of  material's  refractive  index,  density,  effective  AO  coefficient  and  the  velocity  of  the 
acoustic  wave  in  the  material,  is  also  investigated.  By  examining  the  directional  dependent 
velocity,  acousto-optic  coefficients,  and  refractive  index,  the  acousto-optic  FOM  can  be 
calculated  and  plotted  in  all  directions  revealing  the  optimal  crystal  orientation  to  maximize 
coupling  between  the  optical  and  acoustic  waves.  A  finite  element  model  was  developed  to 
corroborate  the  improved  interaction.  The  model  examines  the  diffraction  that  occurs  on  the 
optical  wave  as  it  travels  through  an  acousto-optic  medium.  The  combined  information  gained 
from  Mathematica  and  COMSOL  Multiphysics-based  modeling  is  shown  to  be  an  effective 
means  of  predicating  acousto-optic  device  functionality. 

There  are  two  types  of  photoelasticity,  the  piezo-optic  effect  and  the  elasto-optic  effect. 
Pressure  can  change  the  refractive  index  of  a  material  both  through  the  density  and  through  the 
polarizability.  When  a  material  in  compressed  the  atoms  move  closer  together  and  the  refractive 
index  increases  but  at  the  same  time  the  electrons  are  bound  more  tightly  which  reduces 
polarizability  and  reduces  the  refractive  index.  For  this  reason  and  because  the  two  effects  are  on 
the  same  order  of  magnitude,  pressure  can  increase  or  decrease  the  refractive  index  depending  on 
the  material  [24], 
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Piezo-optic  A Bi}  =  ni]klXkl 


(6.1) 


Elasto-optic  ABtJ  =  pijklxkl  (6.2) 

Where  is  known  as  piezo-optic  coefficient  that  is  defined  by  stress  (Xki)  induced  changes  in 

optical  indicatrix  A Similarly  pykiis  known  as  elasto-optic  coefficient  that  is  defined  by  strain 

(xki)  induced  changes  in  optical  indicatrix  A B^.  The  acousto-optic  figure  of  merit  (FOM)  is  a 

measure  of  the  suitability  of  a  material  to  modulate  the  diffraction  intensity.  The  refractive  index, 

elasto-optic  coefficient,  density  and  acoustic  wave  velocity  are  all  used  in  this  calculation  but  it 

is  the  refractive  index  and  acoustic  wave  velocity  that  are  the  dominant  factors.  The  slower  the 

acoustic  and  optical  waves  in  the  material  the  more  interaction  possible  [24],  There  are  other 

figures  of  merit  related  to  acousto-optic  devices  however  the  FOM  used  in  Equa.(6.3)  is  used 

primarily  for  gauging  the  power  efficiency  of  AO  materials  and  is  not  for  example  used  to 

determine  the  usable  bandwidth  of  a  the  device[25]. 

n6p2  (6.3) 

XI2  =  j- 
pv6 

6.2  Importance 

The  acousto-optic  effect  is  used  in  many  devices  including  modulators,  frequency 
shifters,  filters  and  beam  deflectors.  AO  modulators  are  used  in  many  nanosecond  pulsed  lasers 
for  Q-switching.  The  modulator  is  used  inside  the  laser  cavity  to  remove  the  cavity  feedback 
until  a  large  population  inversion  is  built  up.  Beam  deflecting  modulators  can  be  used  for  laser 
graphics  projectors,  optical  tweezers,  and  optical  switching.  Acoustic  devices  are  also  widely 
used  for  medical  imaging  due  to  the  excellent  penetration  depth  of  acoustic  waves. 


6.3  Recent  Status 
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Recent  research  in  acousto-optics  has  included  higher  frequency  acoustic  interaction  and 
metamaterials.  Though  there  are  many  types  of  metamaterials  the  designation  “meta”  is  typically 
used  to  refer  to  a  material  with  negative  permittivity  and  penneability.  This  configuration  is 
usually  created  by  some  type  of  resonance  or  periodic  condition  in  the  material  and  can  result  in 
the  group  velocity  of  a  wave  in  the  material  being  positive  while  the  phase  velocity  is 
negative[26].  Current  interest  in  acousto-optics  encompasses  a  wide  range  of  research  for  low 
power  consumption,  modeling  techniques,  biological  applications,  integration  with  photonic 
crystals,  integrated  optics,  and  MEMS,  etc.  It  is  highly  desirable  to  give  a  comprehensive  picture 
of  all  these  property  coefficients  in  relation  to  their  point  group  symmetries  and  do  so  by 
demonstrating  such  capability  in  some  important  and  well  established  materials. 

6.4  Approach 

The  necessary  parameters  for  calculating  the  figure  of  merit  of  several  common  materials 
are  listed  in  Table  6.1  through 

Table  6.6.  Table  6.1  and  Table  6.2  give  the  elasto-optic  coefficients  which  have  no  units. 
To  calculate  the  elastic  wave  velocity  in  linear  elastic  materials  (such  as  TiC>2  and  PbMoO/O  only 
the  elastic  stiffness  coefficients  Cjj  and  the  density  p  are  needed;  however  for  piezoelectric 
materials  consideration  of  the  piezoelectric  coupling  coefficients  e,j,  and  the  permittivities  are 
also  required.  The  material  coefficients  listed  are  from  sources  [9,  24,  27-30], 


Table  6.1  Elasto-optic  coefficients  -  of  unitless  -  Part  A 


Pn 

P 1 2 

P 1 3 

Pl4 

P 1 6 

P31 

LiNb03 

-0.026 

0.09 

0.133 

-0.075 

- 

0.179 

Ti02 

-0.001 

0.113 

-0.167 

- 

- 

0.106 

ZnO 

0.222 

0.0999 

-0.111 

- 

- 

0.0888 

PbMoC>4 

0.24 

0.24 

0.225 

- 

0.017 

0.175 

GaAs 

-0.165 

-0.14 

- 

- 

- 

- 

GaP 

-0.151 

-0.082 

- 

- 

- 

- 

Ge 

0.27 

0.235 

- 

- 

- 

- 
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Table  6.2  Elasto-optic  coefficients  -  of  unitless  -  Part  B 


P33 

P41 

P44 

P45 

P61 

P66 

Ref. 

LiNb03 

0.0171 

-0.151 

0.146 

- 

- 

- 

[24] 

Ti02 

-0.064 

- 

0.0095 

- 

- 

-0.066 

[24] 

ZnO 

-0.235 

- 

0.0585 

- 

- 

- 

[31] 

PbMoCE 

0.3 

- 

0.067 

-0.01 

0.013 

0.05 

[32] 

GaAs 

- 

- 

-0.072 

- 

- 

- 

[24] 

GaP 

- 

- 

-0.074 

- 

- 

- 

[24] 

Ge 

- 

- 

0.125 

- 

- 

- 

[30] 

Table  6.3  Elastic  coefficients  -of  units  of  1011  [N/m?] 


Cl  i 

Cl2 

Cl3 

C14 

Cl6 

C33 

C44 

C66 

Ref. 

LiNbOa 

2.03 

0.53 

0.75 

0.09 

- 

2.45 

0.6 

- 

[33] 

Ti02 

2.7143 

1.7796 

1.4957 

- 

- 

4.8395 

1.2443 

1.9477 

[34] 

ZnO 

2.097 

1.211 

1.051 

- 

- 

2.109 

0.4247 

- 

[35] 

PbMo04 

1.09 

0.68 

0.53 

- 

-0.14 

0.92 

0.267 

0.335 

[36] 

GaAs 

1.1877 

0.5372 

- 

- 

- 

- 

0.5944 

- 

[37] 

GaP 

1.412 

0.6253 

- 

- 

- 

- 

0.7047 

- 

[38] 

Ge 

1.2835 

0.4823 

- 

- 

- 

- 

0.6666 

- 

[39] 

2 

Table  6.4  Piezoelectric  coupling  coefficients  -  of  units  of  [C/m“] 


en 

en 

eis 

e22 

e3i 

e33 

Ref. 

LiNbOa 

- 

- 

4.1607 

2.442 

0.8661 

3.7 

[24] 

ZnO 

- 

- 

-0.352501 

- 

-0.35706 

1.56416 

[24] 

GaAs 

- 

0.154544 

- 

- 

- 

- 

[24] 

GaP 

- 

0.1 

- 

- 

- 

- 

[29] 

Table  6.5  Relative  dielectric  pennittivity  coefficients  -  of  unitless 


£11 

£33 

Ref. 

LiNbOa 

43.6 

29.16 

[28] 

ZnO 

8.5446 

10.204 

[28] 

GaAs 

12.459 

- 

[28] 

GaP 

11.1 

- 

[28] 
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Table  6.6  Refractive  index  coefficients  and  density  -  density  of  units  of  [kg/m3] 


nn 

n33 

Ref. 

P  [kg/m3] 

Ref. 

LiNb03 

2.232 

2.156 

[9] 

4644 

[9] 

Ti02 

2.584 

2.872 

[9] 

4260 

[9] 

ZnO 

2.015 

1.998 

[9] 

5606 

[9] 

PbMo04 

2.2584 

2.3812 

[9] 

6920 

[9] 

GaAs 

3.37 

- 

[9] 

5340 

[9] 

GaP 

3.35 

- 

[9] 

4130 

[9] 

Ge 

4.0 

- 

[30] 

5330 

[30] 

Typically  AO  materials  are  chosen  based  on  both  their  availability  and  their  figure  or 
merit.  Some  of  the  most  commonly  used  materials  are  Lithium  Niobate  (LiNbOs),  Lithium 
Tantalate  (LiTaOs),  Lead  Molybdate  (PbMo04).  paratellurite  (TeOo),  fused  silica  (SiOi),  Rutile 
(Ti02),  Zinc  Oxide  (ZnO),  and  Gallium  Arsenide  (GaAs).  Many  other  materials  are  also  used 
that  have  high  figures  of  merit.  This  paper  will  consider  Lithium  niobate,  Rutile,  Zinc  Oxide 
(Wurtzite),  Lead  Molybdate,  Gallium  Arsenide,  Gallium  Phosphide,  and  Germanium  due  to  the 
availability  of  data.  The  developed  capability  can  be  readily  extended  to  advanced  materials 
upon  the  availability  of  needed  physical  parameters. 

Plotting  of  the  directionally  dependent  material  properties  is  more  readily  done  in  the 
spherical  coordinate  system.  Using  spherical  coordinates  the  direction  ( Z[  )  relative  to  the 
Cartesian  system  of  Zi,  Z2,  and  Z3  is  the  angle  cp  from  the  Z\  axis  on  the  Z1-Z2  plane  and  the 
angle  theta  6  from  the  Z3  axis  as  depicted  in  Figure  6.1.  Additionally  the  radius  (i.e.  length  of  the 
Z[  vector)  represents  the  magnitude  of  the  effect. 


105 


Figure  6.1  Conversion  scheme  of  Cartesian  coordinate  system  to  spherical 
coordinates 


The  elasto-optic  effect  is  a  4th  rank  tensor  property  that  can  be  fully  described  with  a  6x6 
matrix  using  Voigt  notation.  To  analyze  longitudinal  directional  dependence  of  elasto-optic 
effect  (  p'iUi  )  about  an  arbitrary  direction  (/-direction),  the  6x6  “a-matrix”  is  used.  Other  effects 
such  as  the  refractive  index  can  be  evaluated  by  way  of  the  3x3  transformation  “a-matrix”.  Both 
the  "a"  and  “a”  matrices  are  made  up  of  the  directional  cosine  elements  ay  as  seen  in  the 
equations  below. 


a31  =  Sin#  Cos <p 
a12  =  Sin#  Sirup 
a13  =  Cos  # 
a2i  =  Cos#  Cos  <p 
a22  =  Cos#  Sin0 
a23  =  —Sin# 
a31  —  Sin  (p 
a32  •—  Cos (p 
a33  =  0 


(6.4) 
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all 

al2 

al3\ 

a21 

a22 

a23  I 

(6.5) 

a31 

a32 

a33/ 

'  ah 

ah 

ah 

2a12a13 

2al3all 

2a11a12 

ah 

ah 

ah 

2a.22a23 

2^23^21 

2^21^22 

ah 

ah 

ah 

2a32a33 

2a33a31 

2a31a32 

a2la31 

a22a32 

a23a33 

a22a33  +  a23a32 

a2ia33  +  a23a31 

a22a3i  +  a21a32 

a31all 

a32al2 

a33a\3 

a±2  a33  +  a13a32 

ai3a3i  +  aua33 

ana32  +  a12a31 

■alla21 

a\2a22 

a13a23 

ai2a23  +  a13a22 

a13a21  +  flUa2  3 

^11^22  +  ni2a21 

The  value  of  p'im  can  be  calculated  using  Equa.  (6.7),  by  multiplying  a  by  the  material  matrix  p 
(whose  components  are  defined  by  its  crystallographic  axis)  and  again  by  the  transpose  of  the  a 
matrix  and  extracting  the  top  left  element  (1,1). 

p'iiii(0 ’  (p)  =  a  P  aT  DXll  (6.7) 

In  the  following  text  the  direction  of  tensile  strain  is  designated  as  i-  (\  -)  direction  and  the 
physical  meaning  of  the  p[  X11  is  the  amplitude  of  longitudinal  effective  elastooptic  effect.  A  plot 
of  Pallid,  <p)  is  a  graphical  representation  of  its  directional  dependence,  p'mi  is  plotted  as  a 
function  of  theta  and  phi,  for  6  from  0  to  7r  and  <p  from  0  to  2n. 

The  refractive  index  (represented  as  a  3x3  matrix)  is  found  in  a  similar  manner  as  to  the  elasto¬ 
optic  effect  but  now  using  the  “a”  matrix. 

n'u{8 ,  (p)  —  a  n  aT  [1,1]  (6.8) 

The  sound  velocity  in  a  homogeneous  anisotropic  medium  has  multiple  solutions  for  any  given 
wave  propagation  orientation  and  thus  requires  a  different  method  of  evaluation  as  compared  to 
the  aforementioned  elasto-optic  and  refractive  index.  Typically  there  is  one  longitudinal  wave 
with  vibration  direction  parallel  to  propagation  and  two  transversal  shear  waves.  In  a  non- 
center-symmetric  material  each  of  the  wave  components  may  be  also  accompanied  by  oscillating 
polarizations  coupled  through  piezoelectric  effect.  Acousto-optic  materials  can  be  either  linear 


elastic  or  piezoelectric.  The  velocity  of  linear  elastic  materials  is  determined  by  orientation, 
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density,  and  the  elastic  constants,  piezoelectric  materials  however  must  also  consider  the 
piezoelectric  coupling  matrix  and  permittivity.  Only  in  a  few  circumstances  are  the  waves 
considered  to  be  pure  that  is  with  polarization  either  perfectly  parallel  or  perpendicular  to  the 
propagation  directions  and  are  more  typically  either  quasi-longitudinal  of  quasi-transverse 
especially  in  materials  of  low  symmetry. 

The  modified  Christoffel  equation  used  to  find  the  wave  velocities  is  shown  below.  The 
tensor  C*  represents  the  linear  elastic  portion  and  is  a  product  of  the  inverse  density  1/p,  the 
elastic  tensor  Cijki  and  the  directional  cosines  N;  Nj.  The  Q  Ck  tensors  and  constant  C  represent 
the  velocity  adjustment  resultant  from  the  piezoelectric  contribution.  emij  is  the  piezoelectric 
coupling  tenor  and  smj  the  permittivity.  It  is  important  to  note  that  s  here  is  not  the  relative 
permittivity  but  rather  the  absolute  permittivity  s  —  s0sr  in  F/m,  the  velocity  is  v,  Slk  is  the 
Kronecker  delta  ( 8ik  =  1  if  i  =  k  and  zero  otherwise ),  and  Uk  is  the  amplitude  of  the  lattice 
displacement  in  k-direction. 


9  \ 

+  c  v28ik)uk  =  0 

(6.9) 

C  —  pNmemjNj 

(6.10) 

1 

dik  ~  ~cijklNjNi 

P 

(6.11) 

Q  —  ^memij^l 

(6.12) 

Nt  =  SinO  Coscp 

(6.13) 

N2  =  Sind  Sincp 

(6.14) 

N3  =  Cosd 

(6.15) 

For  clarity  these  tensor  components  for  the  general  (triclinic)  case  are  written  explicitly  below 
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Cn  =-(c11/V12+c66/V22  +  c55JV32  +  2c56N2N3  +  2c15N3N1  +  2c16N1N2)  (6.16) 

P 

C22  =~(c66N12+c22N22  +  c44N32  +  2  c24N2N3  +  2c46N3N4  +  2c26N1N2)  (6.17) 

P 

C33  =  —{c33N1  +c44N2  +  c33N3  +  2c34N2N3  +  2c35N3N1  +  2c45N1N2)  (6.18) 

1 

623  =  ~(.C56^1  +C24^2  +  C34^3  +  (C23  +  C44)^2^3  +  (C36  +  C45)N3N1  +  (C25  +  C46)N1N2 )  (6.19) 

613  =  -{c^N-^+c^N-22  +  c35N32  +  (c36  +  c43)N2N3  +  (c13  +  c55)N3N4  +  (c14  +  c56)N1N2')  (6.20) 

P 

Cl2  ~  -(C16^12  +  C26^22  +  C45^32  +  (C25  +  C46)^2^3  +  (C14  +  C56)N3N1  +  (C12  +  C66)N1N2)  (6.21) 

P 


Ci  =  e11Al12+e26Af22  +  e35N32  +  (e25  +  e36)N2N3  +  (e15  +  e31)N3N1  +  (e16  +  e21)N1N2  (6.22) 

62  =  el6^1  +e22^2  +  e34^3  +  (e24  +  e32)^2^3  +  (el4  +  e36)N3N1  +  (e12  +  e26)N1N2  (6.23) 

63  =  ei5^i2+e24W22  +  e33N32  +  (e23  +  e24)N2N3  +  (e13  +  e35)N3N1  +  (e14  +  e25)N1N2  (6-24) 

C  =  p^nlV,2  +  £22JV22  +  £33JV32  +  2s23N2N3  +  2£13tV3lV1  +  2e12N1N2)  (6.25) 


The  method  of  calculation  is  more  easily  seen  by  converting  the  Christoffel  Equation  to  matrix 
form  and  inserting  the  above  relations  to  form  the  matrix  A  below.  The  velocities  are  now  simply 
the  square  root  of  the  eigenvalues  of  the  matrix  A  and  the  corresponding  wave  polarizations  are 
from  the  eigenvectors  of  A.  The  eigenvalues  of  A  always  exist  but  the  eigenvectors  may  not 
always  be  determinant  for  all  directions. 


Cn 

Cl2 

6"l3 

c- 12 

CiC2 

c3Ci 

A  = 

6*12 

6*22 

6*23 

+ 

C,C2 

r  2 
l2 

C2C3 

-6*23 

6l3 

6*33- 

-C3Ci 

c2Cs 

r  2 
l3 

(6.26) 


6.5  Results  and  Discussion 

Mathematica  has  been  used  to  provide  a  three  dimensional  visualization  of  the  elasto- 
optic  effect  for  a  crystal  of  any  point  group  symmetry.  One  can  select  the  point  group  and  then 
from  a  library  of  material  data  associated  with  that  group.  The  directionally  dependent  Figure  of 
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Merit  for  four  chosen  materials  is  shown  in  Figure  6.2  below.  All  results  in  the  figure  are  in 

15  3 

terms  of  10'  '  [s  /kg].  Note  that  the  negative  Cartesian  coordinates  are  for  reference  and  do  not 
indicate  a  negative  figure  of  merit 


M2(B,  cp)  = 


p  vL{e,(p )3 


(6.27) 
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Lithium  Niobate  -  Trigonal 


Rutile  -  Tetragonal 


* 


Zinc  Oxide  -Hexagonal 


Lead  Molybdate  -  Tetragonal 


Figure  6.2  Directional  dependence  of  elasto-optic  figure  of  merit  M2 

calculated  for  several  example  materials  representing  different 
crystal  systems 


Others  have  published  on  the  directional  dependence  of  the  elasto-optic  effect  [40]  and 
the  piezo-optic  effect  [41]  but  only  for  a  few  materials  and  not  all  have  considered  the  figure  of 
merit  directional  dependence.  Mathematica  was  used  to  calculate  the  figure  of  merit  by 
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combining  the  effects  of  density  elasto-optics,  acoustic  velocity  and  refractive  index  (Figure  6.2). 
It  is  typical  for  FOM  calculations  to  show  Mi  as  divided  by  the  FOM  of  isotopic  fused  silica 
(M2=  1.5 1x10'  [s  /kg])  thus  giving  a  unitless  relative  value. 

The  figure  of  merit  is  typically  used  to  compare  or  gauge  the  effectiveness  of  AO 
materials;  however  it  can  also  be  used  to  determine  the  ideal  material  orientation  for  maximum 
acoustic/optic  coupling.  The  elasto-optic  matrix  is  certainly  the  dominant  factor  in  detennining 
idea  coupling;  however  as  can  be  seen  in  the  FOM  calculation  the  refractive  index  and  the 
velocity  are  directionally  dependent  they  can  also  have  a  significant  impact  on  coupling.  The 
results  of  this  evaluation  are  consistent  with  those  from  Chang  [30]  with  relative  M2  of  GaAs 
(M2=69),  GaP  (M2=29.5),  Ge  (M2=482).  Our  results  for  GaAs,  GaP  and  Ge  are  respectively  64.3, 
30.5,  and  510.3  demonstrating  good  agreement  and  verification  of  the  methodology  developed. 
Figure  6.3  shows  the  FOM  of  GaAs  in  the  Z1-Z2  plane  and  a  close  match  to  the  literature  value 
reported  for  the  [110]  direction.  Furthennore,  this  work  permits  a  full  view  of  the  directional 
dependence  in  any  direction  interested,  a  powerful  feature  not  previously  available. 

The  three  dimensional  representations  (Figure  6.2)  can  be  used  to  effectively  determine 
the  maximum  coupling  orientation.  Using  this  information  the  ideal  material  orientations  have 
been  detennined  for  each  of  the  materials  and  summarized  in  Table  6.7.  The  values  are  listed  for 
0  <  8  <  n/2  and  0  <  cp  <  2n;  there  are  corresponding  maximum  orientations  for  8  >  n/2  only 
as  additional  directions  can  be  easily  generated  based  on  given  symmetry. 
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[010] 


Figure  6.3  Relative  acousto-optic  figure  of  merit  of  Gallium  Arsenide  in  ZpZi 
plane  and  comparison  to  the  literature  value  reported  in  the  [110] 
direction 
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Table  6.7  Maxima  of  Figure  of  Merit  for  each  material  with  associated  property 
in  the  same  direction.  Theta  and  phi  are  in  degrees 


Pii 

n 

v  [m/s]  p  [kg/m A  3] 

M2[1045  s3/kg] 

M2  relative 

(0^) 

LiNb03 

0.340 

2.203 

6487.6 

4644 

10.418 

6.899 

(51.7,90.0) 

(51.7,210.0) 

(51.7,330.0) 

Ti02 

-0.081 

2.782 

10079.5 

4260 

0.694 

0.46 

(34.1,0)  (34.1,90) 

(34.1,180) 

(34.1,270) 

ZnO 

-0.235 

1.998 

6133.5 

5606 

2.716 

1.799 

(0,0) 

0.222 

2.015 

6116.1 

5606 

2.572 

1.703 

(90, any) 

PbMo04 

0.300 

2.381 

3646.2 

6920 

48.909 

32.39 

(0,0) 

GaAs 

-0.244 

3.37 

5381.6 

5340 

105.066 

69.580 

(54.7,45) 

(54.7,135) 

(54.7,225) 

(54.7,315) 

GaP 

-0.204 

3.35 

6651.4 

4130 

48.242 

31.949 

(54.7,45) 

(54.7,135) 

(54.7,225) 

(54.7,315) 

Ge 

0.413 

4.0 

5543.5 

5330 

770.526 

510.282 

(54.7,45) 

(54.7,135) 

(54.7,225) 

(54.7,315) 

A  COMSOL  model  was  developed  to  examine  the  zero  order  diffraction  in  an  acousto¬ 
optic  modulator  and  the  improved  efficiency  that  can  result  from  properly  orienting  the  AO 
crystal  accurately  as  designed.  A  typical  configuration  for  an  AO  modulator  is  seen  in  Figure  6.4. 
The  incident  optical  beam  comes  in  from  the  top  medium  (air)  at  an  angle  alpha  to  the  normal  of 
the  surface  where  reflection,  refraction,  and  diffraction  occur.  A  periodic  strain  change  causes 
diffraction  of  the  beam.  The  acousto-optic  modulator  can  result  in  modulation  of  the  amplitude, 
the  deflection  angle  as  well  as  frequency  shifting  or  wavelength  filtering  [42].  The  setup  of  the 
piezoelectric  and  electromagnetic  domains  for  the  COMSOL  model  is  shown  in  Figure  6.5.  This 
model  uses  an  initial  strain  defined  in  the  mechanical  (piezoelectric)  portion  that  is  then  coupled 
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to  the  optical  model  (RF)  through  the  elasto-optic  effect  on  the  refractive  indices.  The  simulation 
took  Lithium  Niobate  as  an  example  which  has  3m  trigonal  symmetry.  A  periodically  changing 
strain  is  defined  in  the  AO  material;  the  strain  is  oriented  in  the  out-of-plane  direction  but 
changes  spatially  across  the  x  dimension.  Furthermore,  the  optical  polarization  is  in  the  out-of¬ 
plane  direction  meaning  that  this  is  a  configuration  corresponds  to  the  longitudinal  calculations 
as  evaluated  in  the  previous  sections  where  the  optical  and  piezo-acoustic  polarizations  are 
parallel  to  each  other.  When  evaluated  using  COMSOL  model  developed,  one  can  calculate  the 
resulting  strain  distribution  in  the  AO  material  and  the  optical  wave  propagation  in  the 
periodically  poled  structure,  as  illustrated  in  Figure  6.6. 


Periodic  strain  change, 
x-direction 


Figure  6.4  A  typical  example  of  an  acousto-optic  modulator  configuration 
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2400  nm 


400  nm 


Coupling  through 
refractive  index  ny 


Piezo 


o 

-Q 

C 

fD 


Input  port 


air 


LiNb03 


air 


Output  port 


Figure  6.5  The  setup  of  the  piezoelectric  and  electromagnetic  domains  for  the 
COMSOL  FEA  acousto-optic  model 
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Floquet  periodicity 


Notice  that  the  incident  signal  is  refracted  at  the  top  ainmatcrial  interface  and  then  again  at  the 
bottom  and  that  the  wavelength  is  smaller  inside  the  AO  material.  The  simulation  results 
presented  in  Figure  6.6  is  obtained  using  COMSOL  model  for  zero  order  diffraction  only 
considering  only  waves  transmitting  through  the  material.  The  z  component  (out  of  plane)  of  the 
electric  field  is  plotted.  The  angle  of  incidence  is  alpha=0.7  radians. 


Figure  6.6  COMSOL  simulation  of  the  optical  propagation  from  the  top 
section  (air)  into  the  material  and  to  the  bottom  port 
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6.6  Summary 


This  work  demonstrates  the  usefulness  of  examination  of  the  AO  orientation  dependence.  The 
elasto-optic  effects  and  figure  of  merit  can  be  obtained  for  any  crystal  class,  and  the  maximum 
directions  of  interaction  can  be  inspected.  Of  the  materials  investigated  only  two  showed 
maximum  interaction  along  the  principle  Cartesian  axes  demonstrating  the  importance  of  a 
proper  crystal  cut  for  ideal  interaction.  Germanium  gives  the  largest  FOM  among  the  crystals 
evaluated,  while  there  are  many  other  materials  that  give  larger  interaction  are  candidate  to  be 
evaluated  using  the  method.  The  elasto-optic  tensor  itself  does  not  always  provide  accurate  view 
of  the  best  direction  for  AO  interaction,  this  modeling  approach  allows  for  a  more  reliable 
method  of  determining  the  best  acousto-optical  coupling.  The  COMSOL  model  provides 
verification  of  these  results  and  the  possibilities  to  tailor  materials  for  improved  performance. 
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CHAPTER  7:  FUTURE  WORK 


7.1  Introduction 

This  work  describes  many  ways  of  using  and  studying  the  enhanced  electro-optic 
interaction  available  under  piezo-resonance;  yet  there  are  still  many  other  directions  to  which 
this  subject  can  be  extended. 

There  is  certainly  more  work  that  could  be  done  on  resonance  induced  transmission 
enhancement  as  described  in  section  3.7.  Using  a  high  field  signal  generator  will  allow  for  a 
more  complete  study  of  the  possible  improved  optical  transparency  that  may  results  from  proper 
modulation  at  certain  resonant  modes.  Moving  in  to  this  area  may  show  the  potential  of  using 
this  effect  for  optical  shutters. 

Due  to  the  significant  increase  in  coupling  the  sensitivity  of  these  devices  may  be  high 
enough  to  pick  up  radio  frequency  signals  without  active  amplification.  Such  a  configuration 
where  an  antenna  and  perhaps  a  few  passive  components  are  connected  to  the  sample  material 
may  allow  coupling  of  the  radio  wave  to  modulate  the  sample  directly.  This  configuration  would 
have  a  direct  application  for  radio  over  fiber  technology.  Radio  over  fiber  uses  long  optical  fibers 
to  connect  local  radio  communication  and  currently  require  individually  powered  OEO 
conversion  points.  The  development  of  passive  nodes  that  directly  impress  the  radio  signal  on  to 
the  fiber  could  be  very  advantageous. 

Most  work  on  resonance  enhancement  thus  far  has  been  concentrated  below  about  20 

MHz  and  has  been  limited  by  the  test  equipment  used;  it  would  be  interesting  to  investigate  how 

high  enhancement  can  go.  To  do  this  a  different  testing  configuration  will  be  needed.  One 

possibility  is  to  use  a  traveling  wave  electrode  configuration  such  as  those  used  in  commercial 

RF  fiber  optic  modulators.  This  configuration  may  work  well  because  wide  bandwidth 
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modulation  is  possible;  however  adapting  this  configuration  to  resonate  the  modulators  has  been 
found  to  be  difficult.  Another  possibility  is  using  a  microwave  cavity  such  as  used  in  [43-45], 
This  method  is  discussed  in  detail  in  section  7.1. 

There  is  also  some  significant  potential  for  using  this  technology  for  developing  a 
phonon-plasmon  (defined  in  Section  7.3.1)  metamaterial  device.  The  resonant  enhanced  process 
can  produce  a  periodic  charge  distribution  analogous  to  the  surface  plasmon  widely  studied  on 
metallic  surfaces.  This  method  is  discussed  in  detail  in  section  7.2. 

7.1  Resonance  Enhanced  Electro-Optic  Coefficient  by  means  of  a  Microwave  Cavity 

Here  a  rectangular  cavity  under  measurement  is  set  at  transverse  electric  mode  103 
(TE103).  The  sample  is  typically  placed  in  the  location  of  most  positive  slope  as  shown  in  Figure 
7.1.  When  the  sample  is  placed  at  a  point  of  monotonic  slope  change  there  is  an  induced 
refractive  index  change  in  the  sample  (Huang  et  al.  [44]).  Additionally  a  laser  was  used  and 
coupled  through  a  small  hole  at  both  ends  of  the  waveguide,  thus  a  microwave  field  is  induced 
on  the  sample  and  the  optical  output  can  be  measured.  This  refractive  index  change  will  cause 
the  optical  laser  pulse  to  shift  in  frequency  which  is  the  same  as  a  narrowing  or  widening  the 
pulse  in  the  time  domain  which  is  determined  by  measuring  its  full  width  at  half  max  (FWHM) 
prior  and  after  traversing  the  waveguide.  By  adjusting  the  power  of  the  network  analyzer  and 
monitoring  the  photodetector  voltage  a  linear  relationship  change  be  observed  between  optical 
transmission  intensity  squared  (  I  )  and  the  shifted  pulse;  the  slope  of  which  is  proportional  to 
the  EO  coefficient.  This  experiment  could  be  modified  by  applying  electrodes  to  the  sample  and 
modulating  at  the  resonant  frequency  of  the  sample.  This  technique  may  allow  for  an  observed 
increase  in  the  coupling  measurable  in  the  microwave  range  which  could  have  significant  impact 
on  microwave  electro-optic  communications  systems.  There  are  two  drawbacks  to  this  approach. 
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The  first  is  that  only  a  few  frequencies  can  be  tested,  i.e.  only  the  resonant  frequencies  of  the 
waveguide.  The  second  issue  is  that  the  use  of  a  sample  with  parallel  plate  electrodes  inside  the 
waveguide  may  disturb  field  distributions  that  are  difficult  to  predictor  demand  high  field 
strength.  There  has  been  some  work  to  model  the  waveguide  with  inserted  sample  using  finite 
element  analysis  done  by  the  author  and  it  is  included  in  the  appendix  of  this  dissertation.  The 
appendix  focuses  on  a  method  of  determining  the  microwave  anisotropic  permittivities  of 
samples.  Determining  microwave  permittivity  is  key  to  an  accurate  modeling  of  the  electroded 
electro-optic  sample  in  the  waveguide  and  to  the  characterization  of  the  high  frequency  response 
of  the  sample. 
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Figure  7.1  Microwave  cavity  electro-optic  measurement  technique  [44] 
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7.3  Acousto-optic  Metamaterials 


7.3.1  Plasmonics 

Plasmons  are  electron  density  waves  and  can  be  in  the  form  of  surface  or  bulk  plasmons. 
Surface  plasmons  can  exist  on  the  interface  between  a  metal  and  dielectric;  this  is  a  sufficient 
condition  but  not  necessary.  The  wave  must  be  on  the  interface  of  two  materials  with  complex 
permittivity  of  opposite  sign,  the  sum  of  their  permittivity  must  be  negative  [46].  The  plasmon 
frequency  oop  is  the  point  where  the  electron  gas  of  a  plasmon  has  strong  vibrations.  Higher 
electron  density  results  in  this  frequency  also  being  high.  The  plasmon  energy  (in  electron  volts) 
is  calculated  from  this  frequency  Ep  —  ha>p  and  at  this  energy  the  external  field  vibration  is  so 
fast  that  the  electrons  in  the  metal  cannot  follow  and  the  metal  stops  being  reflective. 

7.3.2  Metamaterials  and  Cloaking 

Metamaterials  in  the  most  general  sense  are  manmade  or  engineered  materials  that  are 
non-homogenous.  Most  commonly  this  term  applies  to  negative  index  materials;  however  that  is 
just  one  type  of  metamaterial.  As  described  by  Smith,  D.  R.,  J.  B.  Pendry,  et  al.  [47] 
conventional  materials  have  properties  that  derive  from  their  constituent  atoms  while  the 
properties  of  metamaterials  are  from  their  constituent  units.  Metamaterials  are  used  to  influence 
the  wave  interaction  with  a  material,  as  long  as  the  size  and  spacing  is  small  compared  to  the 
wavelength  such  that  the  wave  cannot  tell  the  difference  between  an  engineered  group  of  objects 
or  a  homogeneous  medium.  Another  class  of  these  materials  is  known  as  transformative  optics. 
These  materials  have  an  anisotropic  refractive  index  that  changes  spatially  in  the  material. 

Metamaterial  devices  can  be  designed  to  create  a  “cloaking”  effect,  or  obstructing  the 
wave  propagation  altering  the  perception  of  the  “normal”  detection.  The  Electromagnetic  wave  is 
caused  to  bend  around  the  object  and  render  it  invisible  to  a  far  away  observer.  The  resonant 
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enhanced  devices  discussed  in  this  dissertation  are  known  to  increase  the  transmission  of  the 


device  and  could  potentially  be  used  in  a  similar  way  to  the  aforementioned  metamaterials  to 
provide  dynamically  adjustable  cloak.  Typical  cloaking  device  configurations  [48,  49]  disguise  a 
dielectric  cylinder  with  a  larger  cloaking  material  outside  of  the  cylinder.  When  a  plane  wave 
hits  the  cloak  the  wave  bends  around  and  reforms  on  the  other  side  as  if  the  cylinder  and  cloak 
were  never  there.  Were  the  cloak  not  in  place  the  wave  would  scatter  causing  an  observer  on  the 
opposite  side  to  see  a  shadow.  Alternatively  the  resonance  enhanced  technology  may  be  able  to 
achieve  this  type  of  interaction  with  the  added  benefit  of  dynamic  control.  As  seen  in  the  figure 
below  the  current  configuration  uses  a  vibrating  crystal  plus  external  optics  to  control 
transmission.  Because  the  transmission  is  induced  dynamically  by  the  piezoelectric  resonance  it 
can  be  easily  turned  on  and  off  causing  the  object  to  go  from  visible  to  invisible.  The  figure 
demonstrates  the  concept  of  the  scattering  of  the  incident  wave.  Turning  the  cloak  on  results  in  a 
minimum  scattering  cross  section  (SCS),  when  scattering  is  small  an  observer  on  the  opposite 
side  of  the  object  (relative  to  the  applied  wave)  would  not  see  the  shadow  of  the  object. 
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Figure  7.2  Proposed  configuration  for  achieving  piezo-resonance  induced  optical 
cloaking 


7.3.3  Proposed  Phonon-Plasmon  Meta-Device 

To  study  this  effect  and  move  the  research  in  the  direction  of  plasmonics  and 

metamaterials  a  structure  is  proposed.  This  structure  will  make  use  of  surface  acoustic  waves  as 

seen  in  Figure  7.3.  In  this  configuration  coplanar  interdigized  electrodes  are  used  to  induce  a 

surface  wave  by  way  of  the  piezoelectric  effect.  A  high  frequency  signal  is  applied  to  one  set  of 

electrodes  and  the  other  is  grounded.  The  induced  elastic  wave  will  travel  across  the  surface  from 

left  to  right  in  the  figure  and  a  periodic  strain  will  be  induced.  The  periodic  strain  results  in  local 

stress  through  the  fourth  rank  elastic  stiffness  tensor 
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(7.1) 


The  stress  is  converted  to  a  polarization  by  the  direct  piezoelectric  effect  (Equation  (3.3)). 


Interdigitated  electrodes 


Surface  acoustic  wave  and 
periodic  charge  distribution 


Figure  7.3  Proposed  surface  acoustic  wave  device  to  launch  wave  and  induce 
a  periodic  polarization 

The  induced  periodic  charge  distribution  may  be  a  vehicle  to  allow  propagation  and 
enhancement  of  the  plasmon  wave.  This  will  require  a  surface  charge  of  sufficient  amplitude  and 
proper  periodicity.  A  proposed  test  scheme  for  experimentally  verifying  this  interaction  [50]  is 
shown  in  Figure  7.4.  A  signal  generator  is  used  to  control  the  SAW  device  and  induce  the 
periodic  charges.  A  tunable  laser  is  used  to  launch  the  plasmon  on  the  device  surface  and  a 
spectrometer  picks  up  the  signal.  A  vibrometer  is  also  integrated  in  this  experiment  which  can 
be  used  to  verify  the  existence  of  the  surface  wave. 


125 


Ultra-High  Frequency 
Vibrom  eter 


SAW  device 


Figure  7.4  Proposed  experiment  setup  for  Phonon-Plasmon  measurement 

A  preliminary  simulation  of  this  device  has  been  developed  to  test  the  scenario  for  inducing  a 
reliable  periodic  polarization.  This  model  is  based  on  a  similar  model  in  the  COMSOL 
Multiphysics  example  library  [51].  The  simulation  is  done  in  two-dimensions  on  the  xy  plane  (z 
dimension  is  out  of  plane)  as  seen  in  Figure  7.5.  The  sample  used  is  PMN-30%PT  using  the 
same  material  properties  as  given  in  Table  4.3.  Three  sets  of  interdigitized  electrodes  are  on  the 
sample  surface,  the  other  area  to  the  right  of  this  is  an  88  pm  long  section  where  the  wave 
excitation  will  occur.  It  is  assumed  that  this  device  is  mounted  on  another  material  thus  the 
bottom  boundary  is  a  fixed  constraint. 
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Figure  7.5  Finite  Element  Analysis  setup  for  2D  SAW  model 

The  device  was  excited  with  a  lVolf  peak,  786  MHz  signal  and  the  results  are  shown  in  Figure 
7.6.  The  physical  displacement  (in  black)  on  the  surface  is  sinusoidal  as  expected  with  maximum 
of  about  20  micrometers.  The  polarization  in  the  y  dimension  (in  blue)  is  in  phase  with  the 
displacement  with  a  maximum  of  about  2  nC/m  (or  0.2  pC/cnT)  which  is  a  relatively  low  level 
of  polarization.  The  wavelength  of  the  standing  wave  is  approximately  5  pm. 
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SAW  device  at  786  MHz 


y  dimension  displacement 
y  dimension  polarization 


Figure  7.6  Simulation  of  the  y  dimension  surface  displacement  and  resulting 
polarization  in  a  SAW  device  excited  at  786  MHz 


7.4  Summary 

The  results  in  the  research  involved  in  this  dissertation  have  opened  up  many  new 
possibilities  for  further  study,  of  which  a  few  have  been  discussed  in  this  chapter.  Enhancing  the 
electro-optic  interaction  at  microwave  frequencies  as  discussed  in  section  7.2  is  a  very  important 
area  of  research.  Fiber  optic  communications  is  a  vital  area  for  transferring  data  which  depends 
on  electric  modulation  to  transfer  the  information.  It  is  this  link  between  electronics  and  optics 
that  is  often  the  bottle  neck  in  a  system,  which  shows  the  need  for  electro-optic  modulators  that 
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y  dimension  polarization  (C/m  ) 


can  modulate  at  high  frequencies.  The  electro-optic  effect  decreases  significantly  in  the 
microwave  range  which  results  in  lower  EO  coefficients  and  thus  requiring  larger  devices  that 
need  more  power  to  operate.  Energy  is  one  of  the  biggest  issues  being  studied  today  and  this 
includes  energy  efficiency.  In  this  way  the  enhanced  electro-optic  effect  may  be  able  to  make  a 
significant  impact  on  a  major  issue  of  our  future  digital  communications. 

Section  7.3  also  discusses  a  very  import  issue  as  cloaking  is  a  technology  directly 
relatable  to  national  defense.  The  study  of  plasmon-phonon  coupling  will  likely  play  a  critical 
role  in  greater  understanding  of  resonance  enhancement.  Using  this  gained  knowledge  one  is  a 
step  closer  to  acousto-optic  metamaterials  that  offer  greater  control  of  the  optical  transmission 
necessary  for  optic  (or  acoustic)  cloaking. 
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CHAPTER  8:  CONCLUSION 


This  work  has  explored  and  integrated  many  areas  of  research  toward  the  central  goal  of 
enhanced  sensitivity  for  microwave  optoelectronic  devices.  This  enhanced  sensitivity  through 
piezoresonance  coupling  has  implications  for  electro-optic  modulators  and  other  devices 
including  metamaterial  cloaking.  The  applications  for  this  technology  can  easily  range  from 
commercial  telecommunications  and  consumer  electronics,  to  military  defense.  Electro-optic, 
mechanical  vibration  and  admittance  experiments  have  verified  the  high  degree  of  coupling 
available  as  a  result  of  synchronized  crystal  vibration.  The  addition  of  a  high  frequency 
vibrometer  measurement  has  shown  the  ability  of  induced  local  field  due  to  piezoelectric 
polarization  in  the  sample  to  influence  the  optical  signal  propagation  even  after  the  crystal  begins 
to  become  mechanically  clamped.  FEA  models  have  been  constructed  and  implemented  for 
modeling  the  crystal  vibration,  admittance  spectrum,  linear  electro-optic  effect,  photoelastic 
effect,  and  optical  wave  propagation.  These  models  are  consistent  with  observed  experimental 
results  and  have  provided  new  insights  into  the  dynamic  nature  of  the  induced  periodic 
displacive  current  in  a  resonating  sample;  this  gives  evidence  for  the  proposed  current 
distribution  model  proposed  in  [4],  Advanced  time  domain  simulations  have  also  given 
verification  of  the  possibility  of  broad  bandwidth  enhancement  resulting  from  AC  resonance 
biasing,  which  was  shown  experimentally  in  [6],  Models  have  verified  the  influence  of  periodic 
ferroelectric  domain  structure  on  enhancement  in  microwave  devices.  Models  have  also  been 
used  to  verify  the  potential  of  using  piezoresonant  enhancement  in  optical  modulators  and  their 
advantages  over  normal  and  electrically  resonant  phase  modulators  of  identical  size, 
composition,  and  configuration. 
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The  ultra-high  frequency  vibrometer  measurements  have  also  yielded  interesting  results. 
The  longitudinal  displacement  values  have  helped  to  reveal  the  complex  nature  of  the  samples 
vibration.  The  high  frequency  displacement  has  been  shown  to  decrease  as  the  crystal  becomes 
more  constrained  while  the  electro-optic  signal  remained  relatively  high.  The  acceleration  of  the 
sample  surface  gives  a  similar  response  showing  strong  peaks  even  when  the  displacement  has 
decreased.  This  gives  an  interesting  perspective  suggesting  that  the  surface  acceleration  may  play 
a  significant  role  in  enhancement  and  may  perhaps  lead  to  an  interesting  investigation  in  the 
future. 

In  addition  to  resonant  effects  the  acousto-optic  interaction  was  studied  to  facilitate  a 
novel  solution  to  more  efficient  sensors.  Optimizing  AO  devices  is  more  than  just  the  photo 
elastic  relationship  and  this  work  has  displayed  the  formation  of  pertinent  material  parameters 
and  equations  for  identifying  the  maximum  interaction  directions  for  materials  of  any  point 
group  symmetry.  The  maximum  interaction  directions  were  shown  to  not  always  be  along  the 
Cartesian  principle  axes.  These  maximum  directions  typically  occur  where  the  acoustic  velocity 
is  minimized  and  the  refractive  index  is  maximized  (yielding  minimum  light  velocity),  the  longer 
the  two  waves  take  to  co-propagate  in  the  material  the  more  they  interact  and  couple.  The  Finite 
Element  model  of  acousto-optic  effect  and  maximum  interaction  directions  found  by  the  Figure 
of  Merit  model  give  consistent  conclusions. 

Overall  this  dissertation  demonstrates  a  wide  array  of  methods  for  enhancement/increased 
sensitivity  in  electro-optic  devices.  Piezoelectric  resonance  enhancement  has  implications  far 
beyond  the  implementation  of  electro-optic  modulators  operated  at  low  frequency  vibration 
modes.  Synchronization  of  acoustic  phonons  with  optical  photons  can  open  significant  potential 
for  low  voltage,  high  sensitivity  devices  with  wide  bandwidth  operation. 
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APPENDIX 


Appendix  McIntosh,  R.,  A.  Bhalla,  et  al.  (2010).  Dielectric  Anisotropy  of  Ferroelectric  Single 
Crystals  in  Microwave  C-Band  by  Cavity  Vectorial  Perturbation  Method. 

Advances  in  Electroceramic  Materials  II,  John  Wiley  &  Sons,  Inc.:  75-88. 


A.l  Abstract 

Integrating  the  numerical  simulation  of  electromagnetic  field  in  a  perturbed  cavity  with 
the  microwave  measurement  via  swept  frequency  technique  using  a  network  analyzer,  a 
numerically  enhanced  cavity  vectorial  perturbation  (NECVP)  method  is  developed  and  reported 
in  this  paper.  The  NECVP  method  is  capable  of  resolving  anisotropic  dielectric  properties  of 
various  dielectric  and  ferroelectric  single  crystals  (sr  value  in  a  wide  range  from  10  to  10  has 
been  measured).  The  configuration  used  in  the  present  study  is  designed  to  measure  properties  in 
microwave  3-6  GHz  (part  of  the  IEEE  C-band).  Continued  signal  monitoring  and  automatic 
calculation  of  complex  permittivity  are  made  via  a  Lab  VIEW  interface  which  aids  the  data 
collection  process.  Numerical  simulation  of  the  cavity  field  is  carried  out  using  a  finite  element 
analysis  software  package  (COMSOL).  The  reliability  of  the  deduced  dielectric  permittivity 
(and  to  a  lesser  extend  the  dielectric  loss)  by  the  NECVP  method  is  found  to  be  quite  high 
limited  only  by  the  resolution  of  the  numerical  simulation  conducted,  which  becomes  more 
demanding  when  perturbation  of  the  microwave  field  is  significant.  Additionally  the  directional 
dependence  of  dielectric  permittivity  of  a  given  sample  with  arbitrary  shape  can  be  obtained 
using  a  sample  rotation  technique. 

A.2  Introduction 

Cavity  perturbation  method  has  been  widely  used  to  study  the  dielectric  and  magnetic 
parameters  in  the  microwave  frequency  region.  The  dispersive  and  dissipative  terms  of  the 
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materials  are  directly  related  to  the  change  in  the  resonant  frequency  and  the  quality  factor  of  the 
cavity  from  the  respective  empty  cavity  values  [52],  Conventional  microwave  cavity  perturbation 
techniques  have  been  known  as  fast  and  convenient  methods  for  evaluating  gigahertz  dielectric 
permittivities  of  materials  that  are  typically  isotropic,  of  low  sr,  and  small  (compared  to  the 
wavelength)  in  size.  However,  a  cavity  perturbation  technique  that  accurately  evaluates 
anisotropic  dielectric  properties  especially  of  those  highly  polarizable  materials  (e.g., 
ferroelectric  materials  sr  »20)  at  GHz  frequencies  has  not  been  available. 

The  resonant  frequency  fr  and  the  quality  factor  Q  of  a  rectangular  cavity  waveguide,  for  a  given 
standing  wave  TEion  mode  in  the  microwave  region,  are  expressed  by  the  following  equations: 
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where  a,  d,  so,  po,  N,  Rs,  and  r\  are  respectively  the  width  of  the  waveguide,  length  of  the 
waveguide,  permittivity  of  free  space,  permeability  of  free  space,  mode  number,  surface 
resistance  of  the  cavity,  and  intrinsic  impedance.  Inserting  a  sample  into  the  cavity  causes  a  shift 
in  the  resonant  frequency  f.  and  a  change  in  the  quality  factor  Q  of  the  waveguide.  This 
perturbation  is  dependent  on  the  relative  volumes  of  the  cavity  Vc  and  of  the  sample  within  the 
cavity  Vs,  the  permittivity  sr  of  the  sample  and  thus  the  electric  field  concentration  in  the  cavity 
Ec  and  in  the  sample  Es  as  shown  in  the  general  perturbation  equation  below:  [53] 
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A  determination  of  the  complex  permittivity  by  the  perturbation  technique  thus  is  dependent  on 
both  these  changes  and  the  integration  of  the  electric  field  over  the  volumes  of  the  sample  within 
the  cavity.  For  a  sample  with  length  parallel  to  the  electric  field  direction,  assuming  a  small 
perturbation  of  the  field,  the  above  relation  can  be  easily  simplified  and  used  to  find  the  complex 
permittivity  [54],  Under  these  conditions  the  electric  fields  in  the  sample  (£’5)  and  in  the  cavity 
(. Ec )  are  approximately  equal;  one  can  then  derive  the  following  expressions  for  the  real  and  the 
imaginary  parts  of  the  complex  permittivity: 
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Compared  with  the  Hakki-Coleman  post  resonant  method  which  requires  typically 
centimeter  diameter  samples,  [55]  the  cavity  perturbation  technique  offers  flexibility  that  permits 
measurement  of  mm  size  samples  of  essentially  any  shape. 

Considerable  errors  however,  using  the  conventional  microwave  perturbation  method  by 
Equations  (A.4)  and  (A. 5),  may  result  in  the  values  of  complex  dielectric  pennittivity  when  the 
assumption  of  small  perturbation  is  invalid.  Santra  and  Limaye  [56]  described  an  approach  using 
finite  element  method  to  assist  in  evaluating  complex  permittivity  of  arbitrary  shape  and  size  and 
demonstrated  the  application  of  such  method  to  measure  isotropic  dielectric  materials  with 
moderate  permittivity  (s’~22).  In  the  current  study  the  integration  of  numerical  finite  element 
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analysis  with  microwave  measurement  is  extended  to  account  for  directional  dependence  of 

3 

samples  and  to  ferroelectric  materials  with  sr’  up  to  10  . 

A.3  Experimental 

A  rectangular  waveguide  was  used  for  all  experiments  reported  in  this  paper.  The 
dimensions  (in  centimeters)  are  b=4.74  in  width,  a=2.21  in  height,  and  d=18.2  in  length.  All 
samples  were  measured  at  modes  TE103  and  TE105  with  cavity  resonant  frequencies  of  4.01GHz 
and  5.19  GHz  at  room  temperature,  respectively.  Using  this  method  samples  are  placed  in  the 
geometric  center  of  the  waveguide  where  the  electric  field  strength  is  maximum  for  the  odd 
modes.  The  waveguide  also  has  custom  holes  in  four  of  its  faces  to  allow  insertion  of  the  sample 
and  a  hollow  metal  tube  welded  to  each  of  the  holes  to  prevent  energy  loss  (See  Figure  A.l). 


Figure  A.l  Configuration  and  dimension  of  the  microwave  cavity  used  for  the 
NECVP  study 

For  a  typical  rectangular  sample  there  are  four  orientations  of  interest  illustrated  in  Figure 
A. 2(a)  through  Figure  A.2  (d),  hereafter  referred  to  as  positions  PE,  PF,  HE,  and  HF,  where  P 
and  H  designate  whether  length  direction  of  the  sample  is  parallel  (P)  or  horizontal  (H)  to  the 


electric  field,  E  and  F  designate  whether  the  sample’s  edge  (E)  or  face  (F)  is  inline  with  the 
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direction  of  propagating  EM  wave.  The  designations  are  as  follows:  (a)  Position  PE:  Length 
parallel  to  E  field,  Edge  inline  with  wave  propagation,  (b)  Position  PF:  Length  parallel  to  E  field, 
Face  inline  to  wave  propagation,  (c)  Position  HE:  length  horizontal  to  E-field,  Edge  inline  with 
wave  propagation,  (d)  Position  HF:  Length  horizontal  to  E-field,  Face  in  line  with  wave 
propagation.  Additionally  using  this  method  a  sample  oriented  with  its  longest  axis  perpendicular 
to  the  electric  field  (as  in  Positions  HE  &  HF)  can  be  rotated  to  investigate  the  permittivity  at 
different  angles. 


Figure  A. 2  Orientation  of  samples  relative  to  microwave  resonant  cavity 

Both  low  permittivity  and  high  permittivity  samples  were  studied  in  the  work,  some  were 
chosen  for  their  well  known  dielectric  properties  at  microwave  frequency,  some  were  studied  to 
demonstrate  the  capability  of  the  method.  The  types  of  the  samples  and  their  dimensions  are 
summarized  in  Table  A.  1 . 
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Table  A.l  Summary  of  sample  dimensions  used  in  microwave  measurement 


Sample 

Shape 

Dimensions  (mm) 

Alumina  Ceramic  Substrate  (99.6%  AFO3) 

Bar 

L=0.63,  W=2.64,  H=16.25 

Coming  0080  glass 

Bar 

L=1.19,  W=1.21,  H=10.26 

Fused  Silica 

Tube 

Dout=7.97,  Dm=5.96,  L=206 

Pyrex  Glass 

Rod 

D=3.92,  L=125 

Teflon  (PTFE) 

Rod 

D=4.8,  L=50 

PMN-PT  single  crystal,  Bridgman  method 

Bar 

L=2.7,W=2.96,  H=4.69//[l  10] 

PZT-5  ceramic,  commercial 

Rod 

D=3.36,  L=3.74 

A  multiphysics  software  package  (COMSOL)  was  used  for  the  finite  element  simulation 
of  the  waveguide  and  the  samples.  All  simulations  were  done  for  the  TE103  and  TE105  modes.  The 
conductivity  of  the  walls  of  waveguide  was  given  in  accordance  with  that  of  copper  and  the  four 
round  ports  were  included  in  the  modeling.  In  this  report  all  the  samples  were  modeled  using 
their  actual  physical  dimensions  but  considered  as  having  isotropic  permittivity  (anisotropic 
modeling  will  be  the  content  of  an  upcoming  report).  A  plot  of  the  perturbed  cavity  electric  field 
(V/m)  with  inserted  sample  of  assumed  dielectric  permittivity  sr=10  is  shown  in  Figure  A.3.  The 
sample  has  dimensions  2.7x2.96x4.69  mm  and  is  positioned  in  the  center  of  the  waveguide  in 
PE  position. 
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Figure  A. 3  Finite  element  simulation  of  the  perturbed  cavity  electric  field 
(V/m),  for  TE103  mode  near  /c=4.003GHz 

By  placing  a  sample  in  the  center  of  the  waveguide  where  \Eq\  is  at  the  maximum  and  by 
assigning  permittivity  of  the  sample  as  a  variable,  a  parametric  sweep  was  performed  to  calculate 
the  fs  as  one  of  the  outputs  of  the  finite  element  analysis.  For  a  given  sample  and  sample 
dimension  the  normalized  resonance  frequency  shift  (fc-fsVfs  can  be  matched  to  the  experimental 
shift  and  thus  to  determine  the  real  permittivity.  Once  the  real  permittivity  is  known  then  the 
imaginary  part  is  resolved  by  an  additional  parametric  sweep  to  match  the  normalized  change  in 
Q,  in  form  of  (1/Qs)-(1/Qc)-  For  the  given  waveguide  the  computation  time  is  less  than  an  hour 
using  a  HP  workstation  (HP  Z80  22  GB  RAM),  depending  on  the  spatial  resolution  desired. 
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A.4  Results  and  Discussion 


A.4.1  Limitations  of  the  Conventional  Perturbation  Method 

Figure  A.4  shows  a  plot  of  normalized  frequency  shift  vs.  real  part  of  the  relative 
permittivity  for  a  sample  of  rectangular  shape  in  PE  position.  The  length  of  the  sample  is  such 
that  its  ends  touch  both  sides  of  the  walls  of  the  cavity.  The  plot  is  presented  to  compare  the 
results  of  sr’  by  FEA  simulation  and  by  the  calculation  from  the  conventional  perturbation 
equations  (Equations  (A.4)  and  (A. 5)).  The  two  methods  are  in  reasonable  agreement  for  low 
shifts;  however  they  start  to  diverge  when  the  normalized  resonance  frequency  shift  (fc-fs)/fs 
>0.11  where  their  discrepancy  becomes  >20%  and  the  conventional  evaluation  becomes  invalid 
as  sr’  approaches  100. 
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Simulation,  bar  sample  2. 7x2. 96x22. 05mm,  real  permittivity  Vs. 


Figure  A.4  Comparison  of  calculated  resonance  frequency  shift  as  function  of 
dielectric  permittivity,  by  FEA  and  by  Equations  (4)  and  (5) 

Another  comparison  using  a  shorter  sample  (L=4.69mm)  is  shown  in  Figure  A. 5.  The  percent 

difference  becomes  greater  than  20%  when  the  real  part  of  the  permittivity  is  above  4.  The 

comparison  demonstrates  the  limitations  of  the  simplified  conventional  perturbation  formula 

used  for  samples  shorter  than  the  height  of  the  cavity.  The  conventional  method  in  such  a  case 

underestimates  the  permittivity  of  a  sample  and  is  essentially  invalid  for  sr’>10.  The  reason  for 

these  large  errors  at  higher  permittivity  becomes  clear  when  we  observe  the  electric  field  profile 

upon  increasing  permittivities  of  the  sample  inside  the  cavity  ( 

Figure  A. 6).  For  permittivity  of  1  the  cavity  field  is  well  defined  and  is  sinusoidal;  however  as 

the  permittivity  increases  the  perturbation  to  the  electric  field  intensifies.  Furthermore  as  the 

permittivity  increases  the  electric  field  will  start  to  bend  around  the  sample  and  will  have  less 
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effect  on  the  frequency  shift  causing  a  decrease  in  measurement  sensitivity.  To  increase  the 
sensitivity  for  samples  of  high  permittivity,  longer  sample  dimension  (//  to  the  electric  field 
direction)  and/or  higher  signal  level  should  be  considered. 


Simulation,  bar  sample  2. 7x2. 96x4. 69mm,  real  permittivity  Vs. 
normalized  frequency  shift,  Position  PE,  TE103 


Figure  A. 5  Comparison  of  calculated  resonance  frequency  shift  as  a  function 
of  dielectric  permittivity,  by  FEA  and  by  Equ.  (9.4)  and  (9.5). 

The  bar-shaped  sample  (2.7x2.96x4.69  mm3)  is  PE  positioned  and 
with  dimension  shorter  than  the  height  of  the  cavity 
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Simulation,  bar  Sample  2.7rnmx2.96rrrnx4.69rrrn,  Electric  field  squared  for 
various  permittivities,  Position  PE,  TE,03 
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Figure  A. 6  Perturbed  electric  field  intensity  profile  in  the  microwave  cavity  of 
TE103  mode  as  a  function  of  the  real  part  of  the  permittivity  for  a 

"5 

PE  center  positioned  sample  (2.7x2.96x4.69  mm  ) 


In  general,  there  are  several  methods  [53]  for  calculating  the  electric  field  concentration 
in  the  sample;  These  methods  are  useful  but  only  for  certain  sample  types  and  orientations.  The 
electric  field  concentration  in  a  cylindrical  sample  where  the  ends  terminate  at  the  waveguide 
walls  is  fairly  easy  to  determine  however  the  problem  becomes  far  more  difficult  for  a  sample 
that  does  not  terminate  at  the  waveguide  walls.  Furthermore  it  becomes  difficult  to  measure  high 
permittivity  dielectrics  in  the  microwave  region  due  to  the  skin  effect.  The  electric  field  becomes 
less  concentrated  in  the  sample  and  thus  sensitivity  decreases. 


A.4.2  Verification  of  NECVP  Using  Low  Permittivity  Materials 

Table  A. 2  summarizes  the  complex  permittivity  (expressed  in  terms  of  sr’,  sr”,  and  tan8) 

obtained  using  the  procedure  described  earlier  (the  NECVP  method)  for  alumina  ceramic  sample 

of  thin  rectangular  shape,  measured  in  all  four  standard  positions.  The  table  shows  that  the  real 
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permittivity  of  the  sample  is  about  8.4  at  ~4GHz  which  is  in  close  agreement  of  reported  values 
(about  9  at  1MHz)  [57].  The  signal  levels  for  positions  HE  and  HF  were  too  low  to  measure 
which  is  sometimes  a  problem  for  low  loss  materials  with  length  dimension  positioned 
perpendicular  to  the  electric  field.  The  imaginary  permittivity  is  about  5x10"  and  the  land  loss 
factor  is  6x10  4  for  positions  PE  and  PF  at  ~4GHz.  The  reported  [57]  values  are  sr”=0. 003 -0.02 
and  tan  8=0.0003-0.002  (at  1MHz). 


Table  A. 2  Alumina  complex  permittivity  measured  by  NECVP  method  near  4.01GHz  (TE103) 
and  5.19  (TEi05)  GHz 


Alumina  Ceramic  Substrate  (99.6%  AEO3)  -  Bar 


Position 

Mode 

£r’ 

+/- 

Sr” 

+/- 

tan  8 

+/- 

PE 

TE103 

8.452 

0.009 

0.005 

0.007 

0.0006 

0.0008 

TE105 

8.382 

0.008 

0.015 

0.016 

0.0018 

0.0019 

PF 

TE103 

8.452 

0.009 

0.005 

0.007 

0.0006 

0.0008 

TE105 

8.359 

0.015 

0.007 

0.0086 

0.0008 

0.0010 

HE 

TE103 

8.376 

0.046 

- 

- 

- 

- 

TE105 

7.900 

0.281 

- 

- 

- 

- 

HF 

TE103 

7.736 

0.047 

- 

- 

- 

- 

TE105 

7.521 

0.105 

- 

- 

- 

- 

Table  A. 3  summarizes  the  complex  permittivity  obtained  using  NECVP  method  for  a  rectangular 
sample  of  Corning  0080  glass.  The  same  material  of  larger  disk  shape  was  also  tested  with  the 
post  resonant  technique  (See  Table  A. 4)  and  the  results  using  both  methods  agree  reasonably 
well. 

Table  A. 5  shows  the  results  obtained  using  NECVP  method  for  a  fused  silica  in  tubular 
form.  The  complex  permittivity  for  bulk  fused  silica  at  5  GHz  was  reported  [58]  to  be  ~3.8  and 
the  tanS  of  the  value  between  0.2  to  0.4x10  4.  The  results  obtained  show  somewhat  higher  real 
part  of  the  permittivity  and  also  higher  tanS;  however  in  the  same  time  both  PE  and  HE  positions 
yield  consistent  results.  Thus  the  results  are  likely  to  be  reliable  taking  into  consideration  that 
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surface  treatment  is  often  applied  onto  fused-silica  tubing  to  decrease  surface  activity  or  inhibit 
UV  absorption,  which  may  account  to  an  increase  in  the  dielectric  loss. 


Table  A. 3  Corning  0080  glass  complex  pennittivity  measured  by  NECVP  method  near  4.01GHz 
(TE103)  and  5.19  (TE105)  GHz 


Coming  0080  -  Bar 

Position 

Mode 

£r’ 

+/- 

£r” 

+/- 

tan  8 

+/- 

PE 

TEios 

6.151 

0.012 

0.088 

0.010 

0.0143 

0.0016 

TEios 

6.103 

0.012 

0.093 

0.039 

0.0153 

0.0064 

PF 

TEiq3 

6.170 

0.014 

0.063 

- 

0.0101 

- 

TEios 

6.100 

0.030 

0.051 

0.0425 

0.0084 

0.0121 

HE 

TEios 

5.923 

0.021 

0.125 

- 

0.0211 

- 

TEios 

5.780 

0.054 

0.098 

0.152 

0.0169 

0.0216 

HF 

TEios 

5.390 

0.020 

0.125 

- 

0.0232 

- 

TEios 

5.288 

0.054 

0.125 

- 

0.0236 

- 

Table  A.4  Corning  0080  glass  complex  permittivity  measured  by  post  resonant  technique 


Sr’ 

£r” 

tan  8 

3  GHz 

6.715 

0.1477 

0.0220 

6  GHz 

6.750 

0.1256 

0.0186 

Table  A. 5  Fused  silica  tube  complex  pennittivity  measured  by  NECVP  method  near  4.01GHz 
(TE103)  and  5.19  (TE105)  GHz 


Fused  Silica  -  Tube 

Position 

Mode 

Sr’ 

+/- 

sr” 

+/- 

tan  8 

+/- 

PE 

TEios 

4.248 

0.001 

0.0332 

0.0003 

0.0078 

0.0001 

TEios 

4.213 

0.001 

0.0319 

0.0010 

0.0076 

0.0002 

HE 

TEios 

4.244 

0.001 

0.0327 

0.0005 

0.0077 

0.0001 

TEios 

4.197 

0.004 

0.0314 

0.0024 

0.0075 

0.0006 

Table  A. 6  and  Table  A. 7  summarize  the  results  for  a  Pyrex  glass  rod  and  a  Teflon  rod. 
For  the  Corning  Pyrex  7740  glass,  the  measurement  results  are  consistent  with  the  reported 
dielectric  constant  (sr=4.6)  and  loss  tangent  (tan8=  0.004)  at  20  °C  and  1  MHz,  respectively  [59]. 
For  PTFE  Teflon  rod  the  real  part  of  dielectric  pennittivity  measured  at  room  temperature 
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compares  well  with  those  previously  reported  (sr’=2.055  at  300K  and  9.93  GHz)  [60],  Although 
the  NECVP  technique  is  highly  sensitive,  no  further  interpretation  on  the  physical  meanings  of 
the  differences  on  measured  results  between  the  PE  and  the  HE  positions  for  the  presumably 
isotropic  Pyrex  and  the  Teflon  rods  can  be  given  without  rigorous  calibration  of  the  3D  field 
distribution  of  the  cavity,  other  than  noting  here  that  certain  uncertainties  exist  caused  by  reasons 
in  addition  to  dimensional  variations. 

Table  A. 6  Complex  permittivity  of  Pyrex  glass  rod  measured  by  NECVP  method  near  4.01GHz 
(TE103)  and  5.19  (TE105)  GHz 


Glass  -  Rod 

Position 

Mode 

Sr’ 

+/- 

sr” 

+/- 

tan  8 

+/- 

PE 

TE103 

3.9561 

0.0006 

0.0294 

0.0010 

0.0074 

0.0002 

TE105 

3.9248 

0.0012 

0.0271 

0.0019 

0.0069 

0.0004 

HE 

TE103 

4.6386 

0.0022 

0.0291 

0.0028 

0.0062 

0.0006 

TE105 

4.5704 

0.0108 

0.0291 

0.0028 

0.0063 

0.0006 

Table  A. 7  Permittivity  of  Teflon  rod  measured  by  NECVP  method  near  4.01GHz  (TE103)  and 
5.19  (TE105)  GHz 


Teflon  (PTFE)  -  Roc 

1 

Position 

Mode 

Sr’ 

+/- 

HE 

TE103 

2.0557 

0.0002 

TE105 

2.0465 

0.0019 

PE 

TE103 

1.8876 

0.0001 

TE105 

1.8788 

0.0006 

A.4.3  Anisotropic  Dielectric  Property  Evaluation 

One  of  the  challenges  of  the  NECVP  method  is  orienting  the  sample  in  the  cavity  in  a 
reproducible  manner.  A  jig  was  constructed  to  allow  simple  manipulation  of  the  sample.  Two 
pieces  of  string  were  strung  through  the  middle  axis  of  the  waveguide  (perpendicular  to  E  field) 
and  secured  on  either  side  (See  Figure  A. 7  and  Figure  A. 8).  The  sample  under  test  is  pinched 


between  the  two  strings.  When  the  strings  are  given  tension  the  sample  is  secure  and  by 
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simultaneously  rotating  the  top  and  bottom  ends  of  the  string  the  sample  is  rotated  inside  the 
waveguide.  This  setup  ensures  that  the  sample  stays  centered  and  the  string  has  only  a  very  small 
effect  on  the  perturbation.  The  empty  cavity  gives  a  resonant  frequency  of  3.9877519  GHz  +/- 
880Hz  while  the  two  pieces  of  thread  give  3.9877002  GHz  +/-510  Hz.  This  gives  a  normalized 
shift  of  only  0.0000129  +/-0. 0000003  and  is  fairly  negligible  compared  to  typical  shifts  for 
samples  which  are  at  least  an  order  of  magnitude  higher  (0.0001). 

A  rectangular  bar  shaped  sample  of  alumina  (isotropic  dielectric  properties)  was  placed  in 
the  waveguide  in  position  HE  and  rotated  about  its  long  axis  such  that  at  0  degrees  it  is  at 
position  HE  and  at  90  degrees  at  position  HF.  The  sample  was  rotated  counter  clockwise  in  15 
degree  intervals.  Figure  A. 9  shows  the  results  of  this  test  of  normalized  frequency  shift  verse 
rotation.  The  top  view  of  the  cross  section  of  the  sample  is  illustrated  for  rotation  angles  at  0,  45, 
and  90°  positions. 
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Figure  A. 8  Actual  configuration  of  the  cavity  perturbation  setup  with  sample 
rotation 
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90 


0.001 


Figure  A. 9  Resonance  frequency  shift  (fc-fs)/fs  as  function  of  rotation  angle 

for  a  rectangular  shaped  alumina  ceramic  sample.  Arrows  indicate 
electric  field  directions 

Ideally  the  above  experiment  can  be  used  as  a  means  of  calibration  for  the  NECVP 
method  to  account  for  any  additional  disturbances  due  to  the  insertion  of  rotation  mechanisms, 
the  existence  of  surface  imperfection  of  the  given  cavity,  as  well  as  the  openings  present  on  the 
cavity.  In  the  current  report  however  the  simulation  was  performed  over  a  range  of  permittivities 
at  each  rotation  angle  taking  into  consideration  of  sample  dimension  variations  only,  thus  is  not 
regarded  as  a  full  calibration.  Figure  A.  10  shows  the  derived  permittivity  sr’  of  the  rotated 
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alumina  sample,  corrected  for  dimensional  difference,  that  varies  from  8.71  to  9.44,  adjusted 
based  on  the  standard  deviation  of  the  measured  values. 


Alumina,  rotation,  permittivity 


Angle  (degrees) 


Figure  A.  10  Derived  permittivity  as  a  function  of  rotation  angle  for  a 

rectangular  shaped  alumina  ceramic  bar,  after  correction  of  the 
dimensional  variation 

A.4.4  NEC  VP  Method  for  High  Permittivity  Ferroelectric  Materials 

Microwave  electric  field  strength  inside  a  sample  is  inversely  proportional  to  the 
dielectric  permittivity  of  the  sample  thus  as  permittivity  increases  the  signal  level  and  the 
sensitivity  of  the  measurement  decrease.  Shown  in  Figure  A.l  1  is  a  simulation  of  a  ferroelectric 
relaxor  (l-x)Pb(Mgi/3Nb2/3)03-xPbTi03  (PMN-xPT,  x=0.30)  sample  over  a  range  of  resonance 
frequency  shift  ( fc~f$)/fs ■  The  experimental  results  obtained  gave  a  normalized  frequency  shift  of 
0.00222  that  corresponds  to  a  derived  permittivity  sr’//[iio]  ~250.  The  sample  was  pre-poled  along 
the  [110]  direction.  In  comparison,  sr’//[ooi]  in  the  range  of  400  on  a  PMN-PT  crystal  by  a 
transmission  line  method  at  10  GHz  was  reported  [61]. 

Figure  A.  12  shows  results  for  a  PZT-5H  sample  for  which  experimental  results  gave 
0.00149  for  the  frequency  shift  and  resulting  permittivity  of  about  300.  The  permittivity  of  PZT 
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in  the  microwave  region  has  been  reported  with  values  between  300  and  700  depending  on  the 
composition  [62], 

A.4.5  Effect  of  Meshing  on  NECVP  Measurement  Results 

Size  and  level  of  meshing  during  FEA  process  can  often  be  a  significant  factor  impacting 
resolution  and  thus  accuracy  of  results  from  a  simulation.  The  finite  element  analysis  method  is 
used  as  a  probing  tool  to  depict  the  actual  three-dimensional  EM  field  distribution  in  a  cavity  and 
the  sample  within.  The  quality  of  this  estimation  is  a  direct  result  of  a  proper  meshing.  Typical 
mesh  used  in  this  paper  gave  100,000  degrees  of  freedom  which  is  sufficient  for  most  of  the  low 
permittivity  samples.  However  due  to  the  lower  sensitivity  and  steep  change  in  field  distribution 
for  samples  of  high  permittivity  it  can  often  be  necessary  to  have  up  to  500,000  degrees  of 
freedom.  Figure  A.  13  illustrates  the  importance  of  meshing  on  results  where  we  see  a  different 
converging  curve  for  each  meshing  freedom.  In  COMSOL  free  meshing  parameters  are 
designated  1  through  9  where  1  is  the  finest  mesh.  Mesh  9  has  about  2,000  degrees  of  freedom 
and  mesh  1  about  500,000  degrees  of  freedom.  Time  per  solution  for  mesh  9  is  only  3.5  seconds 
but  20  minutes  are  typically  needed  for  mesh  1 .  The  necessary  meshing  size  can  be  determined 
by  refining  until  further  refinement  yields  little  change  in  the  curve,  i.e.,  achieved  a  prescribed 
accuracy.  More  effective  meshing  can  also  be  achieved  by  using  advanced  or  adaptive  meshing 
strategy  to  have  fine  mesh  in  sub-domain  containing  the  sample,  to  reach  a  balance  between  high 
resolution  and  low  calculation  power  consumption. 
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Figure  A.l  1  Simulation  result  on  the  normalized  shift  of  resonance  frequency 
as  function  of  log  s’r .  The  relative  permittivity  of  the  prepoled 
PMN-30PT  crystal  in  [1 10]  direction  was  derived  by  comparing 
the  measured  value  with  the  simulated  frequency  shift,  using  the 
NECVP  method 
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Figure  A.  12  Simulation  result  on  the  nonnalized  shift  of  resonance  frequency 
as  function  of  log  s’r .  The  relative  permittivity  of  the  PZT-5H 
ceramic  was  derived  by  comparing  the  measured  value  with  the 
simulated  frequency  shift,  using  the  NECYP  method 
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Figure  A.  13  Effect  of  Meshing,  Mesh  1  being  the  finest 


A.5  Summary 

A  numerically  enhanced  cavity  vectorial  perturbation  (NECVP)  microwave  dielectric 
measurement  method  has  been  proposed,  investigated,  and  demonstrated.  In  contrast  to  the 
conventional  microwave  cavity  perturbation  method  that  deals  with  the  shift  of  cavity  resonance 
frequency  under  the  assumption  of  invariant  cavity  EM  field,  this  study  shows  that  three- 
dimensional  EM  field  distribution  resolved  by  using  a  finite  element  analysis  method 
(e.g.,COMSOL),  is  essential  for  an  accurate  evaluation  of  dielectric  materials  that  are  highly 
polarizable,  anisotropic,  and  of  arbitrary  shape.  This  paper  reports  measurement  results  of 
microwave  dielectric  permittivity  by  NECVP  method  of  a  wide  range  of  dielectric  and 
ferroelectric  materials  including  alumina,  glass,  fused  silica,  Teflon,  PMN-PT  crystal  and  PZT 


154 


ceramic  samples.  Microwave  C-band  dielectric  permittivity  ranging  from  2  (for  Teflon)  to  300 
for  PZT  ceramic  are  reported.  This  method  is  shown  to  have  high  resolution  for  the  real  part  of 
dielectric  permittivity  s’  and  good  sensitivity  to  the  imaginary  permittivity  s”.  For  the  samples 
tested  the  results  obtained  compare  well  with  existing  dielectric  properties  reported  using 
conventional  microwave  measurement  techniques,  without  the  restrains  on  specific  form,  shape 
and  volume  of  a  given  sample.  Preliminary  results  were  also  reported  on  the  evaluation  of 
anisotropic  dielectric  properties  by  sample  rotation  within  the  cavity.  Optimization  of  desired 
resolution  and  reduced  calculation  time  may  be  achieved  by  adaptive  meshing  during  the  finite 
element  analysis  process.  The  frequency  range  of  the  reported  NECVP  method  can  be  extended 
and  additional  measurement  parameters  such  as  temperature  and  bias  may  be  included  in  the 
future. 
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